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ABSTRACT 

7 

A  numerical  method  for  solving  the  electrochemical 
transport  equations  subject  to  an  arbitrary  set  of  charge 
transfer  boundary  conditions  is  presented.  The  method  is 
based  on  a  two-dimensional  time-dependent  formulation  of 
the  equations  derived  from  dilute  solution  theory.  It  is 
applied  to  a  study  of  corrosion  and  transport  processes 
in  crack-like  regions.  The  results  show  that  the  average 
species  fluxes  normal  to  the  metal-electrolyte  interface, 
the  crack  length,  and  the  crack  aspect  ratio  are  the  most 
important  factors  affecting  the  electrolyte  composition. 
Other  factors,  such  as  the  crack  shape,  the  form  of  the 
boundary  conditions,  and  the  transport  mode,  also  affect 
the  composition;  but,  these  factors  are  of  secondary  im¬ 
portance.. 
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The  transport  of  ions  in  electrolytic  solutions  plays 
an  important  role  in  every  electrochemical  process.  Just 
how  important  is  determined  hy  several  factors,  but  elec¬ 
trode  geometry  is  the  primary  one.  When  the  geometry  is 
simple,  as  in  the  electroplating  of  continuous  steel  strip, 
ion  transport  is  important  only  at  high  current  densities. 
When  the  electrode  geometry  is  more  complex,  as  in  the  plat¬ 
ing  of  bolts  and  nuts  and  baby  shoes,  the  distribution  of 
the  plating  current  is  strongly  affected  by  the  shape  of  the 
part,  and  transport  in  the  bulk  electrolyte  is  important 
even  at  low  current  densities.  For  decorative  or  ornamental 
plating,  the  uniformity  of  the  current  distribution  may  be 
unimportant;  but  for  precision,  parts  plated  to  obtain  spe¬ 
cific  engineering  properties,  control  of  the  plating  thick¬ 
ness  is  necessary  to  control  both  product  quality  and  pro¬ 
duction  costs. 

In  this  work  we  are  concerned  with  ion  transport  as  it 
affects  the  electrochemical  corrosion  of  metals.  The  geom¬ 
etry  considered  is  a  surface  discontinuity  such  as  a  pit, 
crack,  or  crevice  on  an  otherwise  smooth  electrode  surface. 
The  importance  of  this  geometry  in  metal  corrosion  may  be 
seen  by  inspection  of  any  bridge  in  Pittsburgh.  The  attack 
is  relatively  mild  and  fairly  uniform  along  the  expanse  of 
the  main  structural  members,  and  concentrated  around  the 
joints,  brackets,  and  cover  plates.  The  attack  may  be  par¬ 
ticularly  severe  around  bolts  and  rivets.  In  general,  any 
point  where  two  parts  join,  forming  a  crack  or  crevice,  is  a 
point  of  accelerated  corrosive  attack. 

The  rate  of  attack  in  a  pit  or  crevice  is  often  large 
and  may  be  difficult  to  control.  The  overall  corrosion  cur¬ 
rent  Typically  increases  by  more  than  an  order  of  magnitude 
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at  the  onset  of  pitting  corrosion  [l].  Depending  on  pit 
size  and  distribution,  the  current  density  in  the  pits  them¬ 
selves  may  increase  by  mere  than  two  orders  of  magnitude. 
Crevice  corrosion  may  be  associated  with  equally  large  cur¬ 
rents.  This  is  a  major  nuisance  in  polarization  experiments, 
since  specimens  must  be  mounted  with  extreme  care  to  control 
such  corrosion  and  obtain  uniform  current  densities  [2]. 

Why  surface  discontinuities  such  as  cracks  increase  the 
corrosion  rate  is  poorly  understood.  In  part,  this  is  due 
to  not  knowing  the  electrolyte  composition  in  such  a  region. 
Sampling  the  electrolyte  in  a  crack  poses  some  difficult  ex¬ 
perimental  problems,  so  most  of  the  composition  estimates 
have  been  based  on  indirect  observations.  The  work  reported 
by  B.  F.  Brown,  C.  T.  Fujii  and  S.  P.  Dahlberg  [3]  is  a 
notable  exception.  By  the  rapid  freezing  of  specimens,  they 
obtained  samples  of  the  electrolyte  in  propagating  stress 
corrosion  cracks  tested  in  neutral  salt  solutions.  They 
found  the  pH  to  be  substantially  less  than  that  of  the  bulk 
electrolyte  for  all  of  the  metals  tested. 

While  the  rapid  freezing  technique  is  promising,  the 
small  voiurne  to  be  sampled  constitutes  a  formidable  obstacle 
to  the  experimental  study  of  solution  chemistry  in  crack-like 
regions.  Size  imposes  no  restrictions,  however,  on  a  theo¬ 
retical  stoidy  of  such  problems.  There  are  well-developed 
theories  describing  both  ion  transport  and  electrochemical 
reaction  rates.  The  theories  have  been  successfully  applied 
in  a  variety  cf  electrochemical  transport  problems  [4]  and 
in  the  kinetic  analysis  of  many  charge-transfer  reactions  [?]. 
These  facts  suggest  that  the  theory  might  be  successfully 
applied  to  the  study  of  solution  chemistry  in  crack-like 
regions. 
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The  theoretical  study  of  corrosion  and  transport  in 
crack-like  regions  has  received  little  previous  attend ion. 

J.  G.  Hines  [6]  derived  a  formula  for  the  steady-state 
potential  drop  in  a  wedge-shaped  crack  assuming  constant 
electrolyte  resistivity  and  negligible  transport  by  diffu¬ 
sion.  Somewhat  moi.)  sophisticated  analyses  have  been  pre¬ 
sented  by  D.  A.  Vermilyea  and  C.  S.  Tedmin,  Jr.  [7]  and  by 
W„  D.  France,  Jr.  and  N.  D.  Greene,  Jr.  [8],  However,  both 
of  these  analyses  are  based  on  highly  simplified  models  of 
electrolyte  behavior. 

T.  R.  Beck  [9,  1 0 ]  has  presentc-d  a  model  of  the  elec¬ 
trochemistry  of  stress  corrosion  cracking  in  titanium  alloys. 
Beck's  'MKT'  model  is  the  most  comprehensive  theoretical 
analysis  published  to  date.  The  crack  is  treated  as  a  one- 
dimensional  continuum,  and  the  system  is  assumed  to  be  at 
steady  state.  Convective  transport  is  neglected.  The  elec¬ 
trolyte  contains  three  ions,  two  singly-charged  salt  ions 
and  the  hydrogen  ion.  The  reactions  at  the  metal-electro¬ 
lyte  interface  are  based  on  those  observed  experimentally., 

The  transport  equations  are  solved  numerically  using  a*3 
additional  boundary  data  experimentally  observed  values  of 
current  density  and  potential. 

R.  C.  Alkire,  E.  A.  Grens  II,  and  C.  W.  Tobias  [ll] 
recently  presented  an  analysis  which  is  potentially  more 
powerful  than  the  method  used  by  Beck.  Although  the  method 
was  developed  to  study  porous  electrodes,  it  is  directly 
applicable  to  a  pit-like  region  and  with  simple  modification 
may  be  applied  to  a  crack-like  region.  Like  Beck's,  theirs 
is  a  one-dimensional,  quasi-steady  state  analysis.  In  their 
formulation,  however,  there  is  no  intrinsic  limitation  to 
the  number  of  ions  and  no  boundary  data  other  than  reaction 
kinetics  are  required .  The  results  presented  to  date  are 
fragmentary  and  the  work  is  continuing. 

All  of  the  studies  above  are  based  on  a  simplified 
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formulation  of  the  transport  equations,  hut  they  have  been 
extremely  useful  in  helping  to  understand  corrosion  in 
crack-like  regions.  In  particular,  the  model  developed  by 
Beck  has  been  a  powerful  tool  in  developing  an  understanding 
of  the  role  of  electrochemistry  in  stress  corrosion  cracking. 

The  successful  application  of  electrochemical  transport 
theory  to  the  study  of  electrochemical  phenomena,  including 
those  in  crack-like  regions,  and  the  relative  difficulty  in¬ 
volved  in  the  experimental  study  of  such  phenomena  provided 
the  motivation  for  the  present  work.  The  work  consists  of 
two  parts.  In  the  first,  the  transport  equations  are  pre¬ 
sented  and  the  method  developed  for  their  solution  is  dis¬ 
cussed.  The  method  is  based  on  the  use  of  an  Alternating- 
Direction  Implicit  (ADI)  technique  to  obtain  the  time- 
dependent  solution  to  a  two-dimensional  formulation  of  the 
ion  transport  equations.  The  boundary  data  may  be  of  quite 
general  form.  Concentration,  concentration  gradient,  ion 
flux,  or  charge- transfer  boundary  conditions  may  be  used. 

In  the  second  part,  the  solution  technique  is  employed 
in  a  systematic  analysis  of  the  factors  influencing  solution 
chemistry  and  ion  transport  in  crack-like  regions.  These 
results  are  then  extended  using  a  one-dimensional  formula¬ 
tion  of  the  ion  transport  equations.  Finally,  certain  as¬ 
pects  of  stress  corrosion  cracking  are  discussed  in  the 
perspective  of  the  results  obtained  from  »he  transport 
analysis. 
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II.  TRANSPORT  IN  ELECTROLYTIC  SOLUTIONS 


In  this  section  the  equations  describing  mass  transport 
in  electrolytic  solutions  are  presented.  The  equations  are 
based  on  the  theory  of  dilute  electrolytes.  Several  addi¬ 
tional  assumptions  are  then  introduced;  these  ^re  discussed 
and  their  implications  considered.  Boundary  conditions  are 
then  examined,  and,  because  it  plays  a  key  role  in  the  theory, 
the  electrone\itrality  equation  is  discussed  in  some  detail. 

The  section  concludes  with  a  discussion  of  the  behavior  of 
the  field  equations  and  boundary  conditions  under  the  class 
of  coordinate  transformations  used  in  the  numerical  calcu¬ 
lations. 

The  Field  Equations 

The  difference  between  an  electrolyte  and  an  ordinary 
solution  is  that  the  dissolved  species  may  exist  as  charged 
particles  or  ions  rather  than  as  uncharged  atoms  or  mole¬ 
cules.  The  solvent  itself  is  electrically  neutral  but  has  a 
sufficiently  large  dielectric  constant  to  make  the  existence 
of  the  solute  in  the  form  of  ions  energetically  favorable. 

The  solvent  may  be  a  solid,  liquid,  or  gas.  In  this  work, 
however,  it  is  assumed  that  the  solvent  is  water. 

The  first  element  required  to  describe  mass  transport 
is  a  relation  between  maos  flux  and  the  forces  which  induce 
that  flux*  In  this  work,  the  flux  equation  based  on  the 
theory  of  dilute  electrolytes  (the  Nemst-Planck  equation) 
is  used.  In  developing  this  equation  it  is  assumed  that  the 
total  flux  can  be  expressed  as  the  sum  of  the  fluxes  due  to 
migration  in  an  electrostatic  field,  simple  diffusion,  and. 
solvent  convection.  The  equation  is  written 


=  -  DjVCi  +  c iv 
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where 

p 

J;  =  mass  flux  of  ith  dissolved  species  (mole/cm  -see) 

=  valence  or  charge  number  of  ith  species  (equiv/mole) 

p 

mi  =  mobility  of  ith  species  (mole-cm  /joule-sec) 

3>  =  Faraday’s  constant  (96,520  coul/equiv) 

C i  -  concentration  of  ith  species  (mole/cm^) 

^  =  electrostatic  potential  (volts) 

D;  r.  diffusion  coefficient  of  ith  species  (cm^/sec) 

V  =  solvent  velocity  (cm/sec) 

7  =  'del'  operator  =  *  j +  K 

and  j  ,  and  k_  are  the  unit  vectors  in  the  three  cartesian 
coordinate  directions,  x,  y,  and  z. 

Following  the  usual  convention,  the  gradient  of  the 
electrostatic  potential  is  the  negative  of  the  electric 
field  strength.  The  mobility,  nr  ,  may  be  interpreted  as  the 
average  velocity,  in  meters  per  second,  of  the  ith  dissolved 
species  when  acted  on  by  a  force  of  one  newton  per  mole 
regardless  of  the  origin  of  chat  force.  The  mobility  and 
diffusion  coefficient  are  related  by  the  generalized  Nernst- 
Einstein  equation. 


t>  i 


rwa  + 


dilnOiH 


(2) 


where 

R  =  universal  gas  constant  (8.317  joule/mole-°K) 

T  =  absolute  temperature  (°K) 

#£  =  activity  coefficient  of  .ifch  species  (dimensionless) 
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Each  dissolved  species  also  satisfies  an  equation  of 
continuity  or  mass  conservation. 


cTt  —  -V-Jf  +  G i 


where 

t  =  time  (sec) 

G {  =  rate  of  generation  of  ith  species  (mole/cm  -sec) 

The  quantity  G^  in  equation  (3)  may  he  positive,  nega¬ 
tive,  or  zero.  If  the  ith  dissolved  species  enters  a  homo¬ 
geneous  chemical  reaction  as  a  product,  G^  is  positive;  if 
it  enters  as  a  reactant,  G^  is  negative;  otherwise  it  is 
zero . 

If  a  homogeneous  reaction  is  to  he  included  in  the 
analysis,  a  relation  between  the  rates  of  generation  of  the 
participating  species  and  their  concentrations  must  he  found. 
Consider,  for  example,  a  simple  association-dissociation  re¬ 
action  of  the  form 


SL  +  =£  (4) 

Equation  (4)  states  that  one  mole  of  species  1  combines 
with  one  mole  of  species  M  to  form  one  mole  of  species  N. 

The  rates  of  generation  of  the  participating  species  may  he 
expressed  with  sufficient  accuracy  by  an  equation  of  the  form 


~C\  —  -G'n"  Gh  •-=■  karjtCM  - 
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where 

ka  =  forward  rate  constant  ( 1/mole-cin^-sec) 
k*  =  backward  rate  constant  (l/cm^-sec) 


While  the  rate  constants,  k  and  k. ,  ha/e  the  dimensions 

ad 

indicated,  the  units  are  usually  suppressed  when  there  is  no 

ambiguity  in  their  meaning.  The  forward  and  backward  rate 

constants  are  related  to  the  equilibrium  constant,  K  ,  by 

eq 

the  equation 


G'h  -  o 


K* 

kd 


(6) 


In  this  work  the  solvent  is  assumed  to  be  electrically 
neutral  and  non-conducting.  (The  small  conductivity  of  pure 
water  is  due  to  the  presence  of  hydrogen  and  hydroxyl  ions 
resulting  from  the  equilibrium,  H+  +  OH”  =  ^0.)  Thus,  every 
homogeneous  reaction  must  satisfy  the  requirement  of  charge 
neutrality.  In  the  reaction  described  by  equation  (4),  for 
example,  the  charge  numbers  of  species  L,  M,  and  N  satisfy 
the  equation 


Zu  +  ZM  —  /-N 


(7) 


If,  for  the  moment,  the  convective  term  in  the  flux 
equation  is  set  equal  to  zero ,  one  more  equation  is  required 
to  complete  the  set  of  field  equations.  That  equation  is  a 
relation  between  the  electrostatic  potential  and  the  charge 
density.  The  most  accurate  physical  description  is  provided 
by  Gauss'  law.  For  a  medium  with  a  constant  dielectric 
strength,  this  may  be  written  as 
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M  Ni 

_2*  y' 

e  u 

i*~  i 


44 


(8) 


where 

e  =  dielectric  constant  (coul/volt-cm) 

(7.083  x  10~12  for  H20) 

N«,  =  total  number  of  dissolved  species  in  the  electrolyte 

As  an  alternative,  it  may  be  assumed  that  the  charge 
density  is  everywhere  zero.  This  assumption  is  a  good 
approximation  for  electrolytes  that  are  sufficiently  concen¬ 
trated.  The  condition  of  zero  charge  density  or  electro¬ 
neutrality  is  written 

N» 

Ez,C«  -  0  (9) 

<*>t 


It  should  be  emphasized  thax  either  one  of  equations 

(8)  and  (9)  may  be  used  but  not  both.  The  use  of  equation 

(9)  does  not  imply  that  the  electrostatic  potential  satisfies 

Laplace’s  equation  0).  In  this  work  the  electro¬ 

neutrality  equation  is  employed.  The  relative  merits  of  this 
formulation  are  considered  below. 

To  the  equations  already  discussed  must  be  added  those 
describing  the  motion  of  the  solvent.  These  are  a  contin¬ 
uity  equation,  a  momentum  equation,  and  an  equation  of  state. 
The  continuity  equation  for  the  solvent  is 

'it  +  v-(p^  c  do) 


where 
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fi  =  density  of  the  solvent  (gram/cm^) 

For  a  Newtonian  fluid  with  constant  viscosity,  conser¬ 
vation  of  momentum  is  expressed  by  the  Navier-Stokes  equation 

'  VP  +  x  *  vT~v  (11) 

where 

2 

P  =  hydrostatic  pressure  (dyne/cm  ) 

%  =  body  force  per  unit  volume  (dyne/cm^) 

O 

p  =  coefficient  of  viscosity  (dyne-sec/cm  ) 

^  =  material  derivative  (  ) 

There  is  no  a  priori  constraint  on  the  form  of  the  body 
force,  %  ,  in  equation  (11).  Typically,  however,  one  is 
concerned  with  forces  which  are  expressible  as  the  gradient 
of  a  potential,  for  example,  the  gravitational  force.  It 
is  convenient  to  separate  the  body  force  into  two  components: 
that  arising  from  a  non-zero  charge  density  in  the  solution, 
and  that  due  to  all  other  causes.  Thus,  the  body  force  may 
be  written  as 


where 

4r 

j(,  =  body  force  from  all  sources  other  than  charge 
(dyne/cm^) 

7 

The  factor,  10,  appearing  in  the  equation  results  from 
the  conversion  of  energy  units  from  joules  to  ergs.  Usirg 


equation  (8),  equation  (12)  may  also  be  written 

X  =  io'e[7^]v$  +  X*  (13) 

Examination  of  equations  (11)  and  (13)  reveals  that  the 
motion  of  the  solvent  is  coupled  to  the  motion  of  the  solute 
species  through  the  electrostatic  body  force.  Jhis  coupling 
is  fundamental  to  the  description  of  phenomena  such  as  mass 
transport  through  membranes  and  capillaries. 

The  equation  of  state  is  an  expression  relating  the  sol¬ 
vent  density  to  its  temperature  and  pressure.  For  an  aqueous 
electrolyte  at  room  temperature,  subjected  to  small  gradients 
of  temperature  and  pressure,  the  solvent  may  be  assumed  in¬ 
compressible.  Making  this  assumption,  the  equation  of  state 
is 

p  -  jOo  =  constant  (14) 

The  concept  of  electrical  current  is  closely  related  to 
the  ideas  of  mass  transport  in  electrolytic  solutions.  It 
is  frequently  useful,  and  occasionally  necessary,  to  express 
mass  fluxes  in  terms  of  current  densities.  Since  a  current 
is  simply  a  net  flux  of  charge  and  the  solvent  is  non-con¬ 
ducting,  the  current  density  may  be  written  as 

N, 

(15) 

where 

2 

I  =  current  density  (amp/cm  ) 


Equation  (15)  may  be  rewritten  using  equation  (1)  to 
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express  the  fluxes  and  equation  (9)  to  eliminate  the  co¬ 
efficient  multiplying  the  solvent  velocity.  The  result  is 


The  second  term  on  the  right  is  zero  when  the  composi¬ 
tion  of  the  electrolyte  is  identical  at  every  point  or  when 
the  diffusion  coefficients  are  all  the  same.  In  this  case, 
equation  (16)  reduces  to  Ohm's  law. 


(17) 


where 

K  =  solution  conductivity  (1/ohm-cm) 

From  equations  (3),  (9),  and  (15) »  it 

V-I^  =  0 

Equations  (16)  and  (18)  were  obtained 
neutrality  equation  (eqn.  9).  If  equation 
stead,  the  equations  would  have  a  more  complex  form. 

In  the  arlier  discussion  of  electroneutrality  and 
Gauss'  law,  it  was  observed  that  one  of  the  two  should  be 
used,  but  not  both.  Thus  the  electrostatic  potential  does 
not,  in  general,  satisfy  Laplace's  equation.  To  emphasize 
this  further,  equations  (16)  and  (18)  may  be  used  to  obtain 
the  differential  equation  that  the  electrostatic  potential 
does  satisfy. 


=  (^Cz/rruCi) 

may  be  shown  that 

(18) 

using  the  electro- 
(8)  were  used  in- 
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N4 

[Ez'rrtiCilv^  +  +  vi^jjw]  =  o  (19) 

Equations  (1),  (2),  (3),  (9),  (10),  (11),  (13),  and  (14) 
form  the  complete  set  of  field  equations  describing  mass 
transport  in  electrolytic  solutions.  The  assumptions  made 
in  obtaining  them  do  not  seriously  restrict  their  applica¬ 
tion.  In  using  the  equations  here,  however,  several  addi¬ 
tional  assumptions  have  been  made  and  these  are  examined  be¬ 
low.  In  some  cases  they  have  been  dictated  by  necessity; 
in  others,  they  are  made  for  the  purposes  of  simplicity  only. 

It  has  been  assumed  that  the  solution  is  'ideal*  so  the 
individual  ion  activity  coefficients  are  all  unity.  With 
this  assumption  equation  (2)  reduces  to 


D i  rrifRT 


(20) 


The  assumption  implies  that  the  ion  diffusion  coeffi¬ 
cients  are  concentration  independent.  Also,  it  limits  the 
coupling  between  ions  to  that  arising  from  the  action  of  the 
electrostatic  potential.  The  difference  between  a  diffusion 
coefficient  at  extreme  dilution  and  one  at  one  mole  per  liter 
is  typically  about  twenty  to  thirty  percent.  By  proper 
selection  of  the  concentration  at  which  a  diffusion  coeffi¬ 
cient  is  evaluated,  the  variation  in  a  particular  problem 
can  usually  be  made  much  less  than  this. 

Neglecting  ion-ion  interactions  is  more  restrictive 
than  simply  ignoring  changes  in  the  diffusion  coefficients 
with  concentration.  But,  there  is  very  .little  published 
data  on  the  thermodynamic  properties  of  other  than  single¬ 
salt  electrolytes,  so  there  is  no  practical  basis  for 
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evaluating  the  necessary  activity  coefficients.  If  one 
wanted  to  study  such  ion-ion  interactions,  the  practical 
difficulties  notwithstanding,  it  would  be  preferable  to  use 
a  different  formulation  of  the  flux  equation  (see,  for  ex¬ 
ample,  Reference  12). 

Several  assumptions  have  been  made  to  simplify  the 
equations  describing  the  solvent  velocity  and  to  decouple 
these  equations  from  those  describing  ion  transport.  The 
first  is  that  the  solvent  velocity  is  slowly  time-varying. 
For  an  aqueous  electrolyte,  the  rate  at  which  an  ion  con¬ 
centration  reaches  steady  state  is  two  to  three  crdex-s  of 
magnitude  slower  than  the  rate  at  which  the  solvent  velocity 
reaches  steady  state  following  a  step  change  in  a  displace¬ 
ment  boundary  condition.  Since  the  time  scale  is  determined 
by  the  rate  of  diffusion,  this  assumption  introduces  no 
serious  error. 
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It  is  also  assumed  that  the  acceleration  terms  on  the 
left-hand  side  of  the  Navier-Stokes  equation  are  negligible 
with  respect  to  the  viscous  terms  on  the  right-hand  side. 
Since  the  Reynolds  number  is  less  than  ^0~^  for  all  the 
problems  considered  here,  this  assumption  introduces  no 
serious  error. 

To  decouple  the  equations  governing  the  solvent  motion 
from  those  governing  ion  transport,  it  is  assumed  that  the 
body  force  in  the  Navier-Stokes  equation  is  zero.  The  justi¬ 
fication  for  this  does  not  rest  with  the  use  of  the  electro¬ 
neutrality  equation  as  opposed  to  Gauss'  law.  Even  when 
using  the  electroneutrality  equation,  equation  (13)  should 
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be  used  to  express  the  body  force. 

The  electrostatic  body  force  increases  as  the  current 
density  increases  and  as  the  electrolyte  concentration  de¬ 
creases.  It  is  negligible  with  respect  to  the  viscous  force 
when  the  current  density  is  sufficiently  small  and  the 
electrolyte  concentration  is  sufficiently  large.  It  is 
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difficult  to  establish  a  priori  how  small  is  'sufficiently 
small'  or  how  large  is  'sufficiently  large'.  Equation  (13) 
may  be  used  after  the  fact,  however,  to  estimate  the  magni¬ 
tude  of  the  electrostatic  body  force  using  the  value  of  the 
electrostatic  potential  from  the  solution  of  the  ion  trans¬ 
port  equations.  This  has  been  done  for  the  problems  consid¬ 
ered  in  this  work.  The  electrostatic  body  forces  are  negli¬ 
gible  in  every  case. 

Equation  (11)  may  be  rewritten  using  these  assumptions. 
-  VP  +  /jV‘V  -  0  (21) 


The  final  simplification  is  to  restrict  attention  to 
planar  problems.  If  x,  y,  and  z  are  the  three  orthogonal 
axes  of  a  rectangular  cartesian  coordinate  frame,  and  $  is 
any  function  of  the  independent  coordinates,  this  assumption 
may  be  written  as 


bh  — 

6Z 


0 


(22) 


These  are  the  assumptions  made  in  this  study.  They  are 
typical  of  those  usually  made  and  are  less  restrictive  than 
the  assumptions  invoked  in  most  applications  of  the  theory 
(see,  for  example,  Reference  4).  We  have  discussed  them 
both  to  indicate  the  limitations  of  the  theory  used  here  and 
to  provide  a  more  complete  view  of  the  general  theory  than 
was  presented  above  [ 1 2 ] . 

With  the  above  assumptions,  the  field  equations  may  be 
written  in  more  compact  form.  Using  equations  (,14)  and  (22), 
equation  (10)  becomes 


I'H  - 


0 


(23) 
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Equation  (23)  is  satisfied  identically  using  a  stream 
function,  defined  by 


All) 

'tf  =  curl  ( ^JkJ  = 

-ml 

(24) 

Using  equation  (21)  and  the 

stream  function  defined 

by 

equation  (24),  it  follows  that 

-  0 

(25) 

V2P  =  0 

(26) 

where  P  and  /jV24'  satisfy  the  Cauchy-Riemann  conditions 


Combining  equations  (1)  and  (3)  and  using  equations 
(20)  and  (24),  the  concentration  of  each  dissolved  species 
satisfies  the  equation 


dC  t 
dt 


Z-D.QV^  +  ■?/riVQ*V(j)  +  D;V2Q  -  VG^c.rl(^k)  +  Gy 


(28) 


where 


=  dimensionless  potential 


(rt/2<)  =  25.6  millivolts  at  25°C 


(29) 


The  system  of  field  equations  is  thus  reduced  to 
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equation  (25),  Ns  equations  of  the  form  of  equation  (28), 
and  the  electroneutrality  equation  (eqn.  9). 

Boundary  Conditions 

Consider  a  domain,  0  ,  lying  in  the  x-y  plane  and 
bounded  by  a  curve,  (3.  let  (X  and  _t_  be  the  cartesian 
unit  normal  end  tangential  vectors  to  the  curve  (3,  re¬ 
spectively.  Polio-wing  the  usual  convention,  the  unit  normal 
is  positive  outward;  the  tangent  vector  is  positive  in  the 
counterclockwise  direction.  Within  S3  all  the  necessary 
functions  are  assumed  to  be  defined,  continuous,  and  have 
as  many  continuous  derivatives  as  may  be  required. 

It  is  convenient  here  and  in  later  sections  to  use 
subscript  notation  to  denote  both  vector  components  and 
partial  differentiation.  A  subscript  appearing  alone  in¬ 
dicates  a  vector  component.  A  subscript  following  a  comma 
indicates  differentiation  with  respect  to  that  variable. 

For  example,  if  V,  is  a  vector, 

Vrt  +  Vtt_  (30-a) 

v^  +  v^  -  Vw+vu  (30-b) 

With  this  notation  in  mind,  consider  the  boundary  con¬ 
ditions  for  the  stream  function,  9 .  Within  the  domain,  9 
satisfies  the  biharmonic  equation.  Since  the  equation  is 
fourth  order,  two  constraints  must  be  imposed  at  every  bound¬ 
ary  point.  From  the  definition  of  9  (eqn.  24),  the  most 
natural  boundary  data  are 


arid 


V  =  Vxt+VJ 


+  - 


7*' 

<■  n 
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where 


=  velocity  component  normal  to  the  curve  3 

A 

Vt  =  velocity  component  tangential  to  the  curve  p 

In  practice,  it  is  usually  more  convenient  to  use  the 
value  of  the  stream  function  than  its  tangential  derivative. 

A 

When  %  is  known,  the  stream  function  is  easi.ly  obtained  by 
Integration.  If  s  is  the  coordinate  measured  along  (3, 


r 


=  9(0)  +  \&n&) 


(32) 


where  the  integral  is  taken  in  the  counterclockwise  direc¬ 
tion.  Since  the  fluid  velocity  is  independent  of  the  ab¬ 
solute  value  of  ,  the  integration  constant  ^(o)  may  be 
assigned  any  convenient  value. 

In  this  work  our  primary  concern  is  transport  in  a 
crack-like  region.  In  this  case,  the  boundary  coincides 
over  a  portion  of  its  length  with  a  solid-liquid  phase  bound¬ 
ary  and  the  fluid  motion  is  induced  by  the  motion  of  the 
solid  boundary.  The  velocity  of  the  solid  boundary  is  inde¬ 
pendently  specified  and  it  is  assumed  that  the  fluid  velocity 
is,  at  every  point  on  the  solid-liquid  boundary,  the  same  as 
that  of  the  solid  (the  so-caJied  zero-slip  condition). 

Along  a  line  of  symmetry,  the  tangential  velocity  is 
usually  unknown;  if  so,  it  is  impossible  to  prescribe  the 
normal  derivative  of  the  stream  function.  However,  the  tan¬ 
gential  velocity  typically  has  an  extreme  value,  so  the  spec¬ 
ification  of  the  first  normal  derivative  may  be  replaced  by 
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(33) 
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For  problems  involving  mass  transport  in  a  crack-like 
region  there  is  generally  no  precise  definition  of  the  fluid 
velocity  along  that  portion  of  the  boundary  corresponding  to 
the  mouth  of  the  crack.  Here  it  is  necessary  to  construct 
approximate  boundary  conditions.  This  may  be  done  using  the 
lowest  order  polynomial  that  satisfies  the  differential  equa¬ 
tion  and  is  consistent  with  the  other  boundary  data. 

The  boundary  specification  for  the  ion  transport  equa¬ 
tions  is  more  involved  than  that  for  the  stream  function. 

For  a  solvent  containing  N  dissolved  species,  equations  (9) 

s 

and  (28)  constitute  a  system  of  (N_  +  1 )  equations  in 

s 

(N  +  1)  unknowns.  N  of  these  are  parabolic  partial  dif- 
s  s 

ferential  equations  (first  order  in  time  and  second  order  in 
space)  and  one  is  algebraic.  Suppose  that  N  (N  £  N  ) 

C  C  3 

species  have  non-zero  charge  numbers  while  (N  -  N  )  have 

s  c 

zero  charge  numbers.  The  (N  -  N  )  equations  are  then  inde- 

s  c 

pendent;  their  solution  reduces  to  the  solution  of  (N  -  N  ) 

s  c 

separate  differential  equations. 

The  remaining  equations  are  coupled  through  the  electro¬ 
neutrality  equation.  While  these  have  a  linear  form,  they 
are  in  general  quasi-linear  since  the  electrostatic  potential 
is  itself  a  function  of  the  concentrations. 

Besides  the  coupling  of  the  equations,  the  remaining 
system  of  (N  +  1)  equations  is  remarkable  because  of  the 
character  of  the  coupling.  Since  the  electroneutrality  equa¬ 
tion  is  algebraic,  the  Nc  concentrations  are  linearly  depen¬ 
dent.  This  means  that  exactly  N  [rather  than  (N  +  1)] 

O  v> 

conditions  must  be  imposed  at  each  boundary  point,  and  these 
must  be  such  that  one  is  effectively  a  boundary  condition  on 
the  electrostatic  potential. 

With  the  above  remarks  in  mind,  consider  some  of  the 
forms  which  the  boundary  and  initial  conditions  may  assume. 
The  initial  data  must  be  of  rather  simple  form.  The  concen¬ 
trations  of  (N  -  1 )  of  the  ions  may  be  arbitrarily  specified 

V 
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as  functions  of  the  spatial  coordinates.  The  concentration 
of  the  remaining  ion  is  then  fixed  hy  electroneutrality. 
Because  the  differential  equation  for  the  electrostatic 
potential  exhibits  no  explicit  time  dependence  (see  eqn.  19) > 
the  initial  value  of  the  potential  is  determined  by  the  ini¬ 
tial  values  of  the  ion  concentrations  and  the  boundary  con¬ 
ditions. 

The  boundary  conditions  may  take  any  one  of  a  variety 
of  forms.  However,  all  must  satisfy  the  constraint  imposed 
by  electroneutrality.  Those  easiest  to  use  are  the  ones 
that  reduce  to  one  of  the  classical  forms  for  a  second  order 
partial  differential  equation. 

If  s  is  the  coordinate  measured  along  (5,  Dirichlet 
boundary  conditions  are  of  the  form 
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Q  1  =  C,-M 

i  »i,(Ns-i) 

(34-a) 

.  -i  r  a 
--  ZNi^  2/C^> 

9*  0 

( 34— b) 

|ft  =  fa>) 

(34-c) 

A  A 

where  C^(s)  and  0(s)  are 
nate  s. 

arbitrary  functions  of  the 

coordi- 

Neumann  boundary  conditions  are  of  the  form 

vcrol  =  cfco 

l  =  i,(N*-i) 

(35-a) 

u,,-i 


VCtfnJ  =  \ Chef's!  2*7*0  (35— b) 
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V't-nJ  =  $&) 


(35-c) 


where  C^(s)  and  0  (s)  are  arbitrary  functions  of  the  coor¬ 
dinate  s  along  /3. 

While  Dirichlet  and  Neumann  boundary  conditions  are 
easy  to  apply,  the  normal  derivative  of  0  can  seldom  be 
independently  specified,  and  it  is  usually  necessary  to  use 
the  normal  flux.  From  equations  (1),  (20),  and  (29),  the 
flux  of  the  ith  species  is 


Jf  =  -2A'C4V())  -  tyVCf  +•  AC< 


Equation  (36)  may  also  be  written 


(36) 


k 


$  +  J/ 


where 


£ 

i- 

a 


chemical  flux  = 
convective  flux 


pfifHb  -  D(VC( 
=  +  AC, 


(37) 

(38-a) 

(38-b) 


At  a  point  on  a  moving  boundary,  the  solvent  velocity 
is  equal  to  the  velocity  of  the  solid  by  the  so-called  zero- 
slip  condition.  There  is  generally  no  physical  basis  for 
prescribing  the  normal  component  of  the  total  flux  at  such 
a  point.  The  normal  component  of  the  chemical  flux,  however, 
depends  on  the  rates  of  the  reactions  at  the  solid-electro¬ 
lyte  interface  and  can  usually  be  determined.  Thus,  the  flux 
specification  has  the  form 

Jr  a  =  -(ZAC,70  +  i>vQ).-a  -  J£V),  i=  1.N*  (39) 

C)  . (b 

*  /  \ 

where  the  J^(s)  are  arbitrary  functions  of  the  coordinate  s 


measured  along  the  boundary  curve  (i. 

The  flux  specification  expressed  by  equation  (39)  is 
identical  to  a  Neumann  boundary  condition  for  uncharged 
species.  When  the  normal  flux  is  zero,  the  flux  specif ica- 
tion  looks  like  a  homogeneous  linear  boundary  condition  in 
concentration. 


Charge-Transfer  Boundary  Conditions 

One  of  the  aims  of  this  work  has  been  to  describe  the 
behavior  of  real  physical  systems.  Among  other  things,  this 
requires  the  use  of  boundary  conditions  consistent  with  the 
kinetics  of  electrochemical  reactions.  Boundary  conditions 
of  the  type  described  above  do  not  satisfy  this  requirement. 

In  the  pages  that  follow,  the  relationship  between  mass 
flux  and  reaction  current  density  is  discussed.  The  rela¬ 
tionship  between  the  reaction  current  density  and  the  field 
variables  is  then  examined.  Finally,  the  way  this  infor¬ 
mation  is  combined  to  construct  boundary  conditions  is  de¬ 
scribed. 

To  make  these  ideas  explicit,  consider  as  an  example  a 
simple  redox  reaction  of  the  form 


k£' 


(40) 


where 


k 


reduced  substance 
oxidized  substance 

electrode  reaction  valence  =  (ZQ  -  Zr)  >  0 
=  an  electron  released  into  the  solid  electrode 
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The  reaction  expressed  by  equation  (40)  can  take  place 
only  at  an  electrode-electrolyte  phase  boioidary,  since  the 
electrons  generated  in  the  reaction  must  be  removed  from  the 
electrolyte.  By  convention,  the  equation  is  written  with  the 
electrons  (e~)  on  the  right,  and  the  constant  k  is  a  positive 
integer.  When  the  reaction  goes  to  the  right  it  is  said  to 
proceed  in  the  anodic  direction.  When  it  goes  to  the  left  it 
is  said  to  proceed  in  the  cathodic  direction.  Since  k  is 
positive,  the  charge  number  of  the  oxidized  species  is 
greater  than  that  of  the  reduced  species,  and  an  anodic 
current  corresponds  to  a  flow  of  positive  charge  from  the 
electrode  into  the  electrolyte. 

An  example  of  such  a  reaction  is  the  oxidation  of  the 
ferrous  ion  to  the  ferric  ion  according  to  the  reaction 


Fe++  ^  Fc+++  + 


(41) 


It  should  be  emphasized  that  the  reactions  described  by  equa¬ 
tions  (40)  and  (41)  differ  from  ordinary  chemical  reactions 
in  requiring  the  transport  of  charge  from  one  phase  (the 
electrolyte)  to  another  (the  electrode)  across  a  phase  bound- 


For  reactions  more  complex  than  those  described  by  equa¬ 
tions  (40)  and  (41),  it  is  convenient  to  write  the  reaction 
equation  in  the  more  compact  and  general  form 


4 

C  fyM/  '■  kj£-  =  0 


where 


stoichiometric  coefficient  of  ith  species  in 
jth  reaction 
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M/  =  either  chemical  formula  of  ith  species  or 
molecular  weight  of  ith  species 

kj  =  electrode  reaction  valence  of  jth  reaction 

Tj  =  total  number  of  species  entering  into  jth  reaction 

Note  that  the  stoichiometric  coefficients,  Stj  »  are  less 
than  zero  for  the  reduced  species  in  the  jth  reaction  and 
greater  than  zero  for  the  oxidized  species.  The  stoichio¬ 
metric  coefficients  are  usually  assigned  integer  values. 

The  electrode  reaction  valence  satisfies  the  requirement 


When  the  reaction  described  by  equation  (40)  procedes 
at  an  electrode-electrolyte  interface,  the  generation  of  one 
mole  of  oxidized  species,  SQ,  in  the  electrolyte  is  accom¬ 
panied  by  the  consumption  of  one  mole  of  reduced  species, 

Sr,  and  the  passage  of  k  moles  of  electrons  from  the  elec¬ 
trolyte,  through  the  electrode  and  into  an  external  circuit. 
The  overall  process  may  be  viewed  as  consisting  of  three 
simultaneous  fluxes:  a  flux  of  oxidized  species  into  the 
electrolyte,  an  equal  flux  of  reduced  species  out  of  the 
electrolyte,  and  a  flux,  k  times  greater,  of  positive  charge 
into  the  electrolyte.  The  charge  flux  is  just  an  electrical 
current.  Since  the  mass  fluxes  are  proportional  to  the 
charge  flux,  it  is  convenient  to  use  the  reaction  current  to 
describe  the  reaction  rate.  According  to  the  convention 
above,  the  current  is  positive  (anodic)  when  the  reaction 
procedes  to  the  right. 

Thus,  when  the  boundary  /3  coincides  With  an  electrode¬ 
electrolyte  phase  boundary,  the  fluxes  of  the  species  in  re¬ 
action  (40)  may  be  expressed  as 


^  wWill  Mil 
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Jj- n.  s  —  +  (V  k3>) 
Jo- a  I  =  -(Vka) 


»J0  '  W 


(44-a) 


(44-b) 


where 

p 

I  --  reaction  current  density  (amp/ cm'  ) 

The  fluxes  in  equations  (44)  are  measured  relative  to  the 
boundary,  so  the  convective  components  do  not  appear  in  the 
equations . 

If  more  than  one  reaction  occurs,  the  flux  of  each 
species  is  the  sum  of  the  fluxes  from  all  of  the  reactions. 
Thus,  from  equation  (42), 


111  = 


where 


Ij  =  reaction  current  density  of  jth  reaction 

N;  =  total  number  of  reactions  involving  the  ith  species 


Equation  (45)  is  the  desired  relation  between  mass  flux 
and  the  electrochemical  reaction  rates.  It  has  the  same  form 
as  aquation  (39).  When  the  rates  of  all  the  reactions  are 
constant,  the  boundary  conditions  reduce  to  flux  boundary 
conditions. 

To  complete  the  boundary  specification,  it  is  necessary 
to  express  the  reaction  current  density  as  a  function  of  the 
other  field  variables.  Consider,  as  an  example,  the  reaction 
described  by  equation  (40).  The  reaction  current  density  may 
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be  written 

1  ^OL  ~  Ic  (46) 

where 

1^  =  anodic  current  density  (amp/ cm'  ) 

Ic  =  cathodic  current  density  (amp/cnr) 

In  equation  (46),  the  overall  rate  is  expressed  as  the 
difference  between  the  forward  and  backward  reaction  ratas. 
These  may  in  turn  -..xitten 

la  "  Ka  Ca  <??  p[ 3  k  ~  $*)/' Rl]  ( 47-a ) 

Ic  /cCoC>:p[-(i-cA)k^(}v-^VVT]  (47-b) 

where 

Ko.  =--  forward  (anodic)  r&te  constant  (amp-cm/mole) 

K.  =  backward  (cathodic)  rate  constant  (amp-cm/mole) 

Ca.  =  concentration  of  reduced  species  (mole/enP) 

Co  -  concentration  of  oxidized  species  (mole/cm^) 

a  -  charge- transfer  coefficient  (otc*  :  i  ) 

^wss  electrostatic  potential  of  solid  metal 
electrode  (volt) 

-  electrostatic  potential  in  the  electrolyte  at 
the  phase  boundary  (volt) 

Equations  (47)  are  similar  to  those  describing  the  for¬ 
ward  and  backward  rates  of  a  first-order,  homogeneous  chemi¬ 
cal  reaction.  The  forward  rate  is  proportional  to  the  con¬ 
centration  of  the  reduced  species  (the  reactant);  the 
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backward  rate  is  proportional  to  the  concentration  of  the 
oxidized  species  (the  product).  The  rate  constants,  Ka  and 
K  ,  are  intrinsically  positive.  Although  usually  written  as 
constants,  they  are  actually  functions  of  temperature  and 
pressure.  The  exponential  terms  do  not  appear  in  the  equa¬ 
tions  describing  the  rate  of  a  homogeneous  reaction.  They 
appear  in  equations  (47)  because  the  potential  difference  at 
the  electrode-electrolyte  interface  acts  as  a  barrier  to 
charge  transport.  (For  a  detailed  analysis  of  the  rates  of 
charge-transfer  reactions,  see  Reference  5.) 

A  more  convenient  expression  for  the  reaction  rate  can 
be  obtained  by  rewriting  equations  (46)  and  (47).  When  re¬ 
action  (40)  is  at  equilibrium, 

I  -  Io.  ~  Ic  -  0 


and 


where 

superscript  0  =  value  of  parameter  at  equilibrium 

A<$°=  equilibrium  potential  difference  across 
electrode-electrolyte  interface 

The  rate  constants  in  equation  (48)  may  be  eliminated 

a 

by  introducing  the  standard  electrode  potential,  E  ,  The 
standard  electrode  potential  is  defined  as  the  equilibrium 
potential  difference,  relative  to  the  standard  hydrogen 
electrode,  when  the  activity  of  each  species  entering  the 
reaction  is  one.  It  was  assumed  above  that  the  electrolyte 
is  an  ideal  solution.  Invoking  the  same  assumption  here, 
the  activity  of  a  species  is  one  when  its  concentration  is 
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one  mole  per  liter.  Thus,  from  equation  (48), 


r  -  (SM£) 

+  (i&Mct 


Equation  (50)  shows  how  the  equilibrium  potential  varies 
with  concentration.  The  standard  equilibrium  potential  has 
been  tabulated  for  many  charge-transfer  reactions  [ 1 3 ] -  It 
may  be  determined  directly  from  electrochemical  experiments 
or  by  calculation  from  independently  measured  thermodynamic 
data. 

From  equation  (46)  and  the  condition  for  equilibrium, 
the  anodic  and  cathodic  current  densities  are  equal  at  equi¬ 
librium.  This  value  of  the  current  density  is  called  the  ex¬ 
change  current  density  and  may  be  determined  from  either  of 
equations  (47)  and  equation  (48). 


K.C,exp[aln(^)] 


L  - 


U'OO  a  (i -  o<)  o < 

Ka  Kt  C.  Co 


where 


=  exchange  current  density  (amp/cm“) 


Using  equations  (49),  (50),  and  (51),  the  reaction  rate- 
expressed  by  equations  (46)  and  (47)  may  be  written  in  more 
compact  form. 


I J (t& )c x p f a K o V  R T ]  -  (&)e.vpl'-(i-o)K3«7/'RT] 
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where 

7  =  $»  ~  A$0  (53) 

Equation  (52)  is  the  form  usually  found  in  the  litera¬ 
ture.  The  quantity,  7>  defined  by  equation  (53)  is  called 
the  charge-transfer  overvoltage.  It  indicates  the  degree  to 
which  the  reaction  departs  from  equilibrium.  When  the  con¬ 
centrations  are  fixed  at  their  equilibrium  values,  the  sign 
of  the  current  density  is  the  same  as  the  sign  of  the  charge- 
transfer  overvoltage.  Thus,  the  current  densit;*-  is  anodic 
when  the  overvoltage  is  positive. 

Two  specializations  of  equation  (52)  frequently  appear 
in  the  electrochemical  literature  and  should  be  mentioned. 
When  the  concentrations  of  reacting  species  are  everywhere 
equal  to  their  equilibrium  values  and  the  charge- transfer 
coefficient,  a,  is  equal  to  one-half,  equation  (52)  becomes 

I  —  £I0  sinh  [kc^/Rf]  (54) 

The  other  special  case  occurs  when  the  overvoltage  is 
large.  When  the  argument  of  the  exponential  terms  in  equa¬ 
tion  (52'  is  large  with  respect  to  unity,  one  of  the  two 
terms  is  small  with  respect  to  the  other.  For  example,  when 
7  has  a  large  positive  value,  the  second  term  can  be  ne¬ 
glected  and  the  current  density  is  approximately 

1  **  Io(^)eKp[ak^7/RT]  (55) 

If,  in  addition,  the  concentration  of  the  reduced  species  13 
close  to  its  equilibrium  value,  equation  (55)  may  be  written 


7  A  +  Bloat  I) 


(56) 
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where 

A  =-(^>)Ml.) 
2>  =i^]ln(io) 


( 57— a) 

(57-D) 


Equation  (56)  is  the  so-called  Tafel  polarization  law. 
When  ^  is  plotted  as  a  function  of  log(J),  the  constants  B 
and  A  correspond  to  the  slope  and  zero  intercept  of  the  curve 
respectively.  If  the  electrode  reaction  ■valence v  k,  is 
known,  the  kinetic  constants  &  and  I0,  may  be  determined 
from  equations  (57). 

Per  analytic  work,  the  exchange  current  density  defined 
by  equation  (51)  is  inconvenient  because  its  value  depends 
on  the  choice  of  reference  conditions.  Since  the  standard 
electrode  potential  is  defined  with  unit  activities  of  the 
reacting  species,  it  is  convenient  to  use  the  same  conditions 
to  define  a  standard  exchange  current  density.  Thus, 


I 


$ 

o 


K 


(i-a) 


where 


i. 


Ca  *=  1 

cs*=  i 


(58) 


Io  =  standard  exchange  current  density  (amp-cm/raole) 


The  units  of  the  standard  exchange  current  density  are 
chosen  for  later  dimensional  clarity.  As  before,  it  is 
assumed  in  equation  (58)  that  the  electrolyte  is  an  ideal 
solution. 

Taking  as  reference  conditions  unit  concentrations  of 
the  oxidized  and  reduced  species,  and • introducing  the  dimen¬ 
sionless  potential  from  equation  (29),  equation  (52)  becomes 


(59) 
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Io  {Crtexp[aioJ  -  C0exp[“(!“cA)k?]} 


By  analogy  with  the  dimensionless  potential,  £  will  be 
called  the  dimensionless  overvoltage.  Note  that  a  zero  value 
of  the  dimensionless  overvoltage  . '  -'esponds  to  zero  current 
only  in  the  special  case  when  the  L-oncen  orations  of  .both  the 
oxidized  and  reduced  species  are  equal  to  one  mole  per  liter. 

For  the  example  considered,  the  redox  reaction  described 
by  equation  (40),  equations  (44)  and  (59)  constitute  the  de¬ 
sired  boundary  specification.  Since  the  electrode  is  assumed 
to  be  a  good  conductor,  (J)^  is  not  a  function  of  position  but 
may  be  a  function  of  time.  The  values  of  CA»  C0,  and  (^ap¬ 
pearing  in  equations  (59)  and  (60)  are  those  at  the  electrode 
electrolyte  interface. 

Thus,  when  the  boundary  (3  to  the  domain  <0  coincides 
with  an  electrode-electrolyte  interface,  and  the  redox  reac¬ 
tion  described  by  equation  (40)  is  the  only  charge-transfer 
reaction  occurring  at  that  interface,  the  boundary  specifi¬ 
cation  is 

“  *_  n.  &  =  J^*nj  =  J  (6i-a) 

Jl'fll  =  0  L  =?£  0,/t  (61-b) 

where  ^ 

J  —  (^|cA0xp[o(k^]  -  Co cxp[-a-c*)k5]T  (61 -c) 

Equations  (61)  are  a  straightforward  extension  of  the 
flux  boundary  conditions  discussed  above.  They  are  awkward 
to  handle  in  practice,  since  they  are  non-linear  in  the  di¬ 
mensionless  potential.  But  other  than  convenience,  ther^  is 


iWWWW-iii w 'J*1 
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no  limitation  to  their  use. 

To  generalize  these  results,  note  that  equation  (45)  is 
the  general  relationship  between  mass  flux  and  reaction  cur¬ 
rent  density.  If,  for  each  reaction,  the  current  density  in 
equation  (45)  can  be  expressed  as  a  function  of  concentration 
and  potential,  the  boundary  specification  will  be  complete. 

It  should  be  emphasized  that,  while  equation  (59)  is  typical 
of  those  describing  charge-transfer  kinetics,  equations  of 
this  form  do  not  always  provide  an  adequate  description  of 
reaction  kinetics.  An  example  is  the  reaction  discussed  in 
Appendix  B. 

Electroneutrality 

In  the  original  presentation  of  the  transport  equations, 
it  was  observed  that  Gauss'  law  is  more  accurate  than  the 
electroneutrality  equation.  It  was  asserted  that  the  elec¬ 
troneutrality  equation  is  approximately  satisfied  by  elec¬ 
trolytes  which  are  sufficiently  concentrated,  but  no  evi¬ 
dence  was  presented  in  support  of  that  statement.  The  elec¬ 
troneutrality  equation  is,  nevertheless,  employed  throughout 
this  work. 

At  that  point,  we  had  no  basis  for  selecting  one  equa¬ 
tion  over  the  other.  Here,  having  examined  both  the  field 
equations  and  boundary  conditions,  we  are  in  a  better  posi¬ 
tion  to  discuss  the  merits  of  using  the  electroneutrality 
equation. 

By  using  the  electroneutrality  equation,  the  differen¬ 
tial  equation  for  the  electrostatic  potential  is  simplified 
and  reduced  in  order.  Using  equations  (28)  and  (29)  and 
Gauss'  law  (eqn.  8),  the  differential  equation  for  the  di¬ 
mensionless  potential,  (j) ,  is 
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SjUv^)  -  5Dv44)  -  5v(y14))*cur!(H'k_) 

-  Pv4<t>  -  vP*v$  -  y*Q 


(62) 


where 

D 

P 

Q 


6RT/., 


4Cdl 

E  ZtX>iCi 

£ 

E  El  (Pi -  I))  Cl 


(63-a) 

(63-b) 

(63-c) 

(63-d) 


Using  equations  (28)  and  (29)  and  the  electroneutrality  equa¬ 
tion  (eqn.  9),  the  differential  equation  for  the  dimension¬ 
less  potential  is 


P  V2-^  +  vP*v4>  +  VzQi 


0 


(64) 


Equation  (62)  is  a  fourth-order,  parabolic,  partial 
differential  equation.  Equation  (64)  is  second  order  and 
contains  no  explicit  time  dependence.  Both  equations  look 
linear;  neither  is. 

Rather  than  discuss  equations  (62)  and  (64)  in  their 
general  forms,  we  will  limit  our  attention  to  a  special  case 
for  which  the  forms  are  simplified  but  the  essential  features 
are  retained.  Specifically,  let 
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i) 

D  t  =  D  ,  all  i  =  i,N» 

ii) 

^|=0,  all  f  (steady  state) 

(65) 

iii) 

=  constant  (zero  solvent  velocity) 

Using  these  assumptions,  equations  (62)  and 

(64)  become 

V-(CTV(|»  -  SV4(J)  =  0  (66) 

V-(^TV<t»  =■=  O  (67) 


where 


Ns 


E 


•‘-(A't 


(ionic  strength) 


Using  the  same  assumptions,  equation  (28)  reduces  to 


ZiC;VJ(j)  +  ZtVGt*7(J)  +  V*CL  +  (G,‘/l)J  --■=  o  (68) 

It  has  already  been  shown  that  N  equations  of  the  form 
of  equation  (68)  and  the  electroneutrality  equation  consti¬ 
tute  a  system  of  Ns  equations  in  Ng  unknowns.  On  the  other 

hand,  N  equations  of  the  form  of  equation  (68)  and  Gauss’ 
s 

law  comprise  a  system  of  (N  +  1 )  equations  in  (N  +  1 )  un- 

s  s 

knowns.  The  boundary  specification  which  is  complete  when 
using  the  electroneutrality  equation  is  not  complete  when 
using  Gauss’  law. 

This  result  is  also  reflected  in  equations  (66)  and 
(67).  Equation  (66)  is  fourth-order  in  4>  ;  equation  (67) 
is  second-order.  To  obtain  a  solution  to  equation  (66), 
two  constraints  must  be  imposed  at  each  boundary  point. 
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Equation  (67)  requires  one. 

Boundary  conditions  for  the  system  of  equations  using 
Gauss’  law  may  assume  various  forms.  As  was  done  for  the 
system  based  on  the  electroneutrality  equation,  these  may  be 
determined  from  a  detailed  analysis  of  the  original  system 
of  equations.  The  same  result  may  be  developed  heuristically 
from  equations  (66)  and  (67). 

Jf  the  normal  flux  of  each  species  is  specified  every¬ 
where  on  the  boundary,  the  normal  current  density  is  also 
specified  (by  eqn.  15).  If,  in  addition,  the  ionic  strength, 

,  is  constant ,  flux  boundary  conditions  are  equivalent  to 
specifying  the  normal  derivative  of  $  everywhere  on  the 
boundary. 

A  boundary  condition  of  this  form  is  sufficient  for 
equation  (67).  Another  is  required  for  equation  (66).  One 
possibility  is  suggested  by  Gauss'  law;  namely,  specifying 
the  charge  density  in  the  electrolyte  at  every  point  on  the 
boundary. 

For  the  system  of  equations  including  Gauss'  law,  other 
conditions  could  be  imposed  and  it  is  not  necessary  that  they 
be  the  same  at  every  boundary  point.  Nevertheless,  some¬ 
thing  equivalent  to  the  condition  above  must  be  imposed  in 
addition  to  those  used  for  the  system  based  on  the  electro¬ 
neutrality  equation. 

For  that  portion  of  the  boundary  corresponding  to  an 
electrode-electrolyte  interface,  this  introduces  a  serious 
problem.  While  the  ion  fluxes  can  be  determined  from  the 
rates  of  the  charge- transfer  reactions,  the  first  normal  de¬ 
rivative  of  (£>  is  proportional  to  the  charge  density  on  the 
surface  of  the  electrode.  There  is  published  data  on  the 
rates  of  many  electrochemical  reactions,  but  there  are  vir¬ 
tually  no  data  on  surface  charge  density  for  solid  metal 
electrodes.  Lacking  such  data,  we  would  prefer  not  to  use 
Gauss'  law  if  there  is  some  other  justification  for  using 
the  electroneutrality  equation. 
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Equations  (66)  and  (67)  differ  by  the  appearance  of  a 
biharmonic  term  in  equation  (66)  missing  from  equation  (67). 
To  estimate  the  importance  of  this  term,  consider  the  special 
case  when  the  ionic  strength,  CT  ,  is  constant  throughout  the 
electrolyte.  In  this  case,  equation  (66)  may  be  written 


V*  [$  ~  ^V(j)]  -  O 


where 


&  -  ^rafCzfc' 


The  constant  b  is  called  the  Debye  length.  When  it  is 
small  with  respect  to  the  characterist j c  dimensions  of  the 
domain  £  (the  electrolyte),  the  second  term  in  the  brackets 
in  equation  (69)  is  small  with  respect  to  the  first.  The 
biharmonic  term  in  equation  (66)  can  then  be  neglected,  and 
the  electroneutrality  equation  used  to  describe  the  charge 
distribution  in  the  electrolyte. 

In  Figure  1 ,  the  Debye  length  is  plotted  as  a  function 

_5 

of  ionic  strength.  It  is  less  than  10  centimeters  when 
the  ionic  strength  is  greater  than  2  x  10  equivalents 
squared  per  mole  liter.  For  a  1-1  aqueous  electrolyte  (e.g. 

an  aqueous  sodium  chloride  solution) ,  this  corresponds  to  a 

_c 

salt  concentration  of  10  moles  per  liter.  The  minimum  di¬ 
mensions  of  all  the  domains  considered  in  this  work  are  of 
the  order  of  10  J  to  10  centimeters.  The  ionic  strengths 

_p 

are  typically  of  the  order  of  10  to  1 .0  mole  per  liter. 
Under  these  conditions,  the  electroneutrality  equation  does, 
in  fact,  represent  a  valid  approximation  to  Uauss'  law. 

The  use  of  the  electroneutrality  equation  thus  simpli¬ 
fies  both  the  field  equations  and  the  boundary  conditions. 
Because  it  is  also  an  accurate  approximation  over  a  broad 
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range  of  interest,  it  is  possible  to  exploit  these  simplifi¬ 
cations  , 

Coordinate  Transformations 

When  solving  a  boundary  value  problem,  it  is  often  con¬ 
venient  to  choose  a  coordinate  frame  so  the  boundary  coin¬ 
cides  with  a  coordinate  direction.  For  example,  when  the 
boundary  curve  P  is  a  circle,  it  is  usually  convenient  to 
transform  the  field  equations  and  boundary  conditions  to 
polar  coordinates  so  the  boundary  coincides  with  a  constant 
value  of  the  radial  coordinate.  The  advantage  gained  in 
simplifying  the  boundary  conditions  is  offset  in  part  by 
additional  complexity  in  the  field  equations. 

Under  the  class  of  coordinate  transformations  used  here 
the  field  equations  are  particularly  well-behaved.  Consider 
the  coordinate  systems  (x,y)  and  (£,,?) .  We  wish  to  trans¬ 
form  the  field  equations  and  boundary  conditions  from  x-y 
coordinates  to  £,-?  coordinates  in  such  a  way  that  the  do¬ 
main  in  the  x-y  coordinate  frame  is  transformed  onto  the 
rectangle  $  in  the  £,-7  coordinate  frame.  It  is  assumed  that 
such  a  transformation  exists ,  so  we  can  write 


X  -  XU/7) 
4  = 


(71-a) 
(71 -b) 


y 

.  j  where  the  functions  x  and  y  are  continuous,  single-valued, 

y  and  have  as  many  continuous  derivatives  as  may  be  required. 

In  addition,  it  is  assumed  that  x  and  y  satisfy  the  Cauchy- 
||  Riemann  conditions, 


-38- 


X.,  -  (72-b) 

Under  such  a  transformation,  the  field  equations 
(egns.  28,  9,  25,  ana  26)  assume  the  following  form  ir  the 
€,-7  coordinate  frame: 


y(2Clit  =  t  Dc^Ci  -  VCi-cur^k)  +  Vf&i  (73) 

C^Gi  -  0 

(=i 

f,?m>  +  rM$)  +  F^(f)  +  f4v“p  -  o  os) 

V’P  -  0  (76) 

where 

£  =  unit  vector  in  the  positive  4 -direction  (77-a). 

ni  =  unit  vector  in  the  positive  7-direction  (77-h) 

V  =  +  (77-c) 

*•  -  4  ♦  4  <vm> 

T  2  +  $f<  (77‘e) 

an  -  -  rajf  <77’f> 

H*  =  X.\  +  k'4  (77-g) 

Fi  -  M*  (77'h) 

Ft  -  -4HK,,  (rr-D 

F,  -  -4MK,?  (77‘^ 

Fa  «  <b[K,t  +  Mf7]  -  (77-k) 
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The  operators  in  equations  (77-e)  through  (77-f)  are 
defined  by  analogy  to  the  operators  in  cartesian  coordinates. 
The  function  H  ,  defined  by  equation  (77-g),  is  called  the 
warping  function  and  indicates  the  extent  to  which  a  line- 
element  at  a  point  in  the  x-y  coordinate  frame  is  stretched 
or  compressed  by  the  transformation  at  the  corresponding 
point  in  the  4-7  coordinate  frame. 

The  electroneutrality  equation  (eqn.  74)  is  ■unaffected 
by  the  transformation.  Equation  (76)  is  analogous  to  the 
original  equation.  The  conservation  equation  (eqn.  73) 
differs  from  the  original  equation  only  by  the  appearance  of 
the  warping  function  as  a  coefficient  multiplier  in  the  terms 
not  containing  spatial  derivatives.  Only  the  biharmonic 
equation  (eqn.  75)  is  significantly  affected  by  the  coordi¬ 
nate  transformation. 

The  transformation  of  boundary  conditions  is  straight¬ 
forward.  Any  function  of  the  coordinates  x  and  y  may  be 
written  as  a  function  of  the  coordinates  4  and  7  using  the 
transformation  functions  of  equations  (71).  It  should  be 
observed,  however,  that,  for  an  equation  expressing  a  vector 
quantity  as  a  function  of  the  gradient  of  a  scalar  quantity 
or  the  curl  of  a  vector  potential,  the  warping  function  ap¬ 
pears  in  the  transformed  equation.  Thus,  when  transformed  to 
the  4-7  coordinate  frame,  the  flux  equation  becomes 

i  r  - 

~  “  K  [ZtDjC/Vi])  +•  UVCi  +•  Cicurl(H\kJj  (78) 

Similarly,  equation  (24)  becomes 

(/=■>(  curK^k) 


(79) 
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Equations  (73)  through  (79)  may  be  used  in  any  coordi¬ 
nate  frame,  (4>?)»  for  which  the  transformation  functions 
satisfy  the  Cauchy-Riemann  conditions  (eqns.  72),  The  par¬ 
ticular  form  of  the  transformation  functions  is,  of  course, 
dictated  by  the  shape  of  the  domain  <J3  in  the  x-y  coordinate 
frame . 
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III.  METHODS  FOR  3QLVIN&  THE  TRANSPORT  EQUATIONS 

In  this  section,  some  of  the  methods  for  solving  the 
electrochemical  transport  equations  are  discussed.  The 
numerical  technique  used  in  this  work  is  outlined  ::n  the 
next  section.  Here,  consideration  is  limited  to  the  system 
consisting  of  Ng  equations  of  the  form  of  equation  (28)  and 
the  electroneutrality  equation.  It  is  assumed  that  Ck  =  o 
(for  all  i)  and  that  the  solvent  velocity,  IX,  is  indepen¬ 
dently  specified.  For  convenience,  the  equations  are  dis¬ 
played  below. 

St*  =  ZiDiCiV^  +  HA'VCi-Vtj)  +  DiV'Ci  -  VGi-2/  (eo-a) 

C  Z/C/  =  0  ( 80-b ) 

6 

In  several  cases  the  system  of  equations  reduces  to  a 
system  of  linear  equations.  For  these  cases  solutions  can 
be  obtained  by  classical  methods  [ 1 2 ] . 

When  the  dimensionless  potential  is  cr ' ~ tant ,  equation 
(80-a)  becomes 

ac< 

3t  =  D/V2Q  ”  VC (81  -a) 

Equation  (81 -a)  is  the  so-called  convective-diffusion  equa¬ 
tion.  Since  the  ion  concentrations  are  coupled  through  the 
electroneutrality  equation,  the  solution  set  contains  the 
solutions  of  (N  -  1 )  linearly-independent  differential  equa- 
tions. 

If,  in  addition  to  the  dimensionless  potential,  the 
solvent  velocity  is  zero,  equation  (80-a)  reduces  to  the 
equation  describing  transport,  by  simple  diffusion. 
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Once  again,  the  solutions  are  coupled  through  the  electro¬ 
neutrality  equation,  and  the  solution  set  consists  of 
(N  -  1)  linearly- independent  solutions.  Equations  (81) 

u 

have  been  studied  extensively.  They  can  usually  be  solved 
by  the  method  of  separation  of  variables  or  by  the  use  of 
similarity  transforms.  A  complete  discussion  of  solution 
techniques  with  application  to  specific  boundary  value 
problems  is  contained  in  References  [ 1 4 ]  and  [15]. 

Depending  on  the  solvent  velocity,  one  of  equations 
(81)  can  be  used  when  the  potential  gradient  is  small.  Ne¬ 
glecting  the  potential  introduces  an  error,  but  this  can  be 
tested  by  generating  solutions  to  the  complete  set  of  equa¬ 
tions  and  determining  the  extent  to  which  electroneutrality 
is  violated.  If 
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where  the  C.  are  the  values  obtained  from  one  of  equations 
(81),  the  error  is  ligible. 

When  the  electrolyte  composition  is  everywhere  the  same, 
equations  (80)  reduce  to  Laplace’s  equation  for  the  dimen¬ 
sionless  potential. 

V2(|)  =  0  (83) 

As  is  true  for  the  diffusion  equations,  methods  for  solving 
Laplace's  equation  are  well  documented.  This  equation  has 
been  used  in  a  wide  variety  of  electrochemical  problems.  It 
is  accurate  at  short  times  in  many  time-dependent  problems 
and  at  steady  state  in  well-stjrred  electrolytes. 


Equation  [ 86 )  is  just  the  convective-diffusion  equation 
and  may  be  solved  by  the  technique  discussed  above.  Note 
that  species  2  satisfies  the  same  differential  equation 
(eqn.  86),  and  that  After  solving  equation  (86) 


for  C.j ,  C2  can  be  obtained  directly  from  equation  (84-c). 

The  dimensionless  potential  can  then  be  determined  from 
equation  (84-a).  However,  a  simpler  expression  is  obtained 
by  multiplying  equation  (84-a)  by  and  equation  (84-b)  by 
Z^.  Adding  the  resulting  equations  and  using  equation  (84-c) 
to  eliminate  yields 


C,V2(j>  +  VCi*V(j) 


•  _p  1  ~_p’-  7 


0 


(87) 


While  methods  for  solving  Lap] ace ' s  equation  and  the 
diffusion  equation  are  well  documented,  there  is  no  general 
procedure  for  solving  the  transpor4  equations  (eqns.  80). 

For  problems  involving  dependence  cn  one  spatial  variable, 
the  method  outlined  in  Appendix  C  can  sometimes  be  used. 

The  method  is  based  on  treating  the  potential  $  as  an  un¬ 
known  function  in  the  solution  of  equations  (8C-a).  Using 
these  solutions  and  the  electroneutrality  equation,  an  in¬ 
tegral  equation  is  obtained  for  <$.  The  form  of  the  equation 
varies  from  one  problem  to  another  but  is  typically  non¬ 
linear.  The  usefulness  of  the  method  depends  on  the  effort 
required  to  extract  a  solution  from  the  integral  equation. 

No  way  has  yet  been  found  to  extend  the  method  to  two  or 
three  spatial  variables. 

Approximate  solution  techniques  are  of  two  general 
types:  approximate  analytic  methods  and  numerical  methods. 
Perturbation  methods  are  the  most  powerful  of  the  analytic 
ones.  These  are  based  on  the  idea  that  a  small  change  in  one 
of  the  field  variables  or  boundary  conditions  should  cause 
only  small  changes  in  the  solution.  The  equations  for  the 
binary  electrolyte  (eqns.  86  and  87)  provide  the  most  satis¬ 
factory  basis  for  such  methods. 

In  the  simplest  perturbation  method,  the  concentrations 
of  the  two  principal  ions  are  determined  from  equation  (86). 
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The  concentrations  of  the  minor  species  are  determined  from 
equation  (28),  using  the  electrostatic  potential  obtaiiiCd 
from  equation  (87).  The  electroneutrality  equation  is  then 
used  to  obtain  an  improved  estimate  of  the  potential.  There 
are,  of  course,  more  sophisticated  methods  and  they  may  be 
used  iteratively  to  improve  the  accuracy  of  the  approxima¬ 
tions. 

Numerical  methods  have  received  increasing  attention  in 
recent  years.  Much  of  the  work  has  been  based  on  some  spe¬ 
cialization  of  the  ion  transport  equations.  For  example, 
the  work  reported  in  Reference  [ 1 6 ]  is  based  on  the  convec¬ 
tive-diffusion  equation.  Such  approximations  can  be  extreme¬ 
ly  useful,  but  only  that  work  including  the  effects  of  both 
diffusion  and  migration  is  considered  here. 

The  'MKT1  analysis  dev  loped  by  T.  R.  Beck  [ 1 0 ]  is 
basically  a  one-dimensional,  steady-state  formulation  of  the 
ion  transport  equations  using  a  particular  set  of  charge- 
transfer  boundary  conditions.  It  is  assumed  that  the  solvent 
velocity  is  zero  and  that  the  electrolyte  contains  three 
charged  species.  The  equations  are  integrated  numerically 
using  a  Runge-Kutta  technique,  and  the  boundary  conditions 
are  satisfied  by  trial  and  error. 

R.  C.  Alkire,  E.  A.  Grens  II,  and  C.  W.  Tobias  [ll] 
recently  presented  a  more  sophisticated  treatment  of  one¬ 
dimensional  transport.  Developed  for  a  study  of  porous 
electrodes,  their  analysis  is  based  on  a  quasi-steady  state 
formulation  of  the  transport  equations  in  which  the  solvent 
velocity  is  determined  by  the  rate  of  dissolution  of  the 
metal  electrode.  The  equations  are  linearized  and  then  cast 
in  finite  difference  form.  The  resulting  system  of  algebraic 
equations  is  solved  numerically.  An  iterative  procedure  is 
used  to  simultaneously  satisfy  both  the  field  equations  and 
boundary  conditions. 

There  has  been  no  report  of  methods  for  problems  in  two 
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or  three  spatial  variables.  The  methods  sported  in  Refer¬ 
ences  [10]  and  [ 1 1 ]  could  be  generalized  to  electrolytes 
containing  more  than  three  charged  species  and,  without 
much  difficulty,  to  include  explicit  time-dependence.  It 
is  doubtful,  however,  that  the  method  used  in  Reference  [ 1 0 ] 
can  be  extended  to  problems  in  two  or  three  spatial  vari¬ 
ables. 


IV.  COMPUTER  PROGRAM  DEVELOPED  FOR  THIS  STUDY 


Several  specific  goals  were  established  to  guide  the 
design  of  the  computer  program.  The  primary  goal  was  that 
the  program  be  capable  of  generating  time-dependent  solu¬ 
tions  to  the  electrochemical  transport  equations,  using  an 
arbitrary  set  of  charge- transfer  boundary  conditions. 
Secondly,  it  was  desired  that  the  program  be  flexible  and 
have  a  broad  applicability.  While  the  immediate  aim  of  this 
study  was  to  examine  transport  in  crack-like  regions,  the 
program  is  designed  to  solve  problems  in  two  spatial  vari¬ 
ables  rather  than  one.  Also,  the  program  can  be  easily  spe¬ 
cialized  to  solv*'  Laplace's  equation  or  the  convective- 
diffusion  equation  by  the  input  of  appropriate  control  char¬ 
acters. 

Basically,  the  program  calculates  three  functions:  the 
stream  function,  the  dimensionless  electrostatic  potential, 
and  the  concentration  of  each  dissolved  species.  In  every 
case  the  approach  is  the  same.  The  controlling  differential 
equation  is  replaced  by  an  equivalent  difference  equation  at 
a  number  of  preselected  points.  The  resulting  system  of 
linear  algebraic  equations  is  then  solved  by  Gaussian  elim¬ 
ination.  Additional  subroutines  are  required  to  perform 
service  an-3,  control  functions.  The  most  important  of  these 
are  the  subroutines  which  establish  the  size  and  shape  of  the 
physical  region  and  those  in  which  boundary  conditions  are 
calculated  and  applied. 

The  problem  to  be  solved  is  defined  by  the  inpui  data 
under  the  control  of  the  executive  portion  of  the  program. 
First,  the  initial  concentrations  and  the  transport  proper¬ 
ties  of  the  dissolved  species  are  specified.  Up  to  ten 
species  may  be  included  in  a  single  problem.  The  type  of 
problem  to  be  run  (simple  diffusion,  convectj.ve  diffusion, 
or  electrochemical  transport)  is  then  determined  and  the 
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form  of  the  boundary  conditions  is  established.  Finally, 
the  size  and  shape  of  the  physical  domain  are  fixed. 

The  physical  domain  may  be  hyperbolic,  elliptic,  pie¬ 
shaped,  or  rectangular.  Other  shapes  are  possible  but,  as 
the  program  is  now  written,  the  transformation  functions  are 
limited  to  those  listed.  The  way  the  domains  transform  to  a 
rectangle  is  indicated  in  PiDure  2.  The  numbered  points 
transform  to  the  numbered  corners  of  the  rectangle,  and  the 
lines  connecting  them  transform  to  straight  lines.  The  size 
of  the  physical  domain  is  determined  by  two  characteristic 
dimensions,  the  length  and  width.  The  way  these  are  measured 
is  also  t'hown  in  Figure  2. 

One  finite  difference  grid  is  used  for  all  the  calcula¬ 
tions.  The  grid  spacing  is  uniform  in  each  of  the  directions 
£,  and  7*  The  maximum  grid  size  which  may  be  used  is  11  by 
41.  As  an  example,  the  grid  for  a  hyperbolic  domain  is  shown 
in  Figure  3  in  both  the  physical  coordinates  and  the  trans¬ 
formed  coordinates.  Note  that  the  spacings  in  the  and 
7-directions  are  uniform  but  not  equal. 

All  quantities  related  to  the  motion  "f  the  solvent  are 
derived  from  the  stream  function.  The  stream  function  it¬ 
self  is  calculated  using  a  straightforward  finite-difference 
analogue  to  equation  ('75).  The  resulting  system  of  alge¬ 
braic  equations  is  solved  by  Gaussian  elimination  with  back 
substitution  to  provide  error  control.  The  solvent  velocity 
at  each  grid  point  is  determined  from  the  finite-difference 
analogue  to  equation  (79). 

The  pressure  at  any  point  (£,0,7o)  in  'the  solvent  can  be 
determined  by  integration  of  equation  (21).  If  is  any  con¬ 
venient  curve  connecting  the  two  points  (0,0)  and  (£,0)70) 
then 
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The  pressure  is  calculated  at  points  along  the  crack  center- 
line  according  to  equation  (88).  The  appropriate  derivatives 
are  replaced  by  their  finite-difference  approximations  and 
the  integration  is  performed  numerically  using  the  trape¬ 
zoidal  rule. 

From  equations  (73)  and  (74),  the  differential  equation 
for  the  dimensionless  potential,  (j>,  in  terms  of  the  trans¬ 
formed  (%,?  )  coordinates  is 
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The  potential  is  calculated  using  a  straightforward  finite- 
difference  analogue  to  equation  (89)  at  each  grid  point. 

The  resulting  system  of  algebraic  equations  is  solved  by 
Gaussian  elimination. 

Since  there  is  no  explicit  time  dependence  in  equation 
(89) ,  the  potential  is  established  instantaneously.  For 
problems  in  electrochemical  transport,  this  initial  value 
must  be  calculated  before  attempting  to  calculate  ion  con¬ 
centrations  at  an  advanced  time.  When  the  electrolyte  is 
initially  homogeneous,  equation  (89)  reduces  to  Laplace’s 
equation. 

When  charge- transfer  boundary  conditions  are  used,  the 
initial  estimate  of  the  boundary  conditions  is  generally  in 
error.  In  this  case  it  is  necessary  to  repeat  the  calcula¬ 
tion  of  the  potential,  using  the  previous  value  to  obtain  an 
improved  estimate  of  the  potential  difference  across  the 
electrode-electrolyte  interface.  Four  to  six  such  iterations 
are  usually  sufficient  to  reduce  the  error  everywhere  to  less 
than  one  part  in  ten- thousand. 

At  advanced  times,  the  initial  estimates  of  concentra¬ 
tions  and  potential  all  contain  errors  regardless  of  the 
form  of  the  boundary  conditions.  Then  it  is  necessary  to 
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employ  an  iterative  procedure  involving  all  the  field  equa¬ 
tions  and,  in  some  cases,  the  boundary  conditions.  This 
procedure  is  discussed  more  fully  after  first  discussing  the 
method  used  to  advance  the  calculation  of  concentration  in 
time. 

Computation  of  the  ion  concentrations  at  an.  advanced 
point  in  time  is  carried  out  in  two  stages.  The  concentra.- 
tion  of  each  species  is  calculated  using  a  finite-difference 
analogue  to  equation  (73).  The  electroneutrality  equation 
(eqn.  74)  is  then  used  to  modify  these  values  and  obtain  im¬ 
proved  concentration  estimates. 

Initial  estimates  of  the  ion  concentrations  at  a  new 
point  in  time  are  obtained  using  an  ADI  (Alternating-Direc¬ 
tion  Implicit)  technique.  A  technique  of  this  form  was 
selected  because,  unlike  forward-difference  analogues  to 
equation  (73),  such  techniques  are  unconditionally  stable. 
There  is  no  limit  to  the  size  of  the  time  step  beyond  the 
practical  limit  required  to  control  truncation  error.  In 
addition,  ADI  techniques  are  relatively  fast  as  compared  to 
other  stable  finite-difference  approximations,  e.g.,  back- 
ward-difference  methods. 

When  using  an  ADI  technique  with  an  equation  in  two  or 
three  spatial  variables,  the  problem  is  decomposed  into  two 
or  three  problems,  each  one  of  which  exhibits  implicit  de¬ 
pendence  on  only  one  spatial  variable.  In  the  present  case 
the  differential  equation  depends  on  the  spatial  variables, 
and  7.  The  finite-difference  analogue  to  this  equation 
consists  of  two  parabolic  difference  equations.  One  carries 
an  implicit  dependence  in  the  ^-direction,  the  other  in  the 
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^-direction.  In  operation,  the  equations  for  the  ^-direction 
are  solved  first;  then  the  equations  for  the  7-direction  are 
solved.  This  procedure  ie  repeated  at  succeeding  time  steps; 
hence  the  name  Alternating-Direction  Implicit.  At  each  stage 
the  coefficient  matrix  for  the  system  of  algebraic  equations 
is  tri-diagonal. 

A  more  detailed  discussion  of  the  finite-difference 
analogue  used  in  the  solution  of  equation  (73)  may  bn  found 
in  Appendix  A.  A  comprehensive  review  of  numerical  methods 
for  the  solution  cf  parabolic  differential  equations,  in¬ 
cluding  ADI  techniques,  is  presented  in  Reference  [l7j. 

In  earlier  versions  of  the  program  the  elecironeutral- 
ity  equation  was  used  only  to  generate  an  equation  arid  bound¬ 
ary  conditions  for  the  electrostatic  potential.  It  was  not 
used  explicitly  in  the  program.  Although  several  different 
methods  were  tried,  none  were  satisfactory:  all  resulted  in 
systematic  but  uncontrolled  oscillations  in  the  electro¬ 
static  potential.  While  the  behavior  differed  from  one 
method  to  another,  the  oscillations  were  always  accompanied 
by  oscillations  in  the  ion  concentrations  and  systematic 
deviations  from  electroneutrality. 

In  the  present  program,  oscillations  of  the  icn  concen¬ 
trations  and  electrostatic  potential  have  been  brought  under 
control  by  explicit  use  of  the  electroneutrality  equation. 

The  method  used  is  redundant.  The  concentrations  of  indi¬ 
vidual  species  are  initially  treated  as  being  independent, 
and  earn  concentration  is  calculated  using  the  ADI  technique 
discussed  above.  This  yields  one  concentration  estimate.  A 
second  estimate  is  obtained  using  the  electroneutrality  equa¬ 
tion  and  the  first  concentration  estimates  for  all  the  other 
charged  species.  These  estimates  are  then  combined  to  yield 
an  improved  estimate. 

In  a  typical  problem,  the  concentration  cf  a  particular 
ion  may  be  several  orders  of  magnitude  less  than  the  total 
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ion  concentration.  For  this  reason  it  is  necessary  to  weight 
the  two  concentration  estimates  rather  than  take  a  simple 
average.  For  each  ion  a  weighting  function,  g . ,  is  computed 
according  to  the  formula 


.  fa '*i(V(££ 


£-1  m=i 


Cm) 


where 

Oj(  =  weighting  function  for  ith  charged  species 

C?K  =  estimate  of  the  concentration  of  the  ith  species 
at  the  kth  grid  point  obtained  using  the  ADI 
technique 

Np  =  total  number  of  grid  points 
Nc  =  total  number  of  charged  species 

The  concentration  at  each  point  is  then  determined 
according  to  the  formula 


Civ.  =  (l  -  +  tyCjk 


where 


'best'  estimate  of  concentration  of  ith  species 
at  the  kth  grid  point 


N* 

E  ZjC'Jk 
j-s 

,  estimate  of  the  concentration 
of  the  ith  species  at  the  kth 
grid  point  based  on  the  elec¬ 
troneutrality  equation. 

The  precision  of  this  method  is  quite  good.  The  error, 
,  in  satisfying  the  electroneutrality  equation  may  be  de¬ 
fined  as 
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that  is,  the  maximum  value  of  the  ratio  at  any  grid  point 
within  the  electrolyte.  On  a  computer  carrying  the  equiva¬ 
lent  of  eight  decimal  digits  in  single  precision  arithmetic, 

_7 

no  value  of  CEN  greater  than  5  x  10  has  been  observed. 

This  is  the  same  magnitude  as  the  roundoff  error  in  the  cal¬ 
culation  of  etv) . 

When  the  concentrations,  Ci.(t0),  and  the  dimensionless 
potential,  (j>(t0)  are  specified  at  a  point  in  time,  t0  ,  the 
calculations  at  a  new  point  in  time,  ),  are  carried 

out  as  follows.  First,  estimates  of  the  concentrations, 

Ci(  t0+  At- ) ,  are  calculated  using  the  value  of  the  potential, 
<j>(to) ,  to  approximate  its  value  at  the  intermediate  time, 
(t0+£At).  These  estimates,  Ci(  t0t  ht ) ,  are  then  used  to 
estimate  the  potential,  (f>(  t-,+  At.  ) ,  at  the  new  time.  When 
necessary,  as  when  charge-transfer  boundary  conditions  are 
used,  the  estimates  of  quantities  at  ( t*+  )  are  then  used 

to  revise  the  boundary  conditions.  This  sequence  is  repeated 
three  or  four  times  as  may  be  required  to  reduce  the  error  to 
the  desired  level.  In  each  such  sequence,  the  latest  esti¬ 
mate  of  the  potential  at  the  new  time  and  the  value  at  the 
previous  time  are  used  to  estimate  (^(i0t  j  ,\t.  )  in  the  compu¬ 
tation  of  the  concentrations ,  ).  If  more  than  four 

such  iterations  are  required  it  is  usually  preferable  to  re¬ 
duce  the  size  of  the  time  step. 

In  Appendix  A,  it  is  shown  that  the  ADI  method  is  second 
order  correct  in  both  time  and  space  when  the  dimensionless 
potential  is  evaluated  at  the  intermediate  point  in  time, 
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(  to+^At ).  Because  cj>  is  not  specified  a  priori,  iteration 
is  required  to  retain  this  second-order  r.*rrectness.  The 
iterative  procedure  is  also  effective  in  controlling  other 
errors  of  uncertain  magnitude.  When,  for  example,  charge- 
transfer  boundary  conditions  are  used,  the  boundary  condi¬ 
tions  at  the  new  time  are  not  known  at  the  beginning  of  the 
time  step.  Unless  the  solution  technique  leads  to  simul¬ 
taneous  convergence  of  both  the  field  variables  and  the 
boundary  conditions,  the  theoretical  accuracy  achievable 
using  the  ADI  method  cannot  be  realized. 

Boundary  conditions  for  the  transport  equations  are 
imposed  by  straightforward  application  of  the  principles 
discussed  in  Section  II.  When  charge- transfer  boundary 
conditions  are  used,  however,  the  procedure  is  somewhat 
involved  and  this  case  deserves  some  explanation. 

For  each  charge-transfer  reaction  a  kinetic  equation, 
such  as  equation  (59),  must  be  obtained  and  includeu  in  the 
program.  Storage  capacity  has  been  set  aside  to  permit  the 
use  of  up  to  nine  sucn  equations.  The  current  density  is 
calculated  for  each  reaction  and  for  as  many  boundary  points 
as  may  be  desired.  The  calculation  is  always  based  on  the 
best  estimates  of  the  concentrations,  CL(t0+At  ),  and  poten¬ 
tial,  <])(  t»0+  At  )  at  the  advanced  point  in  time  available  at 
the  tine  of  the  calculation. 

The  potential  of  the  metal  (electrode)  may  be  indepen¬ 
dently  specified.  It  may  be  held  constant  or  permitted  to 
vary  with  time.  As  the  program  is  now  written,  however,  it 
is  not  possible  to  impose  a  constant  current  condition  on 
the  electrode  (corresponding,  for  example,  to  a  freely 
corroding  metal). 

The  next  step  is  to  compute  the  normal  flux  for  each 
species  at  each  boundary  point.  The  computation  is  performed 
according  to  equation  (45)  using  tne  values  of  the  current 
density  just  calculated. 


Finally,  the  normal  derivative  of  the  dimensionless 
potential,  (j>iru  is  calculated  at  each  boundary  point  accord¬ 
ing  to  the  equation 

4^  =  ~  U  (^^TOT  +  ^>n.  )  (93) 

where 

(j),a  =  normal  derivative  of  4)  at  a  point  on  the  boundary 

^TOT  =  J.*.^ 

U  =  C  zf  DiCi 
V  =  C  ZiDiCi 

and  where  the  quantities  I Tcr  ,  U ,  and  V  are  all  evaluated 
at  points  on  the  boundary. 

The  current  density,  Ito1  ,  is  taken  as  the  sum  of  the 
partial  current  densities  calculated  earlier.  The  calcula¬ 
tion  is  always  based  on  the  best  estimates  of  the  concentra¬ 
tions,  Cj.(  ta+  At )  at  the  advanced  point  in  time. 

The  estimates  of  the  normal  fluxes  and  the  normal  de¬ 
rivative  of  the  dimensionless  potential  are  used  with  the 
finite-difference  analogue  to  the  flux  equation  to  construct 
the  boundary  conditions  for  each  species.  As  mentioned  a- 
bove,  all  these  estimates  contain  errors  depending  on  the 
stage  of  the  computation.  This  method  of  constructing 
boundary  conditions  has  other  drawbacks  which  are  discussed 
at  the  close  of  this  section. 

The  printed  output  includes  coordinate  data,  data  rele¬ 
vant  to  the  solvent  velocity,  and  transport  data.  The  (£,-?) 
and  (x-y)  coordinates  are  calculated  and  printed  for  each 
grid  point.  The  stream  function  and  its  first  derivatives 
in  the  and  7 -directions  are  also  printed.  The  value  of 
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the  solvent  velocity  at  each,  grid  point  is  presenter,  in  two 
formats.  In  one  the  x-  and  y-components  of  the  velocity  are 
printed;  in  the  other  the  magnitude  and  direction  of  action 
of  the  velocity  vector  are  presented.  The  transport  data  are 
printed  at  time  zero  and  at  the  end  of  each  time  step.  These 
data  include  the  value  of  the  dimensionless  potential  and  the 
values  of  the  concentrations  at  each  grid  point.  The  partial 
current  density  for  3ach  reaction  is  printed  for  every  bound¬ 
ary  point  lying  on  the  solid-electrolyte  phase  boundary. 
Finally,  the  x-  and  y-components  of  the  total  current  density 
are  printed. 

The  program  has  been  run  on  a  Univac  Corporation  Model 
1108  computer.  The  running  time  depends  on  the  number  of 
dissolved  species,  the  mode  of  transport,  and  the  number  of 
time  steps  desired.  For  a  simple  diffusion  problem  with  one 
dissolved  species  and  charge-transfer  boundary  conditions, 
two  to  three  seconds  of  computation  time  are  required  for 
each  time  step  when  using  four  iterations  per  time  step. 

For  an  electrochemical  transport  problem  with  three  dissolved 
species  and  constant  flux  boundary  conditions,  the  computa¬ 
tion  time  is  about  fifteen  seconds  per  time  step,  again  using 
four  iterations  per  time  step.  More  than  half  this  time  is 
required  for  the  computation  of  the  electrostatic  potential. 
The  total  running  time  for  a  wide  variety  of  problems  is 
typically  between  five  and  thirty  minutes. 

Although  the  program  has  been  successfully  employed  in 
the  solution  of  a  variety  of  transport  problems,  it  has  two 
important  limitations.  The  method  used  to  construct  and 
apply  boundary  conditions  limits  the  capability  of  the  pro¬ 
gram.  In  particular,  difficulties  are  encountered  in  prob¬ 
lems  involving  either  diffusion-limited  partial  current  den¬ 
sities  or  extremely  large  total  current  densities.  These 
limitations  may  be  overcome  by  changes  in  the  subroutine 
used  to  compute  the  partial  current  densities.  That  is, 
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rather  than  treating  the  partial  current  density  as  a  con-  j 

stant,  it  might  he  treated  as  a  linear  function  of  both  con-  j 

centration  and  potential,  where  the  linear  approximation  is 
made  about  the  current  values  of  the  field  variables. 

Another  limitation  derives  from  the  use  of  the  current 
equation  to  compute  c{Vi,  the  normal  derivative  of  the  dimen¬ 
sionless  potential.  This  results  in  a  coupling  between 
errors  in  the  estimates  of  cj>)n.  and  the  concentrations.  In 
electrochemical  transport  problems  this  introduces  an  insta¬ 
bility  which  can  be  controlled  only  by  controlling  the  size 
of  the  time  step.  While  this  method  is  successful,  it  is 
undesirable  because  it  limits  full  exploitation  of  the  ADI 
technique.  It  is  felt  that  thi-T  limitation  can  be  overcome 
by  using  the  sum,  CZiJi/Di  ,  rather  than  the  current  equa¬ 
tion,  for  calculating  <j)irx. 
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V.  RESULTS 


Ine  resets  presented  in  this  section  form  the  major 
part  of  an  effort  to  develop  a  systematic  and  coherent  view 
of  the  factors  affecting  solution  chemistry  and  ion  trans¬ 
port  in  crack-like  regions.  The  results  are  limited  to  the 
behavior  of  an  electrolyte  in  a  pre-existent  notch  or  crack 
on  an  otherwise  smooth  electrode  surface.  For  the  sake  of 
clarity  and  simplicity,  attention  is  focused  on  a  single 
reaction  and  it  is  assumed  that  no  homogeneous  chemical  re¬ 
action  occurs  in  the  electrolyte. 

The  reaction  chosen  for  study  is  an  oxide-film  forma¬ 
tion  reaction.  Such  reactions  play  an  important  role  in 
virtually  all  forms  of  metal  corrosion.  The  rate  of  this 
kind  of  reaction  cannot  be  adequately  described  by  a  kinetic 
equation  of  the  form  of  equation  (59)*  Under  suitable  con¬ 
straints,  however,  behavior  similar  to  that  expressed  by 
equation  (59)  nay  be  observed  (see  Appendix  B) . 

It  is  assumed  that  the  oxide  is  formed  by  the  reaction 

M  +  Z  HzO  =t=  M02  +  4H+  4-  4<Z'  (94) 


For  consistency  and  completeness  it  is  necessary  to  include 
oxide  formation  by  direct  reaction  with  the  electrolyte. 
Such  a  reaction  is  represented  by 


MO 4+  +  HzO  MO*  +  £HX 


(95) 


This  reaction  does  not  require  the  transport  of  electrons 
through  the  oxide  film.  The  reaction  procedes  to  the  right 
or  left  depending  only  on  the  concentrations  of  the  dissolved 


species.  The  reaction  rate  is  independent  of  the  electrode 
potential. 

The  metal  oxide  is  assumed  to  have  zero  electronic  con¬ 
ductivity.  All  current  passing  through  the  film  is  conducted 
by  the  migration  of  ions  or  vacancies.  As  a  consequence, 
no  redox  reaction  can  occur  in  the  system  while  the  oxide  • 
film  is  intact.  Thus,  equations  (94)  and  (95)  represent  the 
complete  set  of  possible  reactions  occurring  at  the  oxide- 
covered  metal  surface.  This  assumption  provides  an  accurate 
physical  description  of  the  oxides  of  titanium  and  aluminum. 

The  rate  of  reaction  (94)  can  be  expressed  as 

Ii  =  I,o{exp[^]  -  exp  [~£i]}  (96) 

where 

0  r  , 

-  b  L<fv-  Ti  -  e,  -  ?,]  (97) 


and  where 


2 

I io  =  exchange  current  density  ^  amp/ cm' ) 

(3  =  high  field  conduction  coefficient  (cm) 

&  =  oxide  film  thickness  (cm) 

4>m,  =  dimensionless  potential  of  metal  electrode 

^  =  dimensionless  potential  of  electrolyte  at 

metal-electrolyte  interface 

CF  =  Plade  potential  (dimensionless  equilibrium 
potential  of  reaction  94) 

=  surface  overpotential 
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The  Flade  potential,  £r ,  is  taken  as  the  equilibrium 
potential  of  reaction  (94).  Thus,  eF  is  a  function  of  the 
hydro son  ion  concentration.  It  is  easily  shown  that 


e, 


e; 


+ 


i0+3  [Hf]f 


(98) 


where 

e?  =  standard  dimensionless  equilibrium  potential 
of  reaction  (94) 

K  =  hydrogen  ion  concentration  (mole/cm^) 


The  surface  overpotential,  %  ,  is  determined  by  the  rate 
of  reaction  (95).  It  may  be  determined  from  the  equation 


(99) 


where 

=  exchange  current  density  for  anodic  partial 

p 

electrode  reaction  for  reaction  95  (amp/cm  ) 
Ijp  --  e/. change  current  density  for  cathodic  partial 

p 

electrode  reaction  for  reaction  95  (amp/cnr) 
[Uj  =  concentration  of  species  U  (mole/cm^) 


The  equilibrium  concentrations,  [M0*+]  and  [H+],  in 
equation  (99)  can  be  determined  from  the  solubility  product 
of  the  hydrated  oxide  or  from  the  free  energy  of  formation 
of  the  oxide  according  to  reaction  (95)  [ 1 8] . 

The  kinetic  parameters,  1^  and  (3,  appearing  in  equa¬ 
tions  (96)  and  (97)  have  been  evaluated  from  data  reported 
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fcr  titanium  by  H.  A.  Johansen  and  co-workers  [ 1 9] .  I1rt  has 
been  taken  as  1.0  x  10“  amps  per  square  centimeter  and  (3 

Q 

as  16  x  10  centimeters.  The  standard  Flade  potential,  6®  , 
has  been  taken  as  the  standard  equilibrium  potential  for  the 
formation  of  hydrated  titgnium  dioxide  by  the  reaction  of 
titanium  with  water.  The  value  reported  by  W.  M.  Latimer 
[13]  is  -37.1  (-0.95  volts  S.H.E.  at  25°C).  The  equilibrium 
concentrations  in  equation  (99)  have  been  evaluated  from  data 
reported  by  M.  Pourbaix  [18]  for  the  formation  of  the  hy¬ 
drated  form  of  titanium  dioxide  by  the  reaction  of  hydrogen 
and  titanyl  ions.  The  values  are: 

[Hj]  =  2[M0^+]  =  7.57  x  10-5  moles  per  cubic  centimeter. 


For  this 

■\++ 


A  detailed  study  of  the  rate  of  reaction  (95),  based  on 
the  equations  developed  in  Appendix  B,  indicates  that  the 
rate  is  negligible  with  respect  So  the  rate  of  reaction  (94) 
when  the  solution  pH  is  greater  than  minus  one. 
reason  and  also  to  standardize  the  influence  of  the  MO-1 

-l—L.  _Q 

ions,  the  concentration  of  MO  has  been  fixed  at  10  ?  moles 
per  cubic  centimeter  and  the  rate  of  reaction  (95)  has  been 
set  equal  to  zero  for  the  problems  in  this  section.  To  de¬ 
termine  ,  it  is  still  necessary  to  evaluate  the  constants 
and  I^Q.  There  are  no  data  in  the  literature  from  which 
to  evaluate  these  exchange  current  densities.  The  values 
10“^  and  10"4  amps  per  square  centimeter  have  been  selected 
as  being  reasonable. 

When  the  rate  of  reaction  (95)  is  negligible  with  re¬ 
spect  to  the  rate  of  reaction  (94),  the  rate  of  oxide  film 
growth  may  be  written  as 


(100) 
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where 

b  =  time  rate  of  change  of  the  film  thickness,?! 
(cm/sec) 

W  =  molecular  weight  of  the  oxide  (gram/mole) 

p  =  mass  density  of  the  oxide  (gram/cm^) 

&  =  Faraday's  constant  (coul/equiv) 

The  parameters  W  and  p  have  been  evaluated  from  the  chemi¬ 
cal  formula  and  density  of  titanium  dioxide. 

The  transport  properties  summarised  below  are  based  on 
data  contained  in  References  [20],  [21 ],  and  [22].  None  of 
the  references  is  complete;  the  data  reported  are  from 
measurements  made  in  binary  electrolytes.  The  values  below 
are  based  on  the  reported  values  at  infinite  dilution.  At 
concentrations  of  0.1  mole  per  liter,  ionic  mobilities  are 
typically  about  twenty  percent  less  than  at  infinite  dilu¬ 
tion. 


Dg-*-  =  9.00  x  IQ"5  (cm2/sec) 
DC1-  =  2.00  x  10“5  (cm2/sec) 
DNa+  =  x  10“5  (cm2/ see) 


Reference  Problem 


kt  any  point  in  time,  the  electrolyte  composition  in  a 
crack-like  region  is  determined  by  three  factors:  the  crack 
size  and  shape,  the  boundary  and  initial  conditions,  and  the 
transport  mode.  Each  of  these  is  examined  below  using  re¬ 
sults  obtained  numerically.  To  facilitate  comparison  of  the 
results,  we  begin  with  a  problem  that  is  relatively  simple, 
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but  illustrates  the  general  characteristics  of  transport  in 
crack-like  regions.  It  is  a  problem  of  transport  by  simple 
diffusion  under  constant  flux  boundary  conditions  in  a  hy¬ 
perbolic  crack. 

In  this  problem,  one  ion  is  considered;  the  effects  of 
migration  and  convection  are  ignored.  The  ion  is  the  hydro¬ 
gen  ion  and,  from  the  values  listed  above,  its  diffusion  co- 

-5 

efficient  is  9  x  10  centimeters  squared  per  second.  The 
governing  eifferential  equation  is  the  diffusion  equation, 
recovered  from  equation  (73)  by  setting  the  potential,  ve¬ 
locity,  and  generation  rate  equal  to  zero. 

It  is  assumed  that  the  electrolyte  is  initially  homo¬ 
geneous  and  the  composition  of  the  bulk  electrolyte  does  not 
change  with  time.  The  hydrogen  ion  concentration  at  time 

_7 

zero  and  at  points  along  the  crack  mouth  is  10  moles  per 
liter  (pH  =  ! j .  Along  the  crack  walls,  the  normal  flux  is 

_q 

constant  and  equal  to  3.25  x  10  moles  per  square  centi¬ 
meter.  This  value  is  determined  from  the  rate  of  reaction 
(94),  the  oxide-film  formation  reaction,  U3ing  equation  (96) 
with  (jv  and  equal  to  zero,  a  film  thickness  of  fifty  ang- 

*7 

stroma,  and  a  hydrogen  ion  concentration  of  10  '  moles  per 
liter.  The  normal  flux  is  equivalent  to  a  current  density 
of  314  microamps  per  square  centimeter. 

The  crack  is  the  hyperbolic  crack  shown  in  Figure  6. 

It  is  0.20  centimeters-long  and  0.0?  centimeters-wide. 

The  results  are  summarized  in  Figures  4  and  5.  In  Fig¬ 
ure  4,  the  hydrogen  ion  concentration  at  the  crack  tip  is 
shown  as  a  function  of  time.  The  concentration  is  expressed 
as  a  pH,  the  negative  common  logarithm  of  the  hydrogen  ion 
concentration  in  moles  per  liter.  The  pH  falls  quite  rapid¬ 
ly.  From  an  initial  value  of  seven,  it  is  lees  than  three 
within  one  second.  In  ten  seconds,  it  has  a  value  less  than 
one  pH  unit  larger  than  its  final  or  steady  state  value. 

The  steady  state  is  reached  in  approximately  700  seconds,  or 
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about  twelve  minutes  and  the  crack-tip  pH  is  very  close  to 
one. 

The  change  in  concentration  with  time  is  accompanied  by 
a  significant  change  in  the  shape  of  the  concentration  pro¬ 
file.  This  is  illustrated  in  Figure  5  where  the  concentra¬ 
tion  at  points  on  the  crack  centerline  is  shown  at  three  dif¬ 
ferent  times.  In  order  to  display  all  three  curves  on  a  sin¬ 
gle  graph,  the  data  are  presented  in  terms  of  the  dimension- 
less  parameter  C ,  defined  as 

I  Cam  -  C?  \ 

C;  =  ioTo!  -  cl)  Ooi) 

where 

CjOO  =  concentration  evaluated  at  the  distance  x 
from  the  crack  tip 

=  concentration  evaluated  at  the  crack  tip 

C?  =  concentration  at  the  mouth  of  the  crack  (equal 
to  the  concentration  in  the  bulk  electrolyte) 

The  lower  curve  is  a  typical  short-time  concentration 
profile.  The  concentration  gradient  is  large  in  the  region 
near  the  crack  tip,  approaches  zero  in  the  central  region, 
and  increases  again  near  the  mouth  of  the  crack.  The  cen¬ 
tral  curve,  corresponding  to  a  later  point  in  time,  behaves 
similarly,  but  the  change  in  the  concentration  gradient  is 
not  as  extreme  as  in  the  lower  curve.  The  upper  curve  cor¬ 
responds  to  the  steady  state.  The  gradient  is  small  at  the 
crack  tip  and  increases  continuously  as  the  distance  from 
the  crack  tip  increases. 

The  shape  of  the  concentration  profile  is  strongly  in¬ 
fluenced  by  +he  crack-like  geometry.  At  short  times  the  in¬ 
fluence  of  the  ends  is  small  except  in  the  regions  near  the 


tip  and  the  mouth  of  the  crack.  The  flux  in  the  central 
region  is  primarily  in  the  transverse  direction.  The  small 
gradient  in  the  central  region  of  the  lower  curve  in  Figure 
5  is  almost  entirely  due  to  an  increase  in  the  crack  width 
with  increasing  distance  from  the  crack  tip,  rather  than 
significant  transport  in  the  longitudinal  direction.  As 
the  steady  state  is  approached,  the  flux  in  the  longitudinal 
direction  becomes  much  greater  than  that  in  the  transverse 
direction  except  in  the  region  near  the  crack  tip.  There, 
both  components  remain  cf  the  same  order  of  magnitude. 


Crack  Shape 


The  effects  of  changes  in  th_  irack  shape  are  shown  in 
the  next  series  of  problem;:.  Except  for  the  crack  shapes, 
the  problems  are  identical  to  the  reference  problem.  One 
species,  the  hydrogen  ion,  is  considered.  Transport  is  by 
simple  diffusion.  The  don  flux  normal  to  the  electrode¬ 
electrolyte  interface  is  constant  and  equivalent  to  a  cur¬ 
rent  density  of  ?14  microamps  per  square  centimeter.  The 
bulk-solution  pH  is  seven. 

The  crack  profiles  are  shown  in  Figure  6.  They  are 
rectangular,  hyper oolic,  and  pie-shaped.  All  are  nominally 
0.2  centimeters  long  by  0.02  centimeters  wide.  The  pie¬ 
shaped  crack  is  one  percent  shorter  than  the  others  because 
of  the  coordinate  transformation  used.  The  area  of  metal- 
electrolyte  interface  is  virtually  the  same  for  all  of  the 
cracks;  the  maximum  difference  between  any  two  is  less  than 
five  percent.  The  crack  volumes  vary  by  a  factor  of  two. 

The  influence  of  shape  on  the  behavior  with  time  is  il¬ 
lustrated  in  Figure  7.  The  crack-tip  pH  is  shown  as  a  func¬ 
tion  of  time  for  each  of  the  three  crack  shapes.  For  the 
pie-shaped  crack  the  approach  to  steady  state  is  initi^  ly 
somewhat  faster,  for  the  rectangular  crack  the  approach  is 
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initially  slower,  than  for  the  hyperbolic  crack.  However, 
the  times  to  reach  the  steady  state  are  virtually  the  same, 
about  700  seconds  or  twelve  minutes.  The  steady-state 
crack-tip  pH  is  approximately  one  for  all  three  shapes.  The 
differences  are  more  apparent  when  the  values  are  expressed 
in  terms  of  concentration  rather  than  pH.  For  the  pie-shaped 
crack,  the  steady-state  crack-tip  concentration  is  0.137 
moles  per  liter.  For  the  hyperbolic  crack,  the  value  is 
0.099  moles  per  liter.  For  the  rectangular  crack  it  is 
C.,079  moles  per  liter. 

The  crack  shape  influences  the  form  of  the  concentra¬ 
tion  profile  as  shown  in  Figure  8.  The  steady-state  concen¬ 
tration  along  the  crack  centerline  is  plotted  as  a  function 
of  distance  from  the  crack  tip  for  each  of  the  three  crack 
shapes.  The  ordinate  in  Figure  8  is  C  as  defined  by  equa¬ 
tion  (101).  Thus,  only  differences  in  the  shapes  of  the 
concentration  profiles  are  shown  in  the  figure. 

For  the  pie-shaped  crack,  the  concentration  profile  is 

very  nearly  linear,  with  a  slight  deviation  from  linearity 

* 

in  the  region  close  to  the  crack  tip.  The  concentration 
profile  is  parabolic  for  the  rectangular  crack.  The  profile 
for  the  hyperbolic  crack  lies  between  the  other  two. 

The  influence  of  time  on  the  shape  of  the  concentration 
profile  is  liscussed  above  for  the  case  of  the  hyperbolic 
crack.  The  behavior  is  summarized  in  Figure  5  where  the  con¬ 
centration  on  the  crack  centerline  is  shown  at  three  differ¬ 
ent  times.  Similar  curves  are  presented  for  the  pie-shaped 
ci ...  ;k  in  Figure  9  and  for  the  rectangular  crack  in  Figure  10. 


* 

The  deviation  from  linearity  near  the  crack  tip  is  caused 
by  'chopping  off'  the  tip  of  the  wedge  in  making  the  coor¬ 
dinate  transformation.  It  may  be  shewn  that  tne  profile  is 
linear  over  the  entire  crack  length  only  when  the  cr^ck  tip 
coincides  with  che  p^int  of  intersection  of  the^exfcensions 
of  the  two  straight  sides,  (see  pages  95  end  96) 
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The  timed  are  the  same  as  those  in  Figure  5.  As  in  Figure  5> 
the  data  are  presented  in  terms  of  the  dimensionless  param¬ 
eter  C  in  order  to  display  the  curves  in  a  single  graph. 

Comparison  of  Figures  5»  9,  and  10  reveals  that  the  be¬ 
havior  with  time  is  similar  for  all  three  crack  shapes.  At 
short  times  there  is  negligible  transport  in  the  longitudinal 
direction  over  most  of  the  crack  length.  The  concentration 
gradient  is  large  in  the  region  near  the  crack  tip,  decreases 
to  a  small  value  in  the  central  region,  and  increases  some¬ 
what  near  the  mouth  of  the  crack.  At  intermediate  times  the 
behavior  is  similar,  but  the  change  in  the  gradient  over  the 
length  of  the  crack  is  less  extreme  than  at  the  shorter 
times.  At  large  times  the  concentration  profiles  approach 
steady  state  curves  characteristic  of  the  particular  crack 
shape . 

It  is  stated  above  that  the  concentration  gradient  in 
the  central  region  of  the  crack  at  short  times  is  primarily 
due  to  an  increase  in  the  crack  width  with  increasing  dis¬ 
tance  from  the  crack  tip.  This  is  apparent  from  the  lower 
curves  (elapsed  time  equal  to  0.1 64  seconds)  in  Figures  9 
and  10.  In  the  central  region  of  the  bottom  curve  in  Fig¬ 
ure  9,  the  concentration  decreases  linearly  with  increasing 
distance  from  the  crack  tip.  For  this  crack,  the  width  also 
increases  linearly.  In  the  central  region  of  the  bottom 
curve  in  Figure  10,  on  the  other  hand,  the  concentration  re¬ 
mains  constant.  Since  this  crack  is  rectangular,  its  width 
is  also  constant. 

Two-Dimensional  Effects 

The  problems  discussed  in  this  section  have  been  solved 
using  a  two-dimensional  formulation  of  the  transport  equa¬ 
tions.  For  problems  involving  transport  in  crack-like  re¬ 
gions,  ho r ever,  one  might  expect  transport  in  the  transverse 
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direction  to  play  a  relatively  minor  role.  This  expectation 
is  based  on  simple  dimensional  considerations.  For  a  more 
quantitative  view,  data  from  the  previous  set  of  problems 
may  be  used.  The  rectangular  crack  is  particularly  conve¬ 
nient  for  this  purpose  because  of  its  simple  geometry. 

The  transverse  concentration  gradient  varies  from  zero 
along  the  crack  centerline  to  a  finite  value  at  the  crack 
wall  (the  metal-electrolyte  phase  boundary).  At  the  wall, 
the  transverse  gradient  is  proportional  to  the  normal  flux 
and,  because  of  the  boundary  conditions,  remains  constant 
with  time.  The  longitudinal  gradient  varies  with  both  time 
and  position.  The  ratio  of  the  two  gradients  provides  a 
simple  and  direct  measure  of  their  relative  magnitudes. 

In  Figure  1 1  the  gradient  ratio  is  shown  as  a  function 
of  time  for  three  points  adjacent  to  the  crack  wall.  This 
ratio  is  defined  as  the  magnitude  of  the  longitudinal  con¬ 
centration  gradient  divided  by  the  transverse  concentration 
gradient.  The  points  are  located  at  distances  from  the 
crack  tip  equal  to  one-eighth,  one-fourth,  and  one-half 
times  the  overall  crack  length. 

At  short  times,  transport  in  the  transverse  direction 
is  much  greater  than  in  the  longitudinal  direction,  and  the 
gradient  ratios  are  small  at  all  three  points.  As  time  in¬ 
creases,  transport  in  the  longitudinal  direction  becomes  in¬ 
creasingly  important.  At  the  end  of  one  minute,  the  gradient 
ratios  are  all  of  the  order  of  one,  and  the  gradients  in  the 
longitudinal  and  transverse  directions  are  approximately 
equal.  As  time  increases  further,  the  ratio  at  each  point 
approaches  a  steady  state  value  which  is  a  function  of  dis¬ 
tance  from  the  crack  tip.  Because  the  crack  width  and  nor¬ 
mal  flux  are  both  constants,  the  steady  state  value  at  each 
point  is  proportional  to  the  total  wall  area  lying  between 
that  point  and  a  point  on  the  crack  centerline  at  the  crack 
tip.  The  steady  state  values  at  one-eighth,  one-fourth,  and 
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one-half  are,  therefore,  3*5>  6.0,  and  11.0,  respectively. 
Only  for  points  close  to  the  crack  tip  are  the  transverse 
and  longitudinal  concentration  gradients  of  the  same  order 
of  magnitude. 

Two-dimensional  transport  effects  may  also  be  viewed 
from  another  perspective.  Since  there  is  a  concentration 
gradient  in  the  transverse  direction,  the  concentration 
varies  from  point  to  point  on  a  section  perpendicular  to  the 
crack  axis.  As  a  measure  of  this  variation,  consider  the 

A 

parameter  C(Xo)  defined  as 


coo 


cjXo)  -  cc(yj 
CM  ~  c° 


(102) 


where 

CJx.)  =  concentration  at  a  point  adjacent  to  the  crack 
wall  and  located  at  a  distance  x0  from  the 
crack  tip 

Cc()0  =  concentration  at  a  point  on  the  crack  center- 
line  and  located  at  a  distance  X0  from  the 
crack  tip 

C°  =  concentration  in  the  bulk  electrolyte 

A 

When  C(aq)  is  one,  the  concentration  at  the  given  transverse 

A 

section  is  highly  non-uniform.  When  CU„)  is  zero,  the  con¬ 
centration  is  uniform  across  the  section.  Because  of  the 
assumed  symmetry  about  the  crack  centerline,  the  numerator 
in  equation  (102)  is  equal  to  the  difference  between  the  max' 
imum  and  minimum  concentrations  on  the  given  transverse  sec¬ 
tion. 


'The  functions  for  three  transverse  sections  are 

shown  as  functions  of  time  in  Figure  12.  The  sections  are 
located  at  distances  of  one-eighth,  one-fourth,  and  one-half 
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times  the  overall  crack  length.  The  curves  for  the  sec¬ 
tions  at  one-eighth  and  one-fourth  are  indistinguishable 
and  are  shown  as  a  single  curve.  The  curve  for  one-half 
lies  sl-ghtly  above  the  other  cvr-re. 

Initially,  the  concentration  changes  more  rapidly  at 
points  along  the  wall  than  at  points  away  from  it.  There 
is  a  substantial  variation  in  the  concentration  across  the 

A 

width  of  the  crack  and  the  function  C(X0)  is  close  to  one. 

As  time  increases,  the  concentration  becomes  more  uniform 
and  the  variation  across  the  width  of  the  crack  decreases. 
After  one  minute,  the  variation  is  about  one  percent.  As 
steady  state  is  approached  the  variation  becomes  quite 
email.  It  is  less  than  one  percent  at  all  three  sections. 

These  results  indicate  that  two-dimensional  transport 
effects  are  significant  only  at  short  times  or  at  points 
near  +he  ends  of  a  crack.  The  results  for  hyperbolic  and 
pie-shaped  cracks  are  similar.  In  the  case  of  the  gradient 
ratio,  however,  they  do  not  admit  as  simple  a  geometric  in¬ 
terpretation  as  is  the  case  for  the  rectangular  geometry. 

Crack  Length 

The  effect  of  crack  length  is  illustrated  by  the  next 
series  of  problems.  The  cracks  have  hyperbolic  profiles  and 
are  0.02  centimeters  wide.  The  influence  of  crack  width  is 
considered  elsewhere. 

Consideration  is  limited  -o  hydrogen  ion  transport  by 
simple  diffusion  with  a  constant  normal  flux  at  the  crack 
wall.  The  flux  is  determined  from  the  rate  of  the  film  for¬ 
mation  reaction  (eqn.  94)  with  the  potential  difference  at 

the  metal-electrolyte  interface  equal  to  zero  and  the  hydro- 

-7 

gen  ion  concentration  equal  to  10  moles  per  liter 
(pH  =  7.0).  The  oxide  film  thickness  is  fixed  at  fifty  ang¬ 
stroms.  This  corresponds  to  a  normal  current  density  of  314 
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microamps  per  square  centimeter.  The  hydrogen  ion  cor  co  - 

_7 

tration  in  the  bulk  electrolyte  is  10  moles  per  liter. 

At  this  point,  it  is  convenient  to  introduce  the  notion 
of  a  crack  aspect  ratio.  The  aspect  ratio  will  be  defined  as 
the  ratio  of  crack  length  to  crack  width  where  the  crack 
width  is  measured  at  the  mouth  of  the  crack.  In  this  prob¬ 
lem  set,  the  cracks  have  aspect  ratios  ranging  from  1.25  to 
40.  The  profiles  of  three  are  shown  in  Figure  13.  In  the 
upper  portion  of  the  figure,  th*  cracks  are  drawn  to  the 
same  scale  in  the  length  and  width  directions.  In  the  lower 
portion,  they  are  drawn  with  the  scale  in  the  width  direc¬ 
tion  ten  times  the  scale  in  the  length  direction.  The  cracks 
have  aspect  ratios  of  1.25,  10,  and  40.  The  crack  with  the 
aspect  ratio  of  ten  is  the  one  used  as  a  reference  problem. 
The  profiles  range  from  what  should  probably  be  called  a 
'dimple'  to  a  truly  crack-like  geometry. 

The  results  are  summarized  in  Figures  14,  15,  and  16. 

In  Figure  14,  the  crack-tip  pH  is  plotted  as  a  function  of 
time  for  each  of  the  geometries  shown  in  Figure  15.  Both  the 
time  to  reach  steady  state  and  the  steady  state  concentration 
increase  with  increasing  crack  length.  Even  for  the  shortest 
crack  (the  'dimple'),  the  steady  state  pH  is  less  than  three. 

In  Figure  15,  the  steady- state  crack -tip  pH  is  shown  as 
a  function  of  crack  length.  Note  that  the  scale  on  the  ab¬ 
scissa  is  logarithmic  rather  than  linear.  The  curve  exhibits 
a  slight  negative  curvature  for  small  values  of  crack  length 
but  approaches  linearity  as  the  crack  length  increases.  The 
curvature  at  small  values  of  the  crack  length  is  a  result 
both  of  the  function  plotted  and  the  method  used  to  define 
the  crack  domain.  As  the  crack  length  approaches  zero,  the 
pH  at  the  crack  tip  assymptotically  approaches  the  bulk  solu¬ 
tion  pH.  In  this  case  that  value  is  seven,  The  method  used 
to  define  the  crack  domain  and  if.s  effect  on  these  results  is 
discussed  below. 
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For  values  of  the  crack  aspect  ratio  much  greater  than 
one,  the  curve  becomes  linear.  In  this  range  a  doubling  of 
the  crack  length  reduces  the  crack-tip  pH  by  0.6.  In  tae 
next  section,  it  is  shown  that  the  slope  of  this  curve  is 
independent  of  the  normal  flux. 

The  steady- state  concentration  profiles  for  the  cracks 
in  Figure  13  are  shown  in  Figure  16.  The  data  have  been 
normalized  by  plotting  on  the  ordinate  the  function  0 ,  de¬ 
fined  by  equation  (101).  A  curve  has  been  drawn  through  the 
points  for  the  0.200- centimeter-long  crack.  No  attempt  has 
been  made  to  draw  curves  for  the  other  two  geometries,  since 
this  would  obscure  the  data  themselves.  It  may  be  seen  that 
all  the  data  very  nearly  fall  on  a  single  curve. 

The  data  for  the  0.025-centimeter-long  crack  fall 
slightly  below  those  for  the  other  two  geometries.  Like  the 
curvature  exhibited  at  small  crack  lengths  in  Figure  15, 
this  is  due  to  the  method  used  to  define  the  crack  domain. 
When  transforming  a  crack  onto  a  rectangle  in  the  computer 
program,  the  crack  is  bounded  on  the  solution  side,  not  by  a 
straight  line,  but  by  an  ellipse.  As  a  result,  the  crack 
length  measured  along  the  crack  centerline  from  the  tip  to 
the  bounding  ellipse  is  greater  than  the  length  measured  to 
the  metal  surface.  For  the  longer  cracks  bhe  difference  be¬ 
tween  the  two  is  negligible;  but  for  the  0.025 -centimeter- 
long  craci  the  difference  is  about  ten  percent. 

Boundary  Conditions 

In  the  problems  considered  thus  far,  the  flux  normal  to 
the  crack  wall  (metal-electrolyte  phase  boundary)  has  been 
treated  as  constant.  In  the  next  series  of  problems,  the 
effects  of  other  boundary  conditions  are  examined.  The 
cracks  are  C.200-oenticjetera  long  by  0.0?c-centimeters  wide 
and  hyperbolic  in  profile.  Consideration  is  limited  to 
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hydrogen  ion  transport  by  simple  diffusion. 

The  concentration  dependence  of  the  normal  flux  is  ex¬ 
amined  first.  The  flux  is  determined  from  the  rate  of  the 
film  formation  reaction  (equ.  94)  with  the  potential  differ¬ 
ence  across  the  metal-electrolyte  interface  equal  to  zero, 
and  the  oxide  film  thickness  equal  to  fifty  angstroms.  The 
current  density  is  calculated  from  equation  (96)  using  the 
hydrogen  ion  concentration  at  points  adjacent  to  the  crack 
wall.  Thus,  the  normal  flux  varies  both  with  time  and  posi¬ 
tion  along  the  crack  wall.  The  hydrogen  ion  concentration 

„7 

in  the  bulk  electrolyte  is  10  moles  per  liter. 

The  results  are  summarized  in  Figures  17,  18,  and  19. 

In  Figure  17,  the  crack-tip  pH  is  shown  as  a  function  of 
time.  The  results  of  the  reference  problem  are  included  for 
comparison.  The  curve  for  this  problem  lies  everywhere  a- 
beve  the  one  for  the  reference  problem.  As  these  results 
ehow,  including  the  concentration  dependence  of  the  normal 
flux  in  the  boundary  conditions  reduces  the  flux  and  in¬ 
creases  the  3teady  state  pH.  In  particular,  when  the  pH  at 
the  metal-electrolyte  interface  is  2.5,  the  normal  flux  is 
about  one  and  one-half  orders  of  magnitude  less  than  when 
the  pH  is  seven.  The  magnitude  of  this  effect  may  vary  from 
one  reaction  to  another,  but  the  direction  is  always  the 
same.  An  increase  in  the  concentration  of  a  product  species 
always  decreases  the  reaction  rate. 

The  variation  of  the  steady  state  concentration  along 
the  crack  centerline  is  shown  in  Figure  18.  The  results  ob¬ 
tained  for  the  reference  problem  are  also  shown  for  compari¬ 
son.  The  data  have  been  normalized  by  plotting  as  the  ordi¬ 
nate  the  function  C  defined  by  equation  (101).  The  concen¬ 
tration  profile  for  this  problem  is  very  nearly  the  same  as 
that  for  the  reference  problem.  The  former,  however,  lies 
everywhere  above  the  latter. 

Because  the  flux  normal  to  the  crack  wall  decreases  as 
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the  hydrogen  ion  concentration  increases,  it  is  greater  at 
the  mouth  of  the  crack  than  at  the  crack  tip.  This  is  shown 
in  Figure  19  where  the  normal  flux  is  plotted  as  a  function 
of  distance  from  the  crack  tip  at  steady  state.  For  conven¬ 
ience,  the  flux  is  expressed  as  a  current  density.  The  flux 
at  the  mouth  of  the  crack  is  more  than  eight  times  greater 
than  at  the  crack  tip. 

When  the  longitudinal  concentration  profile  is  compared 
with  one  obtained  for  a  constant  value  of  the  normal  flux, 
as  in  Figure  18,  it  may  be  seen  that  the  effect  of  the  vari¬ 
ation  in  the  normal  flux  is  to  reduce  the  magnitude  of  the 
longitudinal  concentration  gradient  near  the  crack  tip,  and 
to  increase  it  near  the  mouth  of  the  crack.  The  concentra¬ 
tion  profile  for  product  species  obtained  using  charge- 
transfer  boundary  conditions  is  always  found  to  lie  above 
that  obtained  using  constant  flux  boundary  conditions. 

The  same  boundary  conditions  at  the  crack  wall  have 
been  used  to  investigate  the  effect  of  changes  in  the  con¬ 
centration  of  the  bulk  electrolyte.  The  concentration  de¬ 
pendence  of  the  normal  flux  is  considered  but  the  oxide-filin 
thickness  is  fixed  at  fifty  angstroms.  The  potential  dif¬ 
ference  at  the  interface  is  set  equal  to  zero. 

The  results  are  summarized  in  Figure  20.  Here  the  pH 
at  the  crack  tip  is  shown  as  a  function  of  the  pH  in  the 
bulk  electrolyte.  There  are  two  distinct  regions  in  the 
figure.  In  the  first,  the  pH  in  the  bulk  electrolyte  is 
small  and  the  crack-tip  pH  is  approximately  equal  to  the  pH 
in  the  bulk  electrolyte.  In  this  region,  the  concentration 
difference  between  the  tip  and  the  mouth  of  the  crack  is 
much  less  than  the  concentration  at  the  mouth  of  the  crack. 
The  variation  in  the  pH  along  the  length  of  the  crack  is 
negligible  and  the  normal  flux  is  essentially  constant  every¬ 
where  along  the  crack  wall. 

In  the  second  region,  where  the  pH  in  the  bulk  electro- 
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lyte  is  large,  '.he  pH  at  the  crack  tip  ia  virtually  indepen¬ 
dent  of  the  pH  in  the  bulk  electrolyte.  In  this  region,  the 
difference  between  the  concentrations  at  the  tip  and  the 
mouth  of  the  crack  is  much  greater  than  the  concentration  at 
the  mouth  of  the  crack.  Thus,  the  hydrogen  ion  concentration 
in  the  bulk  electrolyte  may  be  treated  as  being  zero,  and  the 
variation  in  the  normal  flux  along  the  crack  wall  is  virtually 
the  same  as  that  shown  in  Figure  19. 

In  the  final  problem  of  this  set,  the  effects  of  both 
film  growth  and  concentration  are  considered.  The  normal 
flux  is  determined  from  the  rate  of  the  oxide-film  formation 
reaction  (eqn.  94)  with  the  potential  difference  across  the 
metal-electrolyte  interface  equal  to  zero.  The  hydrogen  ion 
concentration  in  the  rate  equation  (eqn.  96)  is  taken  as  the 
concentration  in  the  electrolyte  at  points  adjacent  to  the 
crack  wall.  The  initial  film  thickness  is  taken  as  fifty 
angstroms  at  every  point  on  the  boundary.  The  film  thick- 
nese  at  each  boundary  point  is  then  obtained  by  numerical 
integration  of  equation  (100). 

The  results  are  summarized  in  Figures  21  and  22.  In 
the  upper  portion  of  Figure  21 ,  the  crack-tip  pH  ia  shown 
as  a  function  of  time.  For  comparison,  a  curve  showing  the 
response  in  the  absence  of  film  growth  (corresponding  to  the 
upper  curve  in  Figure  17)  ie  also  shown.  In  the  lower  por¬ 
tion  of  the  figure,  the  film  thickness  at  the  crack  tip  is 
shown  as  a  function  of  time. 

At  very  short  times,  there  is  no  sensible  change  in  the 
oxide-film  thickness  and  the  pH  versus  time  curve  coincides 
with  that  obtained  using  a  constant  film  thickness.  With 
increasing  time,  the  thickness  of  the  oxide  film  increases, 
decreasing  the  hydrogen  ion  flux  into  the  crack.  The  con¬ 
centration  at  the  crack  tip  then  becomes  less  (the  pH  great¬ 
er)  in  tue  comparable  problen  with  no  film  growth. 

\».tinued  film  growth  results  in  the  most  striking  feature 
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of  the  pH- versus- time  curve.  At  about  one  hundred  seconds 
the  curve  exhibits  a  minimum  after  which  the  pH  rises  con¬ 
tinuously. 

This  calculation  was  terminated  after  an  elapsed  time 
of  about  three  thousand  seconds  or  fifty  minutes  (calculated 
time,  not  computer  running  time)  so  no  detailed  information 
was  obtained  on  the  behavior  at  large  times.  It  should  be 
observed,  however,  that  the  concentration  in  the  crack  oan- 
not  attain  a  steady  state  value  unless  the  rate  of  oxide 
film  growth  beccmes  zero.  Thus,  the  pH  will  continue  to 
rise  toward  the  bulk  electrolyte  value  (in  this  case  seven) 
unless  some  other  reaction  intervenes  or  the  potential  dif¬ 
ference  across  the  metal-electrolyte  interface  is  reduced  to 
its  equilibrium  value. 

For  the  first  problem  in  this  set,  the  oxide-film 
thickness  was  fixed  at  fifty  angstroms.  Under  this  con¬ 
straint,  the  normal  flux  was  found  to  exhibit  considerable 
variation  along  the  length  of  the  crack.  It  is  shown  in 
Figure  19,  for  example,  that  the  normal  flux  near  the  mouth 
of  the  crack  is  about  eight  times  greater  than  at  the  crack 
tip.  When  the  influence  of  oxide  film  growth  is  considered, 
however,  it  is  found  that  the  no.  mal  flux  at  the  crack  wall 
is  very  nearly  constant.  In  Figure  22,  the  normal  flux, 
expressed  as  a  current  density,  is  shown  as  a  function  of 
distance  from  the  crack  tip.  The  variation  in  the  film 
thickness  along  the  length  of  the  crack  is  also  shown.  The 
difference  in  the  normal  flux  between  the  tip  and  the  mouth 
of  the  crack  is  less  than  ten  percent.  The  greater  thick¬ 
ness  of  the  film  near  the  mouth  of  the  crack  tends  to  com¬ 
pensate  for  the  lower  ion  concentration  in  this  region. 


Convective  Transport 


Thus  far,  our  attention  has  been  restricted  to  problems 
involving  transport  by  simple  diffusion.  The  effects  of 
convection  and  migration  have  been  ignored.  We  now  wish  to 
consider  these  other  modes  of  transport.  The  effects  of 
convective  transport  are  considered  in  the  next  series  of 
problems.  Transport  by  migration  is  considered  in  the  fol¬ 
lowing  problem  set. 

When  a  structural  member  containing  a  crack  is  loaded, 
the  crack  opens  and  the  volume  of  the  crack  increases.  The 
crack  volume  likewise  increases  as  a  crack  grows  or  length¬ 
ens.  In  either  case,  some  of  the  electrolyte  is  drawn  into 
the  crack  from  the  bulk  solution  in  the  same  way  liquid  is 
drawn  into  a  suction  pump.  Our  interest  here  is  in  the 
effect  of  such  flow  on  the  composition  of  the  electrolyte 
in  the  crack. 

To  be  more  explicit,  consider  the  pie-shaped  crack 
shown  at  the  top  of  Figure  23-  If  the  crack  deforms  in  such 
a  way  as  to  remain  pie-shaped,  the  straight  sides  must  remain 
straight.  In  this  case,  it  is  convenient  to  identify  two 
modes  of  crack  deformation.  For  the  mode  identified  as 
’crack  growth'  in  Figure  23,  the  sides  of  the  crack  are 
displaced  vertically  with  no  change  in  tbe  crack  angle.  The 
ratio  of  the  vertical  displacement  of  a  point  on  the  crack 
wall  to  the  horizontal  displacement  of  the  crack  tip  i3  a 
constant  for  all  points  on  the  crack  wall.  For  the  mode 
identified  as  ’crack  opening'  in  Figure  23,  the  crack  angle 
increases  with  no  change  in  the  displacement  of  the  crack 
tip.  In  this  case,  the  displacement  of  a  point  on  the  crack 
vail  is  proportional  to  its  distance  from  the  crack  tip. 

Any  crack  deformation,  satisfying  the  requirement  that  the 
pie-shaped  crack  remain  pie-shaped,  can  be  expressed  as  a 
linear  combination  oi  these  two  deformation  modes. 


A  closed-form  solution  for  the  solvent  velocity  in  a 
pie-shaped  crack  has  been  obtained  using  the  Stokes-flow 
approximation  to  the  Navier-Stokes  equations.  When  the 
crack  angle,  2<x  ,  is  small,  the  solution  may  be  expressed 
approximately  as 


Vk 


V* 


where 


-i[x(a)  +  2d][i  “  (Si)  ] 


(103-a) 


(103-b) 


Uk  =  velocity  component  in  the  radial  (longitudinal) 
direction  (cm/sec) 

14  =  velocity  component  in  the  tangential  (transverse) 
direction  (cm/sec) 

a  =  crack  half  angle  (radian) 

0.  =  crack  length  (cm) 

a  =  time  rate  of  change  of  crack  half  angle 
(radian/sec) 

Q  =  time  rate  of  change  of  crack  length  (cm/sec) 

/l,©  =  coordinate  of  a  point  (cm,  radian) 


The  stream  lines  corresponding  to  the  velocities  in 
equations  (103)  are  shown  in  the  central  portion  of  Figure 
23  for  the  two  modes  of  crack  deformation.  In  the  lower 
portion  of  the  figure,  the  magnitudes  of  the  centerline 
velocities  are  shown  as  functions  of  the  distance  from  the 
crack  tip. 

For  both  deformation  modes,  the  solvent  flow  is  toward 
the  crack  tip.  At  a  given  distance  from  the  crack  tip,  the 
velocity  in  the  radial  direction  is  a  maximum  on  the  crack 
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centerline.  For  the  'crack  opening'  mode,  the  magnitude  of 
the  solvent  velocity  increases  linearly  as  the  distance  from 
the  crack  tip  increases.  For  the  'crack  growth'  mode,  the 
solvent  velocity  is  independent  of  the  distance  from  the 
crack  tip. 

The  results  obtained  for  the  pie-shaped  crack  provide 
a  good  approximation  to  solvent  motion  in  a  crack-like  re¬ 
gion.  For  the  numerical  work,  however,  a  different  geometry 
has  been  used.  This  geometry  is  somewhat  less  restrictive 
and  of  more  direct  physical  significance.  The  configuration 
adopted  is  shown  in  Figure  24.  It  consists  of  a  double 
edge-notched  tensile  specimen  2.40  centimeters  wide.  The 
side  notches  are  0.20-centimeter-long  by  0.02-centimeter¬ 
wide  hyperbolic  notches.  The  deformation  of  the  crack  walls 
is  determined  from  A.  A.  Griffith's  [23]  analysis  for  a 
specimen  containing  deep  hyperbolic  notches.  The  rate  of 
crack  deformation  is  found  to  be  directly  proportional  to 
the  specimen  loading  rate. 

In  this  series  of  problems  consideration  is  limited  to 
transport  of  the  hydrogen  ion  and  the  normal  flux  is  con¬ 
stant  at  every  point  along  the  crack  wall.  The  flux  is  de¬ 
termined  from  the  rate  of  the  film  formation  reaction 
(eqn.  94)  with  the  potential  difference  at  the  metal-elec¬ 
trolyte  interface  equal  to  zero  and  the  hydrogen  ion  con- 

_7 

centration  equal  to  10  moles  per  liter.  The  oxide-film 
thickness  is  fixed  at  fifty  angstroms.  Under  these  condi¬ 
tions  the  current  density  normal  to  the  crack  wall  is  314 
microamps  per  square  centimeter.  The  hydrogen  ion  concen- 

_7 

tration  in  the  bulk  electrolyte  is  taken  as  10  moles  per 
liter. 

The  solvent  velocity  is  calculated  using  a  prescribed 
specimen  loading  rate  and  the  initial  specimen  geometry. 

The  crack  dimensions  are  assumed  to  remain  constant  with 
time.  In  reality,  while  the  length  of  the  crack  is  not 
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affected  by  the  loading,  the  crack  width  increases  with 
time.  The  effect  of  changes  in  the  crack  width,  ignored 
in  this  set  of  problems,  is  examined  below. 

In  Figure  25,  the  magnitude  of  the  solvent  velocity 
along  the  crack  centerline  is  shown  as  a  function  of  dis¬ 
tance  from  the  crack  tip.  The  nominal  net-section  strain 
rate  is  10”^  centimeters  per  centimeter  per  second.  The 
velocity  is  for  all  intents  and  purposes  a  linear  function 
of  distance  from  the  crack  tip.  The  velocity  is  zero  at 
the  crack  tip  and  a  maximum  at  the  mouth  of  the  crack.  At 
the  mouth  of  the  crack,  the  centerline  velocity  is  3*34  x 
lO”5  centimeters  per  second.  Since  the  rate  of  deformation 
of  each  point  along  the  crack  wall  is  a  linear  function  of 
the  nominal  strain  rate,  the  velocity  at  any  other  strain 
rate  may  be  determined  by  simple  proportioning  from  this 
curve. 

With  the  nominal  net-section  strain  rate  equal  to  10”^ 
centimeters  per  centimeter  per  second,  the  difference  be¬ 
tween  the  results  obtained  with  convective  transport  and 
those  obtained  in  the  reference  problem  is  barely  percepti¬ 
ble.  With  the  strain  rate  increased  to  10~4  centimeters  per 
centimeter  per  second,  a  difference  is  readily  apparent. 

The  results  obtained  at  the  latter  strain  rate  are  summarized 
in  Figures  26  and  27. 

In  Figure  26,  the  crack-tip  pH  is  shown  as  a  function 
of  time.  The  solid  curve  corresponds  to  the  results  for 
the  convective-diffusion  problem  at  a  nominal  strain  rate 
of  10~4  centimeters  per  centimeter  per  second  (maximum 
solvent  velocity  of  3.34  x  1C centimeters  per  second). 

A  broken  curve,  illustrating  the  behavior  observed  for 
the  reference  problem,  is  included  for  comparison.  The 
results  for  the  two  problems  are  virtually  identical  for 
times  less  than  seventy  seconds,  about  one-tenth  of  the 
time  required  to  reach  steady  state.  For  times  greater  than 
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seventy  seconds,  the  crack-tip  pH  for  the  convective  trans¬ 
port  problem  is  somewhat  greater  than  that  for  the  reference 
problem.  At  steady  state,  the  ion  concentration  at  the  crack 
tip  for  the  convective  transport  problem  is  about  twenty-five 
percent  less  than  that  for  the  reference  problem. 

The  steady  state  concentration  along  the  crack  center- 
line  is  shown  as  a  function  of  distance  from  the  crack  tip 
in  Figure  27.  The  results  of  the  convective  transport  prob¬ 
lem  and  the  reference  problem  are  both  shown  in  the  figure. 
The  data  for  the  reference  problem  have  been  normalised  ir; 
the  same  way  as  in  previous  figures.  The  function  plotted 
as  the  ordinate  is  C  defined  by  equation  (101).  The  data 
for  the  convective  transport  problem  have  been  normalized  in 
the  same  way  except  that  the  value  C(0)  in  the  denominator 
in  equation  (101)  has  been  taken  as  the  value  obtained  in 
the  reference  problem.  In  other  words,  in  normalizing  the 
data  both  sets  have  been  divided  by  the  same  number. 

The  effect  of  the  solvent  motion  is  to  reduce  the  mag¬ 
nitude  of  the  concentration  gradient  and  thereby  to  reduce 
the  concentration  at  every  point  within  the  crack.  The  in¬ 
fluence  on  the  concentration  gradient  is  the  greatest  at  the 
mouth  of  the  crack  where  the  solvent  velocity  is  a  maximum. 
The  influence  decreases  as  the  distance  from  the  crack  tip 
decreases  and  the  solvent  velocity  decreases. 

Crack  Width 


As  mentioned  above,  the  effect;  cf  changes  in  the  crack 
geometry  with  time  was  ignored  in  generating  the  results 
summarized  in  Figures  26  and  27.  A  substantial  change  in 
the  crack  dimensions  may  occur,  however,  in  times  of  the 
order  of  the  time  required  to  reach  steady  state,  particu¬ 
larly  at  the  higher  strain  rates.  For  example,  when  the 
nominal  net-section  strain  rate  is  constant  and  equal  to 
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10-4  centimeters  per  centimeter  per  second  and  when  the 
specimen  configuration  is  that  shown  in  Figure  24,  the 
crack  width  increases  from  an  initial  value  of  0.02  centi¬ 
meters  to  0.26  centimeters  in  700  seconds.  Thus,  in  the 
time  required  to  reach  steady  state  the  crack  width  would 
have  increased  by  a  factor  of  thirteen. 

The  significance  of  such  a  change  may  be  judged  from 
the  results  of  a  set  of  problems  run  for  hyperbolic  cracks 
of  various  widths.  The  cracks  are  all  0.20  centimeters 
long.  Transport  is  by  simple  diffusion.  The  boundary  con¬ 
ditions  are  identical  to  those  used  in  the  reference  problem 
and  in  the  convective-diffusion  problem. 

The  results  are  summarized  in  Figure  28.  The  steady 
state  concentration  at  the  crack  tip  is  expressed  in  moles 
per  liter  and  plotted  as  a  function  of  crack  width.  For 
comparison  with  the  convective-diffusion  problem  discussed 
above,  the  steady  state  concentration  for  the  crack  0.20- 
centimeters  long  by  0.02-centimeters  wide  loaded  at  a  rate 
of  10"4  centimeters  per  centimeter  per  second  is  shown  in 
the  figure.  The  times  required  for  such  a  crack  to  attain 
the  same  width  as  the  widths  of  the  cracks  included  in  this 
problem  set  are  also  shown. 

The  steady  state  concentration  at  the  crack  tip  is  re¬ 
duced  by  a  factor  of  about  one-half  each  time  the  crack 
width  is  doubled.  By  comparison,  had  the  crack  width  in  the 
convective-diffusion  problem  been  allowed  tc  increase  rather 
than  being  fixe^  at  0.02  centimeters,  it  would  have  doubled 
in  about  sixty  seconds  or  less  than  one-tenth  of  the  time 
required  to  reach  steady  state. 

These  results  indicate  that  geometry  exerts  a  much 
stronger  influence  on  composition  than  does  solvent  motion. 
This  conclusion  is  strengthened  by  the  results  mentioned  a- 
bove  for  a  loading  rate  of  10“^  centimeters  per  centimeter 
per  second.  At  this  loading  rate,  the  effect  of  convection 
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was  barely  perceptible.  The  difference  between  the  concen¬ 
tration  calculated  for  this  problem  and  that  obtained  for 
the  reference  problem  was  everywhere  less  than  half  of  one 
percent.  Nevertheless,  for  times  of  the  order  of  the  time 
to  reach  steady  state,  the  width  of  this  crack  should  have 
more  than  doubled,  reducing  the  concentration  by  more  than 
one  half. 

Migration 

In  the  problems  discussed  thus  far,  transport  by  migra¬ 
tion  has  been  ignored.  The  results  are,  therefore,  strictly 
valid  only  when  the  electrostatic  potential  is  everywhere 
dimini shingly  small.  One  indication  of  the  range  of  appli¬ 
cability  of  these  results  is  provided  by  the  next  problem. 

A  broader  view  is  provided  by  the  analytic  results  presented 
in  the  next  section. 

The  problem  considered  here  is  one  of  electrochemical 
transport  in  a  rectangular  crack  0.20-centimeters  long  by 
0.02-centimeters  wide.  The  choice  of  the  rectangular  crack 
shape  is  based  on  the  results  presented  above  on  the  effect 
of  crack  shape  and  to  facilitate  a  later  comparison  with  an 
analytic  solution  to  a  similar  problem. 

The  electrolyte  consists  of  three  ions:  the  sodium 
ion  (Na+),  the  chloride  ion  (Cl“),  and  the  hydrogen  ion  (H+). 
The  bulk  electrolyte  is  a  0.6-mole  per  liter  sodium-chloride 
solution  with  a  pH  of  seven.  The  initial  composition  of  the 
electrolyte  is  the  same. 

At  the  crack  wall,  the  normal  fluxes  of  both  the  sodium 
ion  and  the  chloride  ion  are  zero.  The  hydrogen  ion  flux  is 
constant  and  has  the  same  value  as  in  the  reference  problem. 
That  is,  the  value  is  determined  from  the  rate  of  reaction 
(94)  with  the  potential  difference  across  the  metal-electro¬ 
lyte  interface  equal  to  zero,  the  pH  equal  to  seven,  and  the 
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oxide-film  thickness  equal  to  fifty  angstroms.  Therefore, 
the  normal  current  density  is  31 4  microamps  per  square  cen¬ 
timeter. 

The  behavior  with  time  is  summarized  in  Figure  29.  Two 
curves  are  shown.  The  first  shows  the  variation  in  the 
crack-tip  pH  with  time.  Perhaps  the  most  striking  aspect  of 
this  curve  is  its  similarity  to  the  curves  obtained  in  ear¬ 
lier  problems.  It  exhibits  no  new  or  unusual  features. 

The  second  curve  in  Figure  29  illustrates  the  behavior 
of  the  electrostatic  potential,  measured  at  the  crack  tip, 
with  time.  Since  the  electrostatic  potential  is  zero  at  the 
mouth  of  the  crack,  the  potential  at  the  crack  tip  is  numer¬ 
ically  equal  to  the  potential  drop  along  the  length  of  the 
crack.  Note  that  this  potential  drop  is  always  positive. 

That  is,  it  is  always  of  a  sense  to  carry  some  portion  of 
the  total  current. 

At  short  times,  the  gradient  of  the  hydrogen  ion  con¬ 
centration  is  small,  and  nearly  all  of  the  current  is  carried 
by  the  migration  of  ions  in  the  electrostatic  field.  The 
electrostatic  potential  remains  close  to  its  initial  value 
of  about  nine  millivolts.  As  the  hydrogen  ion  concentration 
reaches  a  value  of  about  one  percent  of  its  steady  state 
value,  the  potential  begins  to  decrease  rapidly.  As  the  hy¬ 
drogen  ion  concentration  increases  further,  an  increasing 
fraction  of  the  total  current  is  carried  by  hydrogen  ion  dif¬ 
fusion  and  the  electrostatic  potential  continues  to  drop.  At 
steady  state,  a  substantial  fraction  of  the  total  current  is 
carried  by  diffusion  and  the  electrostatic  potential  stabi¬ 
lizes  at  about  1.5  millivolts.  This  is  approximately  one- 
fifth  cf  the  initial  value. 

Two  unusual  features  of  the  potential  versus  time  curve 
should  be  mentioned.  At  short  times,  the  computed  value  of 
the  potential  rises  before  it  starts  to  fall  again,  intro¬ 
ducing  a  'bump'  in  the  curve.  This  'bump'  is  an  artifact 
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caused  by  truncation  error  in  calculating  the  boundary  con¬ 
ditions.  It  does  not  represent  the  true  behavior  of  the 
system,  so  the  solid  curve  is  not  drawn  through  this  data. 

At  steady  state,  the  calculated  values  are  erratic.  Varia¬ 
tions  in  the  computed  values  as  large  as  eight  percent  may 
be  observed.  This  behavior  also  derives  from  the  method 
used  to  construct  boundary  conditions.  As  indicated  in  Sec¬ 
tion  IV,  it  is  planned  to  modify  the  routines  used  in  the 
computation  of  the  boundary  conditions.  The  planned  modi¬ 
fication  should  eliminate  both  of  these  defects. 

The  influence  of  migration  on  the  behavior  of  the  hy¬ 
drogen  ion  may  be  judged  by  comparing  the  pH-versus-time 
curve  in  Figure  29  with  the  upper  curve  in  Figure  7.  That 
curve  was  obtained  using  the  same  crack  geometry  and  bound¬ 
ary  conditions  as  the  present  problem,  but  with  transport 
by  simple  diffusion  only.  The  two  curves  are  virtually 
identical.  There  is  a  difference  between  the  two,  but  the 
difference  is  extremely  small  when  the  concentrations  are 
expressed  in  terns  of  pH. 

The  effect  of  migration  is  shown  more  clearly  in  Figure 
30.  Here  the  steady-state  hydrogen  ion  concentration  along 
the  crack  centerline  is  shown  as  a  function  of  distance  from 
the  crack  tip.  In  addition  to  the  reeults  for  this  problem, 
the  results  obtained  fc  ■*  transport  by  simple  diffusion  are 
reproduced  from  Figure  10.  The  effect  of  migration  is  to 
reduce  the  hydrogen  ion  concentration  everywhere  within  the 
crack.  However,  the  reduction  is  quite  small.  In  the  pres¬ 
ent  case,  transport  by  both  migration  and  diffusion,  the 
concentration  Is  about  three  percent  less  than  in  the  case 
of  transport  by  diffusion  only. 

The  influence  of  the  potential  gradient  on  the  salt 
ions  iB  summarized  in  Figures  31  and  32.  In  Figure  31 ,  the 
changes  in  the  concentrations  of  all  three  ions  are  shown  as 
functions  of  time.  In  every  case,  the  quantity  plotted  on 


the  ordinate  is  the  difference  between  the  concentration  at 
the  crack  tip  and  the  initial  concentration  in  the  electro¬ 
lyte.  Alternately,  since  the  initial  concentration  of  each 
ion  is  the  same  as  its  concentration  in  the  bulk  electrolyte, 
the  ordinate  may  be  viewed  as  the  total  concentration  differ¬ 
ence  over  the  length  of  the  crack. 

The  curves  have  several  interesting  features.  Perhaps 
the  most  obvious  is  the  way  changes  in  the  ion  concentrations 
are  coupled  by  the  requirement  of  electroneutrality.  The  in¬ 
crease  in  the  hydrogen  ion  concentration  is  accompanied  by  an 
increase  in  the  concentration  of  the  negative  salt  ion  (Cl") 
and  a  decrease  in  the  concentration  of  the  positive  salt  ion 
(Na+).  In  addition,  the  hydrogen  ion  reaches  a  stable 
steady-state  value  in  less  time  than  does  the  positive  salt 
ion.  The  ions  do  not  respond  with  a  single  time  dependence. 
This  is,  perhaps,  not  surprising.  The  diffusion  i  efficients 
of  the  three  ions  are  unequal,  that  of  the  hydrogen  ion  being 
the  greatest  and  that  of  the  positive  salt  ion  the  least. 
However,  this  behavior  is  in  contrast  to  that  of  the  binary 
electrolyte  where  the  system  responds  with  a  single  time- 
dependence. 

The  difference  in  the  transport  properties  of  the  two 
positive  ions  is  also  reflected  in  the  behavior  of  the  nega¬ 
tive  salt  ion.  At  first,  the  concentration  of  the  negative 
ion  increases  with  time.  Then,  a  point  is  reached  after 
which  the  concentration  of  the  hydrogen  ion  remains  essen¬ 
tially  unchanged,  while  the  concentration  of  the  positive 
salt  ion  continues  to  decrease  under  the  influence  of  the 
electrostatic  potential.  The  concentration  of  the  negative 
ion  then  peases  through  a  maximum  and  subsequently  decreases. 
Finally,  the  concentration  of  the  negative  salt  ion  approaches 
a  steady  state  value  from  above. 

The  steady-state  concentration  profiles  are  chown  in 
Figure  32.  There  are  three  curves,  one  for  each  ion. 
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Plot  ted  on  the  ordinate  is  the  difference  between  the  con¬ 
centration  at  a  point  on  the  crack  centerline  and  the  con¬ 
centration  in  the  bulk  electrolyte.  The  abscissa  is  dis¬ 
tance  from  the  crack  tip.  The  concentration  profiles  are 
all  very  nearly  parabolic.  As  required  by  electroneutrality, 
the  sum  of  the  ordinates  is  everywhere  equal  to  zero. 

The  spatial  variation  of  the  electrostatic  potential  is 
illustrated  in  Figure  33*  The  initial  and  steady  state 
values  along  the  crack  centerline  are  shown  as  functions  of 
distance  from  the  crack  tip.  The  shapes  of  the  curves  are 
nearly  the  same.  The  initial  profile  is  parabolic  and  the 
final  profile  is  very  nearly  parabolic  as  well.  However, 
the  magnitude  of  the  potential  at  steady  state  is  much  less, 
five  times  smaller,  than  it  is  immediately  following  immer¬ 
sion  of  the  specimen. 

In  Figure  74.  the  current  density  measured  along  the 
crack  centerline  is  shown  as  a  function  of  distance  from  the 
crack  tip.  Since  the  normal  flux  at  the  crack  wall  is  con¬ 
stant,  this  curve  is  invariant  with  time.  The  curve  is 
linear  over  most  of  the  crack  length.  Near  the  mouth  of  the 
crack,  however,  there  is  a  moderate  deviation  from  linearity 
caused  by  two-dimensional  transport  effects. 

At  elapsed  times  which  are  small  with  respect  to  the 
time  required  to  reach  steady  state,  essentially  all  of  the 
current  is  carried  by  the  migration  of  ions  in  the  electro¬ 
static  field.  On  the  other  hand,  at  times  of  the  order  of 
the  time  to  reach  steady  state,  a  substantial  portion  of  the 
current  may  be  carried  by  ion  diffusion.  There  is  no  unique 
method  for  determining  the  fraction  of  the  total  current 
carried  by  diffusion  as  opposed  to  the  fraction  carried  by 
migration.  But,  we  can  obtain  a  measure  of  this  fraction  in 
a  way  that  is  reasonable  and  makes  some  physical  sense. 

One  way  to  estimate  the  fraction  of  the  total  current 
carried  by  diffusion  is  simply  to  partition  the  total  flux 


of  each  ion  into  a  diffusion  flux  and  a  migration  flux,  and 
consider  the  ratio  of  the  current  carried  by  the  diffusion 
fluxes  to  the  total  current.  Thus  we  might  define 


where 


(105-a) 

(105-b) 


The  partitioning  of  the  current  might  also  he  reasonably 
approached  from  another  point  of  view.  At  steady  state,  the 
total  flux  of  both  salt  ions  must  everywhere  be  equal  to 
zero.  This  follows  from  the  boundary  conditions  and  conser¬ 
vation  of  mass.  Thus,  the  current  must  all  be  carried  by 
hydrogen  ions,  and  we  might  consider  the  fraction  of  the 
total  current  carried  by  hydrogen  ion  diffusion,  namely 


006) 


where 

io  ~  -^h-LWCV  (107) 

These  two  measures  of  the  fraction  of  the  total  current 
carried  by  diffusion,  P^  and  P2,  are  shown  in  Figure  34. 

The  fractions  have  been  evaluated  from  the  steady  state  data 
shown  in  Figures  32  and  33,  and  are  plotted  as  functions  of 
distance  from  the  crack  tip.  The  value  of  P^  is  everywhere 
less  than  P~,  but  both  measures  indicate  that  a  greater 


-89- 


fraction  of  the  total  current  ie  carried  by  diffusion  as  the 
distance  from  the  crack  tip  increases.  Also,  both  measures 
indicate  that  a  greater  portion  of  the  total  current  is 
carried  by  diffusion  than  is  carried  by  migration. 


-giri-mi-i 
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VI.  DISCUSSION 

The  results  presented  in  the  last  section  covered  a 
broad  range  of  geometric  and  kinetic  variables.  The  effect 
of  each  was  examined  separately  and  in  some  detail.  In  this 
section  a  more  general  perspective  is  adopted.  First,  based 
on  the  numerical  results  of  the  last  section,  the  character¬ 
istics  of  transport  in  crack-like  regions  are  summarized. 
Consideration  is  then  centered  on  the  extent,  if  any,  to 
which  this  characterization  is  of  more  general  validity. 

Like  other  computer  programs,  the  program  developed  in 
this  study  is  capable  of  providing  detailed  answers  to  par¬ 
ticular  questions.  It  is  ill-suited  to  and  inefficiently 
employed  in  a  parametric  analysis  of  a  wide  range  of  a  large 
number  of  variables.  Therefore,  it  will  be  convenient  in 
this  discussion  to  draw  upon  the  one-dimensional  steady-state 
transport  analysis  developed  in  Appendix  C.  While  its  formu¬ 
lation  constrains  the  application  of  this  analysis,  it  has 
the  important  advantage  that  many  useful  results  may  be  ex¬ 
pressed  in  a  simple,  closed  form. 

From  the  results  of  the  previous  section,  the  following 
general  observations  can  be  made: 

1 )  Time  -  The  initial  rate  of  change  of  the  crack-tip  pH  is 

quite  large.  For  times  greater  than  or  equal  to  the 
* 

diffusion  time,  however,  the  time  rate  of  change  of 
the  pH  is  small  regardless  of  the  form  of  the  boundary 
conditions.  Referring  to  Figure  4,  for  example,  with 
a  constant  normal  flux  at  the  crack  wall  sufficient  to 
produce  a  steady-state  crack-tip  pH  equal  to  one,  the 


*  /  2  /  \ 

The  diffusion  time  is  defined  here  as  (a  /D)  where  a  is 

the  crack  length  and  D  is  the  diffusion  coefficient  of  the 

dissolved  species  of  interest. 


crack-tip  pH  is  less  than  five  in  one  millisecond  and 
less  than  two  in  ten  seconds.  These  times  should  he 
compared  to  the  diffusion  time  of  444  seconds  or  the 
time  to  reach  steady  state  of  700  seconds. 

Under  more  complex  boundary  conditions,  the  steady 
state  may  not  be  attained  in  times  of  the  order  of  the 
diffusion  time.  Shown  in  Figure  21,  for  example,  are 
results  obtained  using  charge-transfer  boundary  con¬ 
ditions  with  oxide  film  growth.  Even  here,  the  time- 
rate-of-change  of  the  pH  is  small  for  times  of  the 
order  of  the  diffusion  time  and  larger.  At  the  end  of 
one  hour,  for  example,  the  pH  differs  by  less  than  one 
unit  from  the  value  at  one  minute. 

While  this  observation  has  been  expressed  in  terms 
of  the  pH,  it  applies  equally  to  the  common  logarithm 
of  the  concentration  of  any  dissolved  species.  It  is 
asserted,  not  that  changes  in  concentration  with  time 
are  insignificant,  but  that  the  order  of  magnitude  of 
the  concentration  is  slowly  time-varying  over  most  of 
the  observable  time  domain. 

2)  Two-Dimensional  Transport  Effects  -  The  influence  of 
transport  in  the  transverse  direction  as  compared  to 
the  longitudinal  direction  may  be  judged  from  the  var¬ 
iation  in  the  electrolyte  composition  across  the  width 
of  the  crack.  If  a  variation  of  less  than  one  percent 
is  regarded  as  negligible,  the  effects  of  transport  in 
the  transverse  direction  are  negligible  at  times  great¬ 
er  than  the  diffusion  time  for  all  cracks  having  an 
aspect  ratio  greater  than  three.  At  large  times  and 
in  regions  near  the  ends  of  a  crack,  composition  var¬ 
iations  due  to  transport  in  the  transverse  direction 
may  be  observed,  but  these  are  small  as  compared  with 
the  variations  along  the  length  of  the  crack. 
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3)  Crack  Geometry,  Shape  -  The  shape  of  a  crack  has  a  modest 

but  definite  effect  on  the  behavior  of  the  electrolyte 
contained  within  it.  The  shape  affects  both  the  magni¬ 
tude  of  the  crack  tip  concentration  and  the  shape  of 
the  longitudinal  concentration  profile.  For  example, 
under  identical  constant  flux  boundary  conditions,  the 
steady-3tate  concentration  difference  between  the  tip 
and  the  mouth  of  a  pie- shaped  crack  is  nearly  twice  as 
great  as  for  a  rectangular  crack  having  the  same  length 
and  width.  For  the  same  boundary  conditions,  the 
steady-state  longitudinal  concentration  profile  is 
approximately  linear  for  the  pie-shaped  crack  and  para¬ 
bolic  for  the  rectangular  crack.  The  results  for  a  hy¬ 
perbolic  crack  lie  between  those  for  the  other  two 
crack  shapes. 

4)  Crack  Geometry,  Length  and  Width  -  The  geometric  param¬ 

eters  having  the  greatest  effect  on  electrolyte  compo¬ 
sition  are  the  length  and  the  width  of  the  crack.  The 
results  for  hyperbolic  cracks  with  constant  flux  bound¬ 
ary  conditions  may  be  cited  as  an  example.  To  a  good 
first  approximation,  the  steady-state  concentration 
difference  between  the  tip  and  the  mouth  of  the  crack 
is  found  to  be  proportional  to  the  square  of  the  crack 
length  and  the  reciprocal  of  the  crack  width  (see  eqn. 
Ill ). 

5)  Boundary  Specification  Along  the  Crack  Wall  -  The  form  of 

the  boundary  conditions  along  the  crack  wall  (metal- 
electrolyte  interface)  may  have  a  significant  influence 
on  the  concentration  of  an  ion  as  measured  by  its  aver¬ 
age  or  maximum  value.  It  has  little  effect,  however, 
on  the  shape  of  the  concentration  profile.  This  is 
illustrated,  for  example,  by  the  results  shown  in 
Figures  18  and  19. 
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6)  Composition  of  the  Bulk  Electrolyte  -  The  electrolyte 

composition  in  the  hulk  solution  may  bear  no  direct 
relation  to  that  within  a  crack-like  region.  This  is 
illustrated  by  the  results  shown  in  Figure  20. 

7)  Transport  by  Convection  -  The  solvent  motion  associated 

with  crack  opening  or  'yawning'  has  a  negligible  in¬ 
fluence  as  compared  to  the  effects  of  crack  geometry 
and  transport  by  simple  diffusion. 

8)  Transport  by  Migration  -  Transport  by  migration  in  an 

electrostatic  field  appears  to  have  little  influence 
on  the  concentration  cf  the  product  species.  The 
principal  effect  of  migration  is  to  alter  the  concen¬ 
trations  of  the  salt  ions  in  accordance  with  the  re¬ 
quirement  of  electroneutrality. 

The  first  two  observations  are  perhaps  not  surprising 
since  they  are  consistent  with  simple  dimensional  arguments. 
The  remainder  are  based  solely  on  numerical  results,  those 
presented  in  the  last  section  and  many  others  not  shown  ex¬ 
plicitly.  They  point  to  the  conclusion  that  the  factors 
which  exert  a  primary  influence  on  the  concentration  of  the 
ith  species  in  a  crack-like  region  are:  the  mean  value 

of  the  normal  flux  of  the  ith  species;  a,  the  crack  length; 
and  b,  the  crack  half-width.  Other  factors,  such  as  crack 
shape  or  transport  mode,  appear  to  exert  a  secondary  influ¬ 
ence. 

While  the  above  observations  represent  an  accurate  sum¬ 
mary  of  the  numerical  results,  the  data  upon  which  they  are 
based  is  limited.  The  last  observation  in  particular  is 
based  on  the  results  of  a  single  problem.  There  remains  the 
question,  therefore,  of  whether  these  observations  are  gen¬ 
erally  applicable  to  transport  in  crack-like  regions. 

To  answer  this  question  requires  further  information. 
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Rather  than  generate  additional  numerical  results,  it  is  more 
convenient  to  exploit  the  observation  that  composition  varia¬ 
tions  in  the  transverse  crack  direction  are  small  with  re¬ 
spect  to  variations  in  the  longitudinal  direction.  This 
suggests  that  a  crack-like  region  may  be  treated  with  rea¬ 
sonable  accuracy  as  a  one-dimensional  continuum.  The  reduc¬ 
tion  of  the  problem  to  dependence  on  one  spatial  variable 
results  in  a  simplified  set  of  governing  differential  equa¬ 
tions  as  shown  in  Appendix  C.  The  advantage  of  this  approach 
is  that  closed-form  solutions  may  be  obtained  for  a  number  of 
cases  important  to  our  discussion. 

Simple  Diffusion  -  Dominant  Effects 

Suppose  we  begin  by  looking  at  transport  by  diffusion 
only.  Consider  a  single  species,  S^,  in  a  rectangular  crack 
with  a  constant  flux,  J?,  normal  to  the  crack  wall.  From 
equations  (C-3a ,  -8  and  -9)  in  Appendix  C,  the  steady  state 
concentration  is 


fc-cn  = 


r- 1 , , 
r,kruof->rt; 


(a- a 


(108) 


where 

C(  =  concentration  of  ith  dissolved  species  (mole/cm^) 
C?  =  concentration  of  ith  species  in  the  bulk 

•z 

electrolyte  (mole/cnr) 

J ,n=  normal  flux  of  ith  species  at  the  crack  wall 

p 

(positive  inward)  (mole/cm  -sec) 

p 

D<  =  diffusion  coefficient  of  ith  species  (cm  /sec) 


a  =  crack  length  (cm) 


b  =  crack  half-width  (cm) 

X  =  distance  from  the  crack  tip  (cm);  [o  <=■  x  &  a] 


A  similar  result  may  he  obtained  for  transport  by  dif¬ 
fusion  within  a  pie-shaped  crack.  Prom  equations  (C-3b,  -8 
and  -9)  in  Appendix  C  with  a  constant  normal  flux,  J?,  the 
steady-state  concentration  of  the  ith  species  is 


(C|-c5)  - 


.ON  r 


(' 


b'2  (o.  -  x) 


+  eacvi. 


■Hd;] 


009) 


where 


€  =  a  number  between  zero  and  one,  typically  small 

X  =  distance  from  the  point  of  intersection  of  the 

extensions  of  the  straight  sides  of  the  crack  (cm) 

[ea  ^  x  a  a] 


It  should  be  noted  that  the  normal  flux,  J^,  appearing 
in  equations  (108)  and  (109)  is  taken  as  positive  for  a  mass 
flux  INTO  the  crack.  Also,  in  obtaining  equation  (109),  the 
crack  tip  is  located  a  distance  ea  from  the  point  of  inter¬ 
section  of  the  extensions  of  the  straight  sides  of  the  crack. 
Thus,  in  general,  the  width  of  the  crack  at  the  crack  tip  is 
non-zero  and  the  crack  profile  is  similar  to  that  of  the  pie¬ 
shaped  crack  shown  in  Figure  6. 

Fi*>m  equation  (108),  the  steady-state  concentration  pro¬ 
file  for  a  rectangular  crack  is  parabolic.  From  equation 
(109),  the  steady-state  concentration  profile  in  a  pie-shaped 
crack  is  very  nearly  linear.  The  deviation  from  linearity 
is  significant  only  for  values  of  x  close  to  ea.  when  the 
second  term  in  the  brackets  on  the  right-hand  side  of  equa¬ 
tion  (109)  is  no  longer  negligible.  When  e,  in  equation 
(109),  is  zero,  the  width  of  the  pie-shaped  crack  is  zero  at 
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the  crack  tip,  and  equation  (109)  reduces  to 

(Q-c  /)  =  (iVTt^(a-x) 


(109-a) 


The  steady  state  concentration  is  a  linear  function  of 
distance  from  the  crack  tip. 

The  steady-state  concentration  difference  between  the 
tip  and  the  mouth  of  the  crack  is  obtained  by  setting  x 
equal  to  zero  in  equation  (108)  or  (109-a).  For  a  rectan¬ 
gular  crack, 


For  a  pie-shaped  crack, 


(110-a) 


(110-b) 


When  the 
the  relations 
approximately 


aspect  ratio,  (a/b),  is  much  greater  than  one, 
expressed  by  equations  (110)  may  be  expressed 
by 


AG/ 


(111) 


where 

a  =  a  shape  factor,  and 

a  1 ,  for  a  rectangular  region 

cx  =  2,  for  a  pie-shaped  region 
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Equation  (111)  is  entirely  consistent  with  the  results 
obtained  numerically  for  transport  by  simple  diffusion  with 
constant  flux  boundary  conditions.  The  steady-state  concen¬ 
tration  difference  between  the  tip  and  the  mouth  of  a  crack 
is  directly  proportional  to  the  normal  flux  and  to  the  square 
of  the  crack  length.  It  is  Inversely  proportional  to  the 
diffusion  coefficient  and  the  width  at  the  mouth  of  the 
crack.  The  constant  ex,  viewed  as  a  shape  factor,  will  have 
a  value  between  one  and  two  for  all  crack  geometries  likely 
to  be  found  in  nature . 

Boundary  Conditions 

While  equation  (111)  adequately  summarizes  the  results 
obtained  for  simple  diffusion  with  constant  flux  boundary 
conditions,  it  does  not  explain  the  behavior  observed  with 
charge-transfer  boundary  conditions.  (The  important  feature 
of  such  a  boundary  specification  is  that  the  reaction  rate, 
or  reaction  current  density,  changes  with  changes  in  the 
concentrations  of  the  dissolved  species.  For  the  reaction 
used  in  the  numerical  calculations,  the  hydrogen  ion  was  the 
only  dissolved  species  affecting  the  reaction  rate  and  its 
normal  flux  was  directly  proportional  to  the  reaction  current 
density. )  Since  the  hydrogen  ion  concentration  was  a  func¬ 
tion  of  distance  from  the  crack  tip,  one  would  expect  its 
normal  flux  to  likewise  be  a  function  of  distance  and  this 
was  in  fact  observed.  With  the  normal  flux  varying  from 
point  to  point  along  the  crack  wall,  one  would  further  ex¬ 
pect  the  shape  of  the  concentration  profile  to  differ  from 
that  of  a  similar  curve  obtained  using  constant  flux  bound¬ 
ary  conditions.  This  was  also  observed,  but  the  difference 
was  small  (see  Figure  18). 

At  this  point  our  concern  is  not  with  the  validity  of 
the  numerical  results.  The  question  is,  do  those  results 
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reflect  the  particular  kinetic  formulation  used  in  that 
work,  or  does  the  form  of  the  boundary  conditions  generally 
have  a  small  effect  on  the  shape  of  the  concentration  pro¬ 
file?  Here,  as  above,  the  one-dimensional  transport  analy¬ 
sis  developed  in  Appendix  C  is  useful  in  framing  an  answer. 
Consider  a  single  dissolved  species,  contained  in  a 
rectangular  crack  and  assume  that  transport  is  by  simple 
diffusion.  Prom  equations  (C-8  and  -9)  in  Appendix  C,  the 
steady  state  concentration  is 

^  A 

(Ci-cfi)  =  4{j‘(a-x)  +  -b\[\ (112) 

4  u  g  X  Jo 

where 

J*  =  normal  flux  of  the  ith  species  at  the  crack 
tip  (mole/cm  -sec) 

=  normal  flux  of  the  ith  species  along  the  crack 

wall  at  a  distance,  x,  from  the  crack  tip 

t  2  » 

(mole/cm  -sec) 

=  dummy  variables 

Note  that  both  fluxes  are  positive  when  they  correspond  to 
flows  INTO  the  crack.  Also,  when  the  metal-electrolyte 
interface  is  everywhere  homogeneous,  is  equal  to  J^(0). 

The  normal  flux,  J*,  may  be  treated  two  different  ways. 
It  may  be  viewed  either  as  an  explicit  function  of  the  dis¬ 
tance  x  or  as  a  function  of  concentration  (an  implicit 
function  of  x).  Neither  approach  is  completely  satisfactory. 
When  the  flux  is  treated  as  a  function  of  concentration,  it 
is  highly  desirable  (although,  perhaps,  not  absolutely 
necessary)  that  the  functional  form  used  be  consistent  with 
some  charge-transfer  kinetic  relation.  Normally  then,  the 
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first  step  is  to  specify  a  charge-transfer  reaction  and 
kinetic  equation.  Once  this  is  done,  the  functional  rela¬ 
tion  between  normal  flux  and  concentration  may  be  determined 
and  the  result  used  with  equation  (112)  to  determine  the 
concentration.  Unfortunately,  it  is  usually  not  possible  to 
obtain  a  closed-form  solution  to  the  resulting  integral  equa¬ 
tion,  so  our  first  objective,  that  of  obtaining  a  simple  but 
completely  general  result,  cannot  be  achieved  in  this  manner. 
Nevertheless,  a  closed-form  solution  can  be  obtained  for  one 
important  case  and  this  is  to  be  examined  below.  First,  it 
will  be  useful  to  consider  the  other  alternative. 

When  the  flux  is  treated  as  a  function  of  distance, 
evaluation  of  the  integral  on  the  right-hand  side  of  equa¬ 
tion  (112)  is  straightforward  and  yields  an  explicit  ex¬ 
pression  for  the  concentration.  Such  a  formulation  has  the 
advantage  of  simplicity,  but  explains  neither  how  nor  why 
the  flux  actually  varies  under  a  given  set  of  boundary  con¬ 
ditions.  However,  it  does  provide  a  means  for  establishing 
bounds  cn  the  effect  of  differences  in  the  form  of  the  flux 
specification.  Specifically,  by  treat inw  the  normal  flux  as 
a  function  of  distance,  it  is  possible  to  establish  absolute 
upper  and  lower  bounds  on  the  shape  of  the  concentration  pro¬ 
file. 

To  establish  bounds  of  this  type  it  is  not  necessary  to 

specify  the  form  of  jV(x).  It  will  be  assumed,  however,  xhat 

w 

J_j  is  expressible  as  an  explicit  function  of  x  and  is  every- 

1  t 

where  of  the  same  sign  as  J^.  Following  notation  introduced 
above,  the  shape  of  the  concentration  profile  may  be  conven¬ 
iently  uaracterized  by  the  parameter  C  defined  by  equation 
(101).  Therefore,  replacing  C^(x)  in  equation  (101)  by  its 
value  in  equation  (112),  the  shape  of  the  concentration  pro¬ 
file  is  given  by 


r.\No  -  x)  +  I ( X ;  1 

'■  iv  *  i*  J 


(113) 
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where 


Itx) 


a  «, 

x  n 


(114-a) 


(114-b) 


Several  features  of  equation  (113)  should  be  emphasized. 
First,  it  may  be  Been  that  I(x)  is  a  monotonic  function  of 
x  and  is  bounded  by 


I(X) '  < 

i 


(114-c) 


Also,  since  jV  is  everywhere  of  the  same  sxgn  a3  jt',  l(x)  is 
monotonically  increasing  or  decreasing  according  as  is 
negative  or  positive.  From  this  it  follows  that  Q-  is  char¬ 
acteristically  positive  regardless  of  the  sign  of  the  flux, 
that  is,  without  regard  to  whether  the  normal  flux  is  into  or 
out  of  the  crack. 

The  bounds  we  seek  are  simply  the  upper  and  lower  limits 
of  the  right-hand  side  of  equation  (113).  For  example,  the 
lower  bound  is  approached  as  the  integral  I  becomes  small. 
That  is, 


Lc 


(1^5) 


The  lower  limit  is  approached,  therefore,  when  the  flux  at 
the  crack  tip  is  muon  greater  than  the  flux  at  any  point 
along  the  crack  wall.  This  behavior  may  be  observed  physi¬ 
cally  when  the  composition  of  the  metal  at  the  crack  tip  is 
significantly  different  from  the  composition  along  the  crack 
wall. 

*  t 

When  the  integral  I  becomes  large,  or  the  flux  be¬ 
comes  small,  a  different  behavior  is  observed.  The  limit  as 


-101- 


+ 

becomes  small  is 

L,  =  ^0(6)  -  (,„ 

The  desired  upper  bound  corresponds  to  the  upper  bound  of 
defined  by  equation  (11 6— a) .  Prom  the  inequality 
(114-c)  it  follows  that 


L-u  =  Max(L-i)  --  t 


(116-b) 


The  upper  limit  is  approached  when  the  flux  near  the  mouth 
of  the  crack  is  much  greater  than  the  flux  at  the  crack  tip 
or  other  points  along  the  crack  wall.  This  behavior  may  be 
observed  physically  when  the  concentration  in  the  crack  ap¬ 
proaches  its  equilibrium  value. 

Therefore,  when  the  normal  flux  is  everywhere  of  the 
same  sign,  the  shape  of  the  concentration  profile  lies  be¬ 
tween  the  limits 


rii-'AL  . 
LQ(o)  -  CJJ 


(117) 


In  addition,  the  graph  of  the  function  Ch(x)  is  a  continuous 
curve  and  has  a  slope  which  is  everywhere  of  the  same  sign. 

Similar  limits  may  be  obtained  for  cracks  having  other 
than  a  rectangular  shape.  Such  limits  provide  no  detailed 
view  of  the  conditions  existing  in  cracks,  but  they  serve 
to  emphasize  that  the  concentration  of  a  dissolved  species 
is  of  the  same  order  of  magnitude  nearly  everywhere  within 
a  crack.  In  addition,  these  limits  and  the  general  form  of 
charge-transfer  kinetic  equations  strongly  suggest  that  the 
normal  flux  is  likewise,  nearly  everywhere,  of  the  same 
order  of  magnitude. 

An  indication  of  how  the  concentration  actually  varies 
within  the  limits  described  by  equation  (117)  can  be  obtained 


from  the  solution  of  equation  (112)  using  a  simple  relation 
between  concentration  and  flux.  First,  however,  it  should 
be  noted  that,  for  any  charge-transfer  reaction,  the  func¬ 
tional  relation  between  normal  flux  and  concentration  is 
always  expressible  in  the  general  form 

Jin  =  c %  -  (3/CiP  (118) 

The  coefficients  c and  (3[  may  be  functions  of  the  potential 
difference  across  the  metal-electrolyte  interface  and  the 
concentrations  of  other  dissolved  species.  The  constant  p 
is  typically  a  positive  integer  but  can  be  any  positive  num¬ 
ber.  The  general  solution  of  equations  (112)  and  (118)  is 
beyond  the  scope  of  the  present  work.  Here  it  is  assumed 
that  of i  and  (3t-  are  functions  only  of  the  potential  differ¬ 
ence  across  the  metal-electrolyte  interface.  Further,  -in 
order  to  display  a  closed-form  solution,  the  constant  p  iB 
taken  as  one.  For  our  purposes  then,  equation  (118)  may  be 
written 

J?  =  Piic“-Ci)  (119) 

where 

cP  -  ai/fe 

Equation  (119)  represents  a  substantial  simplification 
of  equation  (118).  Still,  it  provides  an  adequate  descrip¬ 
tion  of  a  broad  class  of  charge-transfer  reactions.  In  this 
class,  for  example,  are  many  dissolution-precipitation  re¬ 
actions.  Before  continuing,  it  may  be  useful  to  examine  a 
specific  reaction  in  order  to  show  the  relationship  between 
equation  (119)  and  the  usual  kinetic  description  of  charge- 
transfer  reactions.  This  will  also  provide  a  basis  for 
evaluating  the  parameters  (3<  and  A  metal-ion  dissolu- 
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tion  reaction  serves  as  a  good  example.  Such  a  reaction  may¬ 
be  written 

M  M2*  +  Z+£T  (120) 

Equation  (120)  may  be  viewed  as  a  specialisation  of 
equation  (40).  For  this  reaction,  the  reduced  species  (M) 
is  insoluble  and  the  charge-transfer  valence  is  the  same  as 
the  charge  (z+)  on  the  oxidized  species  (Mz+,  the  metal  ion), 
letting  the  subscript  i  represent  the  metal  ion  and  follow¬ 
ing  the  notation  of  Section  II,  the  normal  flux  is 

{  6Xp  [d  Z^]  -  ,g  ex  p  (i  ■ -  «)  Z^]}  (121) 

where 

%  =  021-a) 

and 

Io  =  standard  exchange  current  density  (amp/cnr) 

ot  =  charge-transfer  coefficient 

Cf  =  standard  equilibrium  concentration  (mole/cm^) 

(J)m  =  dimensionless  potential  of  the  metal 

=  dimensionless  potential  of  the  electrolyte  at  the 
metal-electrolyte  interface 

6*  =  dimensionless  standard  electrode  potential 

8 

By  convention,  the  standard  conditions  are  taken  with 
equal  to  one  mole  per  liter. 

By  direct  comparison  of  equation  (119)  and  (121), 
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ia  i  r  -I 

or  = 


ft 


=  ( 


iLyc/V1-^ 

c? ) 


(122-a) 

(122-b) 


The  equilibrium  concentration,  C®^,  is  e  function  of 
the  dimensionless  overpctential.  It  increases  as  the  ap¬ 
plied  potential  increases.  The  slope,  (3<,  depends  on  both 
the  exchange  current  density  and  the  equilibrium  concentra¬ 
tion.  It  increases  as  the  exchange  current  density  in¬ 
creases  and  decreases  as  the  equilibrium  concentration  in¬ 
creases.  Prom  a  comparison  of  equations  (122-a)  and  (122-b), 
it  is  clear  that  the  slope,  (3/,  can  be  written  as  an  explicit 
function  of  the  dimensionless  overpotential.  In  many  re¬ 
spects  such  a  form  is  preferable.  However,  the  equation  is 
written  the  way  it  is  to  emphasize  the  relationship  between 
the  slope  and  the  equilibrium  concentration. 

When  in  equation  (119)  is  used  to  express  and 
in  equation  (112),  it  may  be  shown  that  the  concentration  is 


(C<  -  Q 


coch(  >  ,  ^ 


(123) 


where 

x  =  V^Ku- 

»  =  (b>'  JL* 

(v  a 


(124-a) 

(124-b) 

(124-c) 


A  simpler  form  may  be  used  when  the  crack  aspect  ratio 
is  large. 
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(C/  Ci)  (Lsi  CJl  1  co o h  v  a  )  ■* 


(125) 


For  reasons  of  simplicity,  equation  (125)  will  be  used 
through  the  remainder  of  this  discussion. 


To  compare  this  result  with  those  obtained  earlier,  it 
is  again  convenient  to  use  the  function  (*,:  defined  by  equa¬ 
tion  (101).  From  equations  (101)  and  (125), 


j~co?>h(X<i  -  c 

L  coohOo  —  i  J 


(126) 


The  behavior  described  by  equation  (126)  is  summarized 
in  Figure  35.  The  function  C£  is  shown  for  several  values 
of  the  parameter  X.  When  X  is  less  than  one,  the  deviation 
from  the  parabolic  profile  characteristic  of  constant  flux 
boundary  conditions  is  small.  An  increase  of  X  by  one-and- 
one-half  orders  of  magnitude  results  in  a  concentration 
profile  which  is  virtually  flat  over  ninety  percent  of  the 
crack  length.  Only  in  a  region  close  to  the  mouth  of  the 
crack  does  the  concentration  differ  significantly  from  the 
value  at  the  crack  tip. 

The  normal  flux  can  be  determined  from  equations  (119) 
and  (125).  After  introducing  J,  ,  a  non-dimensionalized 
flux,  for  later  notational  convenience,  the  result  is 


rvj 

J,  - 


folcp-q T) 


ooqh(  XVI 

COOX  •  >  V 


(127) 


The  behavior  of  the  normal  flux  is  summarized  in  Figure 
36.  The  dimensionless  flux,  *3^  ,  is  shown  as  a  function  of  £> 
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for  the  same  values  of  X  as  those  used  in  Figure  35.  When  X 
is  less  than  one-third,  the  normal  flux  is  essentially  con¬ 
stant  over  the  entire  crack  length.  Increasing  X  hy  two 
orders  of  magnitude  reduces  the  normal  flux  to  a  value  close 
to  zero  over  most  of  the  crack  length.  Only  in  the  region 
near  the  mouth  of  the  crack  is  the  flux  appreciable.  Com¬ 
parison  of  Figures  35  and  36  shows  that  the  change  in  the 
shape  of  the  concentration  profile  from  parabolic  to  flat  is 
associated  with  a  change  in  the  distribution  of  the  normal 
flux  along  the  length  of  the  crack. 

The  change  from  a  parabolic  to  a  flat  concentration 
profile  occurs  over  a  relatively  narrow  range  of  X  values. 
The  reason  for  the  change  is  shown  in  Figure  37.  The  lower 
curve  is  a  graph  of  the  concentration  ratio,  Q^,  as  a  func¬ 
tion  of  X  where  is  defined  as  the  ratio  of  the  concentra¬ 
tion  difference  between  the  tip  and  the  mouth  of  the  crack 
to  the  concentration  difference  at  equilibrium.  That  is, 


rco^WQ.'i  -  h 
L  cob'mX'l  J 


(128) 


For  small  values  of  X,  the  concentration  ratio,  Q^,  increases 
as  the  square  of  X .  As  X  approaches  one,  begins  to  dev¬ 

iate  from  the  curve  for  small  A  values.  For  values  of  X 
greater  than  10,  is  independent  of  X  and  equal  to  one. 
Therefore,  che  change  from  a  parabolic  to  a  flat  concentra¬ 
tion  profile  occurs  when  the  concentration  in  the  crack  ap¬ 
proaches  its  equilibrium  value. 

The  practical  significance  of  these  results  depends  on 
the  range  of  X  values  likely  to  be  encountered  physically. 

An  estimate  of  this  range  can  be  constructed  from  estimates 
of  the  parameters  in  equation  (124-a).  However,  a  more  con¬ 
venient  form  m ay  be  obtained  by  replacing  (3;,  in  equation 
(124-a)  by  its  value  in  equation  (122-b)  and  setting  the 
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charge-transfer  coefficient  equal  to  one-half.  That  is, 


X*  - 


(129) 


Values  of  the  exchange  current  density,  l“,  are  typi- 

_q  _•*  0 

cally  of  the  order  of  10  to  10  J  amps  per  square  centi- 

-S 

meter.  Ion  diffusion  coefficients  are  about  10  centi¬ 
meters  squared  per  second  at  room  temperature.  The  equi¬ 
librium  concentration  may  vary  over  a  wide  range  depending 

_Q 

on  the  value  of  the  overpotential,  £,  but  the  range,  10 

to  10"1  moles  per  cubic  centimeter,  appears  reasonable.  The 

crack  dimensions  may  also  vary  over  a  wide  range.  However, 

a  crack  length  of  one  centimeter  and  a  crack  aspect  ratio  of 

one  hundred  are  typical.  Faraday's  constant  is  about  10 

coulombs  per  equivalent.  Ion  valences  are  all  of  the  order 

of  magnitude  of  one.  By  convention,  the  standard  equilib- 

rium  concentration  is  10  y  moles  per  cubic  centimeter. 

Using  these  values  in  equation  (129),  the  range  of  X  values 

is  found  to  be  10”  ^  to  10+^. 

A  very  small  or  very  large  value  of  X  characterizes  a 

reaction  which  is  displsced  far  from  equilibrium.  Very  small 

values  correspond  to  large  positive  overpotentials  and  thus, 

to  large  positive  values  of  the  normal  flux,  J?.  Very  large 

values  correspond  to  large  negative  overpotentials  and  thus, 

to  large  negative  values  of  J1?.  For  reactions  close  to  equi- 

1  -2 

librium,  the  value  of  X  is  typically  in  the  range,  10  to 


The  limits,  X  equal  to  zero  and  X  equal  to  infinity, 
correspond  respectively  to  positive  and  negative  values  of 
the  normal  flux.  However,  it  should  be  cautioned  that  the 
sign  of  the  flux  does  not  necessarily  determine  the  shape  of 
the  concentration  profile.  The  sign  of  the  flux  is  deter¬ 
mined  by  the  sign  of  (C?^  -  C?) .  That  is,  the  flux  is 
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positive  or  negative  according  as  the  equilibrium  concentra¬ 
tion  is  greater  than  or  less  than  the  concentration  in  the 
bulk  electrolyte.  On  the  other  hand,  the  shape  of  the  con¬ 
centration  profile  depends  only  on  the  value  of  A.  That 
value  is  in  turn  a  function  of  C®^,  but  it  is  independent  of 
the  concentration  in  the  bulk  electrolyte. 

Taken  together,  the  above  results  indicate  that  treat¬ 
ing  the  normal  flux  as  a  constant  is  an  accurate  approxima¬ 
tion  for  a  wide  range  of  A  values.  More  precisely,  when  X 
is  less  than  or  equal  to  one,  the  error  introduced  by  using 
a  parabolic  approximation  to  Si  is  everywhere  less  than  one 
percent.  In  terms  of  the  kinetic  parameters  of  equation 

r«j 

(121),  the  error  in  Q  is  less  than  one  percent  when 

cxp[ll-ctes]  >  (z”t§fy  (150) 


Of  the  results  thus  far,  two  are  particularly  striking. 
The  first  is  the  bounds  given  by  equation  (117).  The  second 
is  the  close  conformity  of  the  actual  concentration  profile 
to  the  parabolic  approximation  over  a  wide  range  of  A  values. 
These  results  suggest  that  a  satisfactory  approximation  to 
the  concentration  in  a  crack-like  region  may  be  obtained  from 
an  equation  of  the  form  of  equation  (118)  or  (119),  and  the 
average  concentration  calculated  from  equation  (108).  For 

example,  assuming  that  the  crack  aspect  ratio  is  large  and 

n  ^ 
treating  a a  a  constant,  the  average  concentration,  C^, 

calculated  from  equation  (108)  is 


Ct  =  C°i  + 


-tD.  , 

Ji  a 

Sblh 


( 1 3 1 -a ) 


/  \ 

Replacing  in  equation  (119)  by  C^, 


(131-b) 


Equations  (121)  constitute  a  system  of  two  equations  in 

♦  « 

two  unknowns.  Letting  Jk  be  the  value  of  J£,  which  satisfies 
equations  (121),  the  solution  of  this  system  of  equations  is 


C? 

j; 


=  c?  +  >L(c?~cil 

1  '2  +  A* 


Oiler--  cn 

"t-  ^3  - 


(132-a) 


(132-b) 


The  procedure  is  essentially  the  same  when  the  flux  is 
given  by  an  equation  of  the  form  of  equation  (128).  How¬ 
ever,  the  system  of  equations  will  usually  be  larger  and  may 
be  non-linear. 

The  agreement  between  equation  (132-a)  and  the  average 
concentration  from  equation  (125)  is  excellent.  Equation 
(132-a)  is  assymptotically  correct  for  both  small  and  large 
A  values.  The  maximum  error  is  about  twelve  percent  at  a  N 
value  of  about  four. 

In  contrast,  equation  (132-b)  is  correct  only  for  small 
values  of  A.  Letting  be  the  average  value  of  the  normal 
flux  calculated  from  equations  (119)  and  (125),  the  ratio, 
is 


IV 

J*/ 


=  A'i  ♦ 


A 


(133) 


This  flux  ratio  is  shown  as  a  function  of  A  by  the  upper 

curve  in  Figure  37.  For  A  values  less  than  one  the  ratio 

is  essentially  one.  For  larger  A  values,  the  ratio  increases 

* 

with  increasing  A,  so  underestimates  the  true  average 
value  by  increasingly  greater  amounts.  This  occurs  because 
the  true  value  of  the  flux  always  remains  different,  from 
zero  in  a  small  region  near  the  mouth  of  the  crack.  As  a 
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result  t  the  true  value  of  the  average  flux  approaches  zero 

more  slowly  with  increasing  X  than  does  J., . 

*  1  _ 

Whether  the  discrepancy  between  Ji  and  is  significant 
is  a  moot  question.  Because  of  the  shape  of  the  concentra¬ 
tion  profile  for  large  values  of  A,  the  average  flux, 

provides  a  less  accurate  measure  of  the  reaction  rate  in  the 

*  * 

crack  interior  than  does  J^.  On  the  other  hand,  provides 
a  poor  estimate  of  the  total  mass  flux  entering  or  leaving 
through  the  mouth  of  the  crack. 

This  analysis  is  consistent  with  the  numerical  results 
and  places  them  in  a  broader  perspective.  It  shows  that  the 
concentration  profile  is  nearly  parabolic  over  a  wide  range 
of  conditions.  Significant  deviations  from  the  parabolic 
shape  do  occur  for  values  of  X  greater  than  one,  but  for 
values  greater  than  thirty,  the  concentration  profile  is 
virtually  flat .  Orly  for  A  in  the  range  between  one  and 
thirty  does  the  shape  of  the  concentration  profile  differ 
significantly  from  one  of  the  two  limits,  so  the  parameter 
X  may  be  used  as  a  quick  and  simple  indicator  of  the  shape 
of  the  concentration  profile.  Finally,  it  has  been  shown 
that  accurate  estimates  of  the  concentration  and  the  flux  in 
the  interior  of  a  crack-like  region  may  be  obtained  from  the 
assumption  that  the  normal  flux  is  uniform.  The  accuracies 
of  these  estimates  are  virtually  independent  of  the  act\ial 
shape  of  the  concentration  profile. 

The  numerical  results,  the  bounds  summarized  by  equa¬ 
tion  (H7),  and  the  one-dimensional  analysis  presented  above, 
form  a  coherent  and  fairly  complete  description  of  the  inter¬ 
action  between  normal  flux  and  concentration  in  a  crack-like 
region.  The  results  of  the  numerical  calculations  and  those 
of  the  one-dimens.ional  analysis  provide  epecific  examples  of 
how  the  normal  flux  actually  varies  from  point  to  point  along 
the  crack  wall  depending  on  the  reaction  kinetics,  the  ap¬ 
plied  potential,  and  the  composition  of  the  bulk  electrolyte. 


These  results  indicate  that  the  shape  of  the  concentration 
profile  is  very  nearly  parabolic  over  a  wide  range  of  con¬ 
ditions.  In  addition,  the  results  of  the  one-dimensional 
analysis  in  general,  and  equation  (127)  in  particular,  em¬ 
phasize  that  when  the  metal-electrolyte  interface  is  homo¬ 
geneous,  the  normal  flux  is  everywhere  of  the  same  sign. 
Finally,  the  results  most  important  to  an  understanding  of 
transport  in  crack-like  regions  are  the  bounds  summarized 
by  equation  (117).  These  bounds,  together  with  the  fact  that 
the  concentration  is  a  smooth,  continuous  function  of  dis¬ 
tance,  indicate  why  the  form  of  the  boundary  conditions  does 
not  have  a  strong  influence  cn  the  shape  of  the  concentration 
profile.  They  also  provide  a  theoretical  basis  for  the  use 
of  approximations  such  as  equations  (132)  to  estimate  the 
composition  within  a  crack-like  region. 


Convective  Transport 

The  numerical  work  indicated  that  transport  by  diffusion 
plays  a  dominant  role  in  determining  the  concentration  in  a 
crack-like  region.  Transport  by  convection  and  migration 
appeared  to  have  little  effect  on  the  results,  either  at 
short  times  or  long  times.  If  this  is  generally  uruo,  ths 
complexity  of  many  electrochemical  transport  problems  may  be 
substantially  reduced.  If  not,  it  may  still  be  possible  to 
generate  useful  results  using  a  diffusion  analysis.  In  this 
case,  however,  it  will  first  be  necessary  to  determine  the 
range  over  which  the  effects  of  migration  and  convection  are 
negligible  and  to  estimate  bounds  on  the  errors  introduce! 
when  they  are  not. 

The  results  of  the  last  section  indicate  that  convection 
plays  a  minor  role  in  the  case  of  a  hyperbolic  crack  sub¬ 
jected  to  the  opening  mode  of  crack  deformation.  That  con¬ 
vective  transport  is  almost  always  unimportant  may  be  shown 
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from  the  one-dimensional  transport  analysis  in  Appendix  C. 
Consider  first  the  opening  mode  of  crack  deformation  and, 
for  simplicity,  assume  that  the  crack  is  rectangular.  If 
the  crack  walls  move  outward  with  a  uriform  velocity,  It  , 
the  rectangular  region  remains  rectangular,  and  the  average 
longitudinal  velocity,  V ,  is 

where 

b  =  b0  +  vit 

and 

b0  =  value  of  the  crack  width,  b,  at  time  zero  (cm) 
t  =  time  vsec) 

£,  =  dimensionless  distance  from  the  crack  tip  (x/a) 

Note  that  the  average  longitudinal  velocity,  Z/t  is  a  linear 
function  of  distance  from  the  crack  tip.  The  same  result  is 
obtained  for  the  opening  mode  of  deformation  regardless  of 
the  crack  shape. 
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Prom  equations  (C-2  and  -4)  in  Appendix  C,  the  differ¬ 
ential  equation  describint  -nvect.-ive  transport  in  a  rectan¬ 
gular  crack  is 
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where 


r  =  D<£ 

z/a/Dt 


(137-a) 

(137-b) 
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As  the  coefficient  tf*  approaches  zero,  the  convective 
term  approaches  zero,  and  equation  (136)  approaches  the 
differential  equation  governing  transport  by  simple  diffu¬ 
sion.  When  V*  is  greater  than  zero,  omission  of  the  con¬ 
vective  term  introduces  an  error.  For  example,  when  If*  is 
one-tentl:  „  the  error  introduced  by  neglect jng  the  convective 
term  is  about  two  percent;  when  If*  is  one,  the  error  is 
about  fifteen  percent.  By  comparison,  the  maximum  value  of 
V*  in  the  convective-diffusion  problems  of  the  last  section 
was  about  seven.  For  that  value  of  V  ,  the  deviation  of  the 
solution  from  the  solution  for  simple  diffusion  was  about 
twenty-five  percent. 

At  short  times,  the  concentration  at  the  crack  tip  is 
influenced  only  by  conditions  existing  in  a  region  close  to 
the  tip.  But,  by  equation  (134),  the  velocity  there  is  zero. 
Therefore,  equation  (136)  reduces  to  the  equation  describing 
transport  by  simple  diffusion  and,  at  short  times,  the  sol¬ 
vent  motion  has  no  effect  on  the  concentration  at  the  crack 
tip.  This  was  observed  in  the  results  presented  in  the  last 
section. 

At  large  times,  motion  of  the  solvent  affects  the  con¬ 
centration  everywhere  within  the  crack.  But,  from  the  def- 
initi  '  of  V*  and  equations  (154)  and  (135),  it  may  be  shown 
that 


!?/♦!<  P  '  Al 
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(138) 


In  particular,  for  times  of  the  order  of  the  time  to  reach 
the  steady  state, 
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and 


\l/*\  4  <  1  (139-b)  . 

Thus,  for  times  of  the  order  of  the  time  to  reach  the  steady- 
state,  the  error  introduced  by  neglecting  the  convsctive 
term  is  less  than  twenty  percent.  The  error  may  be  greater 
than  this  at  shorter  times,  but  it  is  certainly  less  at 
larger  times. 

For  technical  alloys,  the  crack  deformation  rates  which 
can  be  sustained  for  extended  periods  of  time  are  limited  by 
metal  fracture.  Thus,  for  problems  of  practical  interest, 
the  crack  deformation  rates  will  be  small  and  the  errors 
will  be  less  than  that  corresponding  to  the  absolute  upper 
bound  in  equation  (139-b).  In  this  connection,  it  should  be 
emphasized  that  a  V*  value  of  seven  was  obtained  in  the  nu¬ 
merical  calculations  only  by  ignoring  changes  with  time  in 
the  crack  dimensions. 

For  the  mode  of  crack  deformation  we  have  called  the 
crack  growth  mode,  the  longitudinal  solvent  velocity  is  con¬ 
stant  along  the  length  of  the  crack  and  the  velocity  profile 
is  approximately  parabolic  at  every  cross-section.  As  a  re¬ 
sult,  the  velocity  relative  to  the  crack  tip  is  zero,  and, 
according  to  the  one-dimensional  transport  analysis,  the 
governing  differential  equation  reduces  to  that  for  trans¬ 
port  by  simple  diffusion.  Thus,  when  the  boundary  condi¬ 
tions  are  expressed  as  functions  of  distance  from  the  crack 
tip  (a  moving  coordinate),  the  analysis  of  problems  involving 
the  growth  mode  of  crack  deformation  is  identical  to  that  for 
transport  by  simple  diffusion  in  a  stationary  crack. 

The  above  results  show  that  the  convective  term  may  be 
omitted  in  almost  all  problems  involving  transport  in  a 
crack-like  region.  For  the  opening  mode  of  crack  deformation, 


the  error  introduced  by  neglecting  the  convective  term  does 
not  exceed  twenty  percent  for  times  greater  than  (a  /ll). 

For  the  growth  mode  of  crack  deformation,  no  error  is  intro¬ 
duced  by  neglecting  convective  transport  relative  to  the 
crack  tip,  since  this  contribution  is  in  fact  zero. 

Migration 

The  results  of  only  one  electrochemical  transport  prob¬ 
lem  were  presented  in  the  last  section.  Those  results  were 
consistent  with  the  view  that  diffusion  is  the  dominant  mode 
of  ion  transport.  However,  none  of  the  data  presented  indi¬ 
cates  the  range  over  which  this  observation  may  be  valid. 

It  is  not  clear,  for  example,  that  the  diffusion  equation 
even  constitutes  a  useful  approximation  outside  the  range  of 
very  small  current  densities. 

To  evaluate  the  importance  of  transport  by  migration 
and  to  understand  the  role  of  the  electrostatic  potential  in 
electrochemical  transport  requires  a  more  detailed  analysis. 
Once  again,  the  one-dimensional  transport  analysis  provides 
a  suitable  framework.  It  has  been  possible  to  obtain  a  gen¬ 
eral  solution  to  a  class  of  problems  which  includes  the 
problem  solved  numerically  as  a  special  case.  A  detailed 
discussion  of  this  solution  is  contained  in  Appendix  D. 
Because  the  solution  is  expressible  in  closed-form  and  is 
not  too  complicated,  it  is  ideally  suited  to  our  present 
purposes. 

It  is  assumed  that  the  crack  is  rectangular  and  the 
system  is  at  steady  state.  Transport  is  by  diffusion  and 
migration.  The  electrolyte  contains  three  dissolved  species 
'  with  the  charge  numbers 
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Species  one  and  two  will  be  called  the  salt  ions  and 
species  three  will  be  called  the  product  ion.  By  analogy 
with  the  problem  solved  numerically,  the  species  one,  two, 
and  three  may  be  associated  with  the  ions  Na+,  01“,  and  H+ 
respectively.  The  subscript  notation  is  adopted  here  both 
for  convenience  and  to  emphasize  that  the  results  are  in  no 
way  restricted  to  a  particular  electrolyte. 

The  boundary  conditions  are  similar  to  those  used  in 
the  problem  solved  numerically.  No  restriction  is  imposed 
on  the  normal  flux  of  species  three.  The  normal  flux  of 
species  one  and  two  is  taken  as  zero  at  every  point  along 
the  metal-electrolyte  interface.  In  the  notation  introduced 
above , 

Ji  =  J?  =  0 ,  i  -  1,2  (HI ) 

From  equations  (D-6,  -10,  -11,  and  -12)  in  Appendix  D, 
the  steady  state  solution  is 


Ci  = 

C$/UM 

(142) 

c*  = 

(c;  +  cro) 

(143) 

(c;  -!■  ccuoa  -  Ci/u(x) 

(144) 

4)  - 

|ri[uw] 

(145) 
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The  correctness  of  this  solution  may  be  verified,  by  direct 
substitution  of  equations  (142  -  147)  into  the  flux  equation 
and  the  electroneutrality  equation  (eqns.  C-4  and  -5  in 
Appendix  C).  Note  that  the  dimensionless  potential,  has 
been  set  equal  to  zero  at  the  mouth  of  the  crack  (x  =  a) . 

This  has  been  done  for  convenience  arid  consistency  with  the 
numerical  work,  since  the  absolute  value  of  the  potential  has 
no  effect  on  the  results. 

Several  features  of  the  solution  deserve  special  empha¬ 
sis.  The  concentration  of  species  one,  the  positive  salt 
ion,  is  inversely  proportional  to  the  function  u(x) ;  the 
concentration  of  species  two,  the  negative  salt  ion,  is  di¬ 
rectly  proportional  to  u(x).  This  suggests  that  the  solu¬ 
tion  may  be  unbounded  as  u(x)  approaches  either  zero  or  in¬ 
finity.  It  may  also  be  seen  that  the  diffusion  coefficients 
of  the  salt  ions  do  net  appear  in  the  equations.  In  fact, 
the  steady  state  solution  never  includes  the  diffusion  co¬ 
efficient  of  a  species  whose  flux  is  zero. 

When  the  normal  flux  of  species  three  is  constant  every¬ 
where  along  the  metal-electrolyte  interface, 


Jt  tu>  t  n 

3  —  ,  a  constant 

and  equations  (146)  and  (147)  nay  be  written 


(148) 


U(X)  - 

*  +  2D;(C°  +  Cj)^a~  •  y  j 

(U9) 

qw  = 

t1'1  r  'Mi 

_ j  a  -r- 

lyci-  cnL  -  '  •* 

(150) 

When  the  normal  flux  of  species  three  is  constant,  the 
current  density  normal  to  the  metal-electrolyte  interface, 
In,  is  also  constant  and  may  be  written  as 


-118- 


In  =  (151 -a) 

The  current  density  in  the  longitudinal  direction,  Ij ,  is 

I|  =  [  i  +lb)]  (151-h) 

Ij  is  a  linear  function  of  distance  from  the  crack  tip. 

When  the  crack  aspect  ratio  is  large,  the  current  density 
at  the  mouth  of  the  crack  may  be  one  or  more  orders  of  mag¬ 
nitude  greater  than  1^. 

From  equations  (145)  and  (149) >  the  steady  state  value 
of  the  potential  difference  between  the  tip  and  the  mouth  of 
the  crack,  ,  is 


A<t>» 


4LI-r(C?  + 


r 

C^U1 


f? 

+  -  * 


(152) 


For  comparison,  the  initial  value  of  the  potential 
difference,  a60,  is 


Tf  :i  uv  +  c^iv  L«  v;  (153) 

where,  in  writing  equation  (153)  it  is  assumed  that  the 
initial  composition  of  the  electrolyte  is  everywhere  con¬ 
stant  and  equal  to  the  composition  of  the  bulk  electrolyte. 

For  simplicity  and  clarity,  we  will  focus  attention  on 
the  special  case  when  is  constant.  Also,  for  the  example 
to  be  presented,  the  transport  properties,  crack  dimensions, 
and  boundary  conditions  will  be  the  same  as  those  used  in  the 
electrochemical  transport  problem  solved  numerically.  The 
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normal  current  density,  I  ,  will  be  treated  as  the  indepen¬ 
dent  variable.  This  specialization  permits  a  direct  compar¬ 
ison  with  the  numerical  problem  and  provides  a  feel  for  mag¬ 
nitudes  not  conveyed  when  results  are  presented  in  dimen¬ 
sionless  form. 

Thus,  the  parameters  of  the  solution  will  have  the 
following  values  unless  otherwise  stated. 


=  0.200  cm 
=  0.020  cm 

=  0.6  x  IQ'3  mole/cm^ 

=  1.0  x  10“10  mole/cm^ 
=  1 .33  x  10  cm  /sec 
=  2.00  x  10  cm  /sec 
=  9.00  x  1C'5  cm2/sec 


(154) 


For  the  numerical  problem,  the  normal  current  density, 

I  ,  was  314  microamps  per  square  centimeter  (J^  was  3.25  x 
10~9  moles  per  square  centimeter  per  tecond).  When  In  is 
assigned  this  value  and  the  other  parameters  are  evaluated 
as  shown  above,  the  agreement  between  this  solution  and  the 
numerical  solution  is  excellent.  When  the  appropriate  quan¬ 
tities  are  plotted  in  Figures  30,  32,  33  >  and  34,  the  points 
are  coincident  with  those  shown  for  the  numerical  results, 
with  one  exception.  The  current  density  versus  distance 
curve  does  not  show  the  deviation  from  linearity  exhibited 
by  the  numerical  results. 

The  behavior  of  this  system  over  a  broad  range  of  cur¬ 
rent  densities  is  summarized  in  Figures  38,  39 >  and  40.  In 
Figure  38,  the  crack-tip  pH  (defined  here  as  the  negative 
common  logarithm  of  C^,  evaluated  at  x  equn]  to  zero)  is 
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shown  as  a  function  of  the  normal  current  density.  The  con¬ 
centration  of  species  three  is  calculated  from  equations 
(144)  and  (149)  using  the  values  of  the  parameters  given  in 
equations  (154).  The  point  on  the  curve  corresponding  to  the 
results  of  the  numerical  problem  is  shown  by  a  solid  circle. 
The  pH  which  would  exist  in  the  absence  of  a  potential  grad¬ 
ient  (transport  by  simple  diffusion)  is  shown  by  a  broken 
line. 

As  the  current  density  increases ,  the  concentration  in¬ 
creases  and  the  crack-tip  pH  decreases.  When  the  current 
density  is  small ,  the  potential  gradient  is  small  and  the  pH 
is  essentially  the  same  as  would  be  calculated  ignoring  the 
effects  of  migration.  Here,  the  two  curves  are  coincident. 

As  the  concentration  approaches  the  value  of  the  salt  con¬ 
centration  in  the  bulk  electrolyte  (pH--»0),  the  potential 
gradient  begins  to  affect  the  solution  and  the  curves  sepa¬ 
rate.  At  high  current  densities,  the  curves  become  parallel. 
At  any  given  current  density,  the  diffusion  curve  indicates 
a  pH  value  0.301  less  than  the  true  value.  Thus,  the  effect 
of  migration  at  high  current  densities  is  to  reduce  the  pro¬ 
duct  ion  concentration  to  a  value  one-half  as  large  as  that 
predicted  by  diffusion  alone. 

In  Figure  39,  the  crack  tip  concentrations  of  all  three 
ions  are  shown  as  functions  of  the  normal  current  density. 

The  concentrations  are  calculated  from  equations  (142), 

(143),  (144),  and  (149)  using  the  values  of  the  parameters 
in  equations  (154).  The  solid  points  indicate  the  results 
of  the  numerical  problem. 

When  the  current  density  is  small,  the  concentration 
of  the  product  ion  is  small  compared  to  that  of  the  positive 
salt  ion.  In  this  range,  the  concentrations  of  the  salt 
ions  are  insensitive  to  the  applied  current  or  the  presence 
of  the  product  ion.  The  produe"  ion  behaves  as  though  it 
carried  no  charge.  As  the  cui  ••••■.  c  ruity  increases,  the 
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concentration  of  the  product  ion  increases.  When  it  reaches 
a  value  close  to  that  of  the  salt  in  the  bulk  electrolyte, 
there  are  sensible  changes  in  the  concentrations  of  the  salt 
ions.  The  electrolyte  within  the  crack  becomes  enriched  in 
the  negative  salt  ion  and  depleted  in  the  positive  salt  ion. 
This  trend  continues  with  further  increases  in  the  current 
density.  At  very  large  currents,  the  crack  becomes  severely 
depleted  in  the  positive  salt  ion.  In  this  range,  the  be¬ 
havior  is,  for  all  intents  and  purposes,  the  same  as  in  a 
binary  electrolyte. 

The  behavior  of  the  electrostatic  potential  is  illus¬ 
trated  in  Figure  40.  The  potential  difference  between  the 
tip  and  the  mouth  of  the  crack  is  shown  as  a  function  of  the 
applied  current  density.  Two  values  are  shown.  The  upper 
curve  is  the  initial  value  of  the  potential  difference.  It 
was  constructed  using  equation  (153)  and  the  data  in  equa¬ 
tions  (154).  The  lower  curve  is  the  steady  state  value.  It 
was  obtained  from  equation  (152),  again  using  the  data  in 
equations  (154).  The  solid  points  indicate  the  results  of 
the  numerical  problem. 

When  the  normal  current  density  is  small,  both  values 
of  the  potential  are  proportional  to  the  applied  current, 
but  the  steady  state  value  is  about  one-fifth  of  the  initial 
value.  As  the  current  density  increases,  the  initial  value 
remains  proportional  to  the  current.  In  contrast,  the 
steady  state  value  is  proportional  to  the  current  only  with¬ 
in  a  limited  range.  Beyond  this  range,  the  slope  of  the 
potential-current  curve  decreases  with  increasing  current. 

At  very  large  current  densities  the  ratio  of  the  steady 
state  potential  to  the  initial  potential  approaches  zero. 

By  comparison  with  Figure  39,  it  may  be  seen  that  the  cur¬ 
rent  density  corresponding  to  the  upper  limit  of  the  propor¬ 
tional  behavior  is  the  same  as  that  required  to  raise  the 
product  ion  concentration  to  the  level  of  the  salt  concen- 
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tration  in  the  bulk  electrolyte. 

The  example  above  illustrates  most  of  the  important 
features  of  electrochemical  transport.  It  shows  the  entire 
range  of  behavior  of  the  product  ion,  from  simple  diffusion 
at  low  current  densities  to  behavior  as  one  of  the  ions  in 
a  binary  electrolyte  at  high  current  densities.  It  shows 
the  effects  of  the  electroneutrality  requirement  by  the 
changes  in  the  concentrations  of  the  salt  ions  with  changes 
in  the  applied  current.  Most  importantly,  it  shows  how  the 
electrostatic  potential  changes  with  changes  in  the  applied 
current.  Being  a  specific  example,  however,  it  illustrates 
these  effects  in  relation  to  a  specific  set  of  boundary  con- 
ditions.  It  does  not  indicate  the  range  over  which  the 
effects  may  be  observed,  nor  does  it  provide  a  truly  satis¬ 
factory  basis  for  understanding  them.  For  this,  it  is 
necessary  to  return  to  the  solution  as  expressed  by  equa¬ 
tions  (142  through  145)  and  equation  (149). 

Consider  first  the  behavior  of  the  product  ion  when  the 
normal  current  density  is  positive.  Over  most  of  the  range 
of  positive  current  densities,  the  product  ion  behaves  in  one 
of  two  ways.  Depending  on  the  composition  of  the  bulk  elec¬ 
trolyte,  it  may  behave  at  low  current  densities  as  an  un¬ 
charged  species,  obeying  the  differential  equation  for  trans¬ 
port  by  ordinary  diffusion.  At  sufficiently  large  current 
densities,  it  always  behaves  as  one  species  in  a  binary  elec¬ 
trolyte.  (This  may  be  shown  to  hold  quite  generally,  with¬ 
out  regard  for  the  number  or  charges  of  the  other  charged 
species,  as  long  as  each  of  these  species  has  a  zero  normal 
flux.)  The  transition  from  one  behavior  to  the  other  typi¬ 
cally  occurs  over  a  relatively  narrow  range  of  currei. t  den¬ 
sities. 

All  of  this  is  readily  apparent  when  equation  (144)  is 
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written  in  the  form 


AC, 


Cj(o'.  -  C%  r-,  $>Cari  ~  i;\ (bC,)  ; 


(155) 


where 

(156-a) 


(156-b) 

Comparison  of  equations  M56-a)  and  (110-a)  shows  that  bC^ 
is  just  the  concentration  difference  between  the  tip  and  the 
mouth  of  the  crack  that  exists  when  the  product  species  car¬ 
ries  no  charge.  When  f(oC^)  is  zero,  the  product  ion  behaves 
as  an  uncharged  species.  When  it  is  large,  the  product  ion 
behaves  as  an  ion  in  a  binary  electrolyte. 

Fr.om  equation  (156-b),  the  function  f(^C^)  is  bounded  by 
the  limits 
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From  these  bounds  and  equation  (155) ,  it  is  clear  that  the 


For  rotational  convenience,  equation  (144)  has  been  eval¬ 
uated  at  x  equal  to  zero,  the  crack  tip.  For  a  general 
value  of  x,  it  is  necessary  only  to  replace  the  rigbt-hand 
side  of  equation  (156-a)  by  the  right-hand  side  of  equation 
(108). 
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actual  concentration  difference  between  the  tip  and  the 
mouth  of  the  crack  is  never  less  than  one-half,  nor  greater 
than  one,  times  the  concentration  difference  corresponding 
to  transport  by  difi  3  ion  only.  It  is  also  apparent  that 
the  product  ion  can  actually  behave  as  an  uncharged  species 
only  when  its  concentration  in  the  bulk  electrolyte  is  small 
with  respect  to  the  concentration  of  the  positive  salt  ion. 

The  transition  from  one  mode  of  behavior  to  the  other 
occurs  within  a  relatively  narrow  range.  Suppose,  for  ex¬ 
ample,  that  is  much  less  than  C°.  In  this  case,  the 
error  introduced  by  assuming  that  transport  is  by  diffusion 
is  less  than  ten  percent  when  bC^  is  less  than  (^C°).  On 
the  other  hand,  treating  the  electrolyte  as  a  binary  intro¬ 
duces  less  than  ten-percent  error  when  Sc^  is  greater  than 
(160°).  So,  the  electrolyte  must  actually  be  treated  as  a 
ternary  system  only  over  a  range  of  about  two  orders  of 
magnitude . 

The  behavior  of  the  salt  ions  follows  that  of  the  pro¬ 
duct  ion.  When  f(oC^)  is  small  with  respect  to  one,  their 
concentrations  are  essentially  the  same  as  those  in  the  bvlk 
electrolyte.  When  f(Sc^)  is  very  close  to  one,  the  concen¬ 
tration  of  the  positive  salt  ion  is  near  zero,  while  the  con¬ 
centration  of  the  negative  salt  ion  is  approximately  equal  to 
that  of  the  product  ion.  Clearly,  then,  the  value  of  the 
function  f(&C^)  completely  characterizes  the  behavior  of  the 
system  when  the  current  density  is  positive. 

When  the  current  density  is  negative,  the  behavior  is 
different  one  important  respect.  Since  negative  concen¬ 
trations  are  not  physically  allowable,  there  is  a  lower  limit 
to  the  allowable  values  of  the  normal  current  density.  In 
the  terminology  of  equations  (155)  and  (156),  is  bounded 
from  below.  The  value  of  this  bound  is  obtained  by  netting 
03(0)  equal  to  zero  in  equation  (155).  Thus 
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Prom  equations  ( 1 56— "b)  and  (158),  it  follows  that  for 
negative  values  of  the  normal  current  density 
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Comparison  of  equations  (157)  and  (159)  shews  that  the  func¬ 
tion  f(bC^)  is  always  non-negative  and  is  bounded  by  the 
limits  zero  and  one.  It  provides  the  same  characterization 
of  the  system  and  is  amenable  to  the  same  interpretation  at 
negative  currents  as  at  positive  current  .ensities. 

The  normal  current  density  corresponding  to  the  limit¬ 
ing  value  of  6c5  in  equation  (158)  is  called  the  limiting 
diffusion  current  density  or,  simply,  the  diffusion  current 
density.  It  may  be  determined  from  equations  (151-a), 
(156-a),  and  (158).  The  values  of  the  concentrations  and  the 
electrostatic  potential  at  the  crack  tip  may  be  obtained  from 
equations  (142  through  145)  and  equation  (149),  again,  using 
the  limiting  value  of  Sc^  from  equation  (158).  The  results 
are 
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The  value  of  the  diffusion  current  density  in  equation 
(160-a)  is  consistent  with  simple  dimensional  arguments.  It 
increases  as  the  crack  length  decreases  and  as  the  ion  con¬ 
centration  in  the  bulk  electrolyte  increases.  The-  concentra¬ 
tions  of  the  salt  ions  are  equal  and  approach  zero  as  the 
salt  concentration  in  the  bulk  electrolyte  approaches  zero. 
The  value  of  the  electrostatic  potential,  which  is  equal  to 
the  potential  difference  between  the  tip  and  the  mouth  of 
the  crack,  approaches  zero  as  the  concentration  of  the  pro¬ 
duct  ion  in  the  bulk  electrolyte  approaches  zero.  This  is 
not  surprising  since  the  diffusion  current  density  approaches 
zero  at  the  same  time.  As  the  concentration  of  the  positive 
salt  ion  approaches  zero,  the  electrostatic  potential  at  the 
crack  tip  approaches  minus  infinity.  The  diffusion  current 
density  then  approaches  a  value  twice  as  large  as  the  one 
corresponding  to  transport  by  simple  diffusion. 

It  shoved  be  emphasized  that  the  potential  difference 
between  the  tip  and  the  mouth  of  the  crack  is  unbounded  only 
in  the  case  of  a  binary  electrolyte.  The  presence  in  the 
bulk  electrolyte  of  a  third  ion,  having  the  same  charge  as 
the  product  ion,  eliminates  the  singularity  at  the  crack  tip. 
However,  this  ion  must  be  present  in  relatively  large  concen¬ 
trations  to  have  a  practical  effect.  For  example,  when  the 
concentration  of  the  positive  salt  ion,  C°,  in  equation 
(106-e)  is  equal  to  that  of  the  product  ion,  Cy  the  poten¬ 
tial  difference  corresponding  to  the  limiting  diffusion  cur¬ 
rent  density  is  nine  millivolts.  When  C°  is  equal  to  ICT'rC^, 
the  potential  difference  is  59  millivolts.  When  C°  is  eq  ral 
to  the  potential  difference  is  118  millivolts.  These 

values  should  be  compared  to  eighteen  millivolts,  the  value 


' -t --fr ••r*— M  ,-r^v._  .. 


1 

I 

0 

0 


u 

u 


u 

y 

i 

i 


-127- 

of  the  potential  difference  when  the  normal  current  density 
is  of  the  same  magnitude,  but  of  the  opposite  sign  as  the 
diffusion  current  density  in  the  binary  system. 

Prom  the  definition  of  conductivity  following  equation 
(17),  it  is  clear  that  small  changes  in  composition  produce 
small  changes  in  conductivity.  The  large  potential  changes 
illustrated  above  suggest,  therefore,  that  the  steady  state 
potential  does  not  satisfy  Ohm's  law  at  current  densities 
close  to  the  diffusion  current  density. 

The  same  conclusion,  over  a  broad  range  of  current 
densities,  is  suggested  by  the  results  shown  in  Figure  40. 

The  steady-state  potential  difference  was  obtained  from  the 
solution  of  the  governing  differential  equations,  assuming 
that  the  time  rates-of-change  of  the  ion  concentrations  were 
small  with  respect  to  the  other  terms  in  the  equations.  The 
initial  value  was  obtained  from  the  same  set  of  equations, 
assuming  this  time  that  the  electrolyte  was  homogeneous.  As 
indicated  above,  this  is  equivalent  to  the  requirement  that 
the  current  and  potential  satisfy  Ohm's  law.  While  both 
assumptions  are  valid  in  their  stated  time  domains,  the  Ohm’s 
law  statement  does  not  appear  to  be  valid  at  steady  state. 

Since  Ohm's  law  is  often  used  in  the  solution  of  elec¬ 
trochemical  transport  problems,  it  is  worth  considering 
whether  its  use  is  ever  justified.  First,  however,  it  will 
be  convenient  to  rewrite  equations  (152)  and  (153)  in  terms 
of  the  parameter  SC^  defined  by  equation  ( 1 56-a) .  Prom 
equations  (152)  and  (156-a),  the  steady-state  potential  dif¬ 
ference  between  the  tip  and  the  mouth  of  the  crack  is 

=  ln[t  +  ^%)]  (,61) 

Similarly,  from  equations  (153)  and  (156-a),  the  initial 
value  is 
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The  ratio  of  the  two  values  provides  a  simple  measure 
of  the  agreement  between  Ohm's  law  and  the  actual  steady 
state  potential.  For  simplicity,  we  will  consider  only  the 
case  when  is  small  with  respect  to  C°.  In  this  case,  the 
hulk  electrolyte  is  for  all  intents  and  purposes  a  binary 
salt  solution. 

When  the  current  density  is  small,  the  logarithm  in 
equation  (161)  may  be  approximated  by  its  second  term,  and 
the  ratio  of  the  potential  differences  is 
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The  ratio  is  one  when  the  diffusion  coefficient  of  the  pro¬ 
duct  ion  is  equal  to  the  arithmetic  average  of  the  diffusion 
coefficients  of  the  salt  ions.  So,  Ohm's  law  may  yield  a 
satisfactory  approximation  to  the  steady  state  potential 
when  the  diffusion  coefficients  are  all  about  the  same.  Al¬ 
though  this  condition  was  not  satisfied  in  the  example  prob¬ 
lem,  it  often  is  satisfied  for  problems  not  involving  either 
the  hydrogen  or  the  hydroxyl  ions. 

At  high  current  densities,  the  approximation  based  on 
Ohm's  law  is  never  correct.  It  consistently  overestimates 
the  true  steady  state  potential  by  a  significant  margin. 

This  follows  directly  from  equations  (161)  and  (162).  The 
steady  state  potential  increases  as  the  logarithm  of  the 
concentration  Ohm's  law  indicates  that  it  should  in¬ 

crease  linearly. 

While  Ohm's  law  may  yield  incorrect  estimates  of  the 
electrostatic  potential,  the  major  objection  to  its  use  is 


philosophical  rather  than  practical.  Its  use  carries  the 
implication  that  migration  is  the  primary  mode  of  current 
transport  and  that  compositional  differences  between  the 
crack  and  the  bulk  electrolyte  are  negligible.  In  fact, 
one-half  or  more  of  the  total  current  is  transported  by  dif¬ 
fusion,  and,  as  has  been  shown,  the  compositional  differ¬ 
ences  this  requires  are  often  substantial. 

When  the  ion  with  the  non-zero  flux  has  a  negative 
charge,  it  is  only  necessary  to  replace  u  in  equation  (145) 
by  (1/u)  and  assign  the  opposite  signs  to  the  charge  num¬ 
bers  in  equation  (140).  This  changes  the  sign  of  the  re¬ 
lation  between  the  normal  flux  and  the  current  density  but 
has  no  other  effect.  Therefore,  the  results  presented  here 
apply  to  any  1-1-1  electrolyte  for  which  the  normal  fluxes 
of  two  species  are  zero.  They  may  be  conveniently  summa¬ 
rized  as  follows. 

1 )  When  it  is  assumed  that  transport  is  by  simple  diffu¬ 
sion,  the  resulting  estimate  of  the  product  ion  concen¬ 
tration  is  never  in  error  by  more  than  a  factor  of  two. 

2)  An  estimate  of  the  error  in  the  above  approximation  may 
be  obtained  from  the  ratio,  R  ,  of  the  calculated  con- 
centration  to  the  total  salt  concentration  in  the  bulk 
electrolyte.  Wnen  the  ratio  is  less  than  one-half,  the 
error  is  less  than  ten  percent.  When  the  ratio  is 
greater  than  twenty,  the  error  is  between  1.9  and  2.0. 

In  the  latter  case,  the  electrolyte  essentially  behaves 
as  a  binary  system,  so  the  binary  approximation  will 
often  be  more  useful  in  this  range  than  the  simple  dif¬ 
fusion  analysis. 

5)  The  concentrations  of  the  salt  ions  follow  that  of  the 
product  ion.  The  crack  is  depleted  in  the  stilt  ion 
carrying  the  same  charge  as  the  product  ion,  and  concen- 
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trated  in  the  salt  ion  carrying  the  opposite  charge. 

4)  The  potential  is  described  by  Ohm's  law  at  times  small 

p 

with  respect  to  the  diffusion  time  (a  /D^).  The  steady 
state  potential  may  be  approximated  using  equation  (161) 
with  bCj  replaced  by  its  approximate  value  from  the 
diffusion  analysis.  If  the  concentration  estimate  is 
corrected  as  indicated  above,  the  expression  for  the 
potential  is  exact. 


These  res- tits  are  consistent  with  those  obtained  numer¬ 
ically.  They  indicate  that  the  diffusion  equation  provides 
a  good  estimate  of  the  concentration  of  a  product  species. 
When  the  requirements  of  electroneutrality  are  considered, 
satisfactory  approximations  to  the  salt  ion  concentrations 
may  also  be  obtained. 

Recapitulation 

When  the  separate  observations  and  findings  of  the  last 
two  sections  are  viewed  as  a  whole,  they  form  a  coherent  and 
fairly  complete  picture  of  corrosion  and  transport  in  crack¬ 
like  regions.  They  indicate  that  the  average  concentration 
of  an  ion  which  enters  into  a  surface  reaction  is,  to  a  good 
first  approximation,  given  by 
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Thus,  the  factors  which  have  a  primary  influence  on  the  com¬ 
position  within  a  crack-like  region  are  the  average  values 
of  the  normal  ion  fluxes,  the  crack  aspect  ratio,  and  the 
crack  length.  Other  factcrs,  such  as  the  shape  of  the  crack, 
the  exact  form  of  the  boundary  conditions,  and  the  transport 
mode,  can  and  do  affect  the  composition.  In  general,  however, 
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theee  factors  are  of  lesser  importance. 

Operationally,  the  results  suggest  that  first  order 
approximations  to  the  ion  concentrations  may  he  obtained 
using  a  method  similar  to  that  described  by  equations  (132). 
The  accuracy  of  the  resulting  estimates  may  then  be  deter¬ 
mined  by  comparing  them  with  the  ionic  strength  of  the  bulk 
electrolyte.  For  moat  problems  it  will  be  found  that  the 
electrolyte  can  be  treated  either  as  an  ordinary  solution  or 
as  a  binary  electrolyte.  That  is,  either  the  dissolved 
species  may  be  treated  as  having  zero  charge  and  the  grad¬ 
ient  of  the  electrostatic  potential  ignored,  or  they  may  be 
treated  as  charged  species  moving  in  the  electrostatic  field 
created  by  the  transport  of  the  two  principal  ions. 

Clearly,  situations  will  arise  for  which  neither  of 
these  approximations  is  satisfactory.  In  such  cases,  it 
will  be  necessary  to  use  the  more  sophisticated  solution 
techniques  developed  in  this  work.  When  two-dimensional 
transport  effects  are  known  or  thought  to  be  important  the 
numerical  method  iB  particularly  well-suited.  When  the 
boundary  conditions  are  simple  and  the  time-rate-of-change 
of  the  electrolyte  composition  is  small,  the  steady  state 
form  of  the  one-dimensional  transport  analysis  is  preferred. 
For  problems  involving  time-dependence  or  complex  bounaary 
specifications,  the  numerical  technique  is  the  best  of  the 
existing  methods.  However,  a  numerical  method  based  on  the 
one-dimensional  transport  analysis  would  be  simpler  and  more 
efficient  in  the  long  run. 


***» 
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VII.  CONCLUDING  REMARKS 


The  work  presented  here  was  originally  aimed  at  devel¬ 
oping  an  analytic  framework  for  the  study  of  stress  corro¬ 
sion  cracking.  As  the  work  progressed,  it  became  cleer  that 
it  would  be  necessary  first  to  determine  what  factors  most 
strongly  affect  corrosion  and  transport  processes  in  cracks. 
Only  then  would  it  be  reasonable  to  undertake  the  modeling 
of  such  a  complex  process  as  crack  growth.  While  the  results 
do  not  go  beyond  that  first  stage,  they  do  provide  a  basis 
for  a  discussion  of  some  general  aspects  of  stress  corrosion 
cracking. 

One  of  the  longstanding  rules  of  thumb  in  stress  corro¬ 
sion  cracking  is  'metals  that  corrode  don't  stress  corrode'. 
Like  other  such  rales,  it  is  not  necessarily  true,  but  it 
points  out  that,  by  and  large,  alloys  susceptible  to  crack¬ 
ing  in  a  particular  environment  do  not  actively  corrode  in 
that  environment.  They  typically  have  corrosion  potentials 

in  the  passive  region,  a  range  of  potentials  within  which  a 

* 

metal-oxide  is  stable  in  contact  with  water.  When  there  is 
an  adequate  oxygen  supply,  the  oxide  usually  forms  by  reac¬ 
tion  of  the  metal  with  dissolved  oxygen.  When  the  oxygen 
supply  is  limited,  as  is  probable  in  a  crack  or  crevice,  the 
oxide  forms  either  by  hydrolysis  of  metal  ions  or  by  direct 
reaction  of  the  metal  with  water.  Regardless  of  the  reaction 
path,  the  formation  of  a  metal-oxide  in  the  absence  of  free 
oxygen  is  accompanied  by  the  release  of  hydrogen  ions. 

If  the  electrolyte  within  a  crack  is  depleted  in  oxygen, 
the  pH  should  be  in  the  acid  range.  Evidence  presented  in 
two  recent  papers,  one  by  B.  F.  Brown,  C.  T.  Fujii,  and  E.  P. 
Dahlberg  [3]  and  the  other  by  J.  A.  Smith,  M.  H.  Peterson, 
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and  B.  F.  Brown  [24],  indicates  that  this  is,  in  fact,  the 
case.  Their  data  show  that  for  cracks  in  a  variety  of  ma¬ 
terials  the  pH  is  in  the  range  between  about  1.5  and  3*5. 

In  general,  the  data  indicate  that  the  greater  the  oxide 
stability,  the  lower  the  measured  pH. 

These  pH  measurements  may  be  used  with  the  results  of 
the  last  section  to  obtain  some  general  estimates  of  the 
conditions  in  stress  corrosion  cracks.  For  example,  equa¬ 
tion  (164)  can  be  used  to  estimate  the  normal  current  den¬ 
sity  at  the  metal-electrolyte  interface  due  to  oxide  film 
growth.  Consider  a  crack  1 -millimeter  long  by  0.01-milli- 
meter  wide  and  suppose  that  the  bulk  electrolyte  is  a  neu¬ 
tral,  0.6-mole  per  liter,  sodium-chloride  solution.  Us^ug 
-5 

9  x  lO  centimeters  squared  per  second  as  the  value  of  the 
hydrogen  ion  diffusion  coefficient,  the  average  value  of  the 
normal  current  density  is  in  the  range  3  x  10-7  to  3  x  10"5 
amperes  per  square  centimeter. 

If  the  rate  of  film  growth  changes  significantly  from 
point  to  point  along  the  crack  wall,  the  average  normal 
current  density  is  not  a  good  indicator  of  the  local  reac¬ 
tion  rate.  If,  as  appears  likely ,  the  rate  at  the  crack  tip 
is  much- greater  than  at  points  more  distant,  the  current  den¬ 
sity  in  that  region  will  be  about  two  orders  of  magnitude 
greater  than  the  average  calculated  above.  Whatever  the 

case,  the  current  density  at  the  mouth  of  the  crack  is  in 
—5  —3 

the  range  3  x  10  to  3  r  10  ;  amperes  per  square  centimeter. 

There  is  little  published  data  regarding  the  concentra¬ 
tions  of  ions  other  than  the  hydrogen  ion.  Measurements  re¬ 
ported  in  Reference  3  indicate  the  presence  of  metal  ions 
but  provide  no  estimates  of  their  concentrations.  From  the 
reported  pH  values  and  the  solubility  products  of  the  various 
metal  hydroxides,  however,  it  appears  that  the  metal  don  con¬ 
centrations  are  small  and  that  the  hydrogen  ion  is  the 
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principal  product  ion. 

As  shown  in  the  last  section,  the  importance  of  migra¬ 
tion  may  be  estimated  from  the  re  no  of  the  product  ion  con¬ 
centration  in  the  crack  to  the  ionic  strength  of  the  bulk 
electrolyte.  From  the  measurements  reported  in  Reference  3, 
the  minimum  pH  was  about  1.5,  corresponding  to  a  hydrogen 
ion  concentration  of  about  0.03  moles  per  liter.  The  ionic 
strength  of  the  electrolyte  (0.6-mole  per  liter,  sodium- 
chloride  solution)  was  about  1.2  moles  per  liter.  The  ratio 
of  the  two  is  forty.  Thus,  the  effects  of  migration  and  the 
gradient  of  the  electrostatic  potential  should  be  negligible 
over  most  of  the  crack  length.  It  is  interesting  to  note 
that  in  T.  R.  Beck's  original  'MET'  analysis  [lO],  the  same 
result  was  obtained.  A  large  potential  gradient  in  the  re¬ 
gion  close  to  the  crack  tip  was  later  generated  by  introduc¬ 
ing  a  diffusion-limited  chloride  ion  flux. 

The  results  presented  in  the  last  two  sections  are  based 
on  a  planar  analysis  and  are  strictly  valid  only  when  trans¬ 
port  in  the  specimen  thickness  direction  is  negligible.  In 
most  experimental  work  no  attempt  is  made  to  control  trans¬ 
port  in  this  direction.  A  planar  analysis  should  neverthe¬ 
less  yield  results  consistent  with  experiment  when  the  spec¬ 
imen  thickness  is  large  with  respect  to  the  crack  length. 

When  the  thickness  is  small,  the  analysis  will  fail. 

In  practice,  a  planar  analysis  should  be  approximately 
correct  for  points  closer  to  the  crack  tip  than  about  one- 
half  the  specimen  width.  At  points  more  distant,  composition 
variations  in  the  specimen  thickness  direction  will  be  much 
greater  than  those  in  '■'\e  crack  length  direction.  Corrosion 


The  system  examined  in  Reference  24  is  an  exception  to 
this  rule.  The  metal  (AISI  4340  steel)  may  crack  at  active 
corrosion  potentials.  Although  the  crack  is  acidified  by 
what  appears  to  be  a  hydrolysis  reaction,  the  hydrolyzed 
species  is  soluble. 
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processes  on  the  crack  faces  will  then  be  little  different 
from  those  on  the  sides  of  the  specimen.  Therefore,  finite- 
width  specimens  should  behave  as  if  the  crack  length  were 
constant  whenever  the  crack  length  is  greater  than  about 
one-half  the  specimen  width. 

This  view  is  consistent  with  the  outstanding  success 
obtained  in  correlating  susceptibility  to  cracking  and  crack 
growth  rates  with  applied  stress  intensity.  For  pre-cracked 
specimens  of  a  fixed  alloy  composition,  tested  under  constant 
conditions  in  a  given  electrolyte,  there  appears  to  be  a 
critical  stress  intensity,  ^IgCC,  below  which  cracks  do  not 
propagate.  This  critical  value  is  independent  of  the  crack 
length  or  the  specimen  configuration  [25].  In  the  last  sec¬ 
tion  it  was  sho.  i  that  convective  transport  is  signiiicant 
only  at  very  high  loading  rates.  In  the  usual  range  of  load¬ 
ing  rates,  the  major  effect  of  the  applied  stress  on  ion 
transport  is  its  effect  on  the  crack  aspect  ratio.  For  a 
pre-cracked  specimen  such  as  that  shown  in  Figure  24,  the 
crack  aspect  ratio  depends  on  its  value  under  no  load  and 
the  stress  intensity.  By  using  a  fatigue  crack  as  a  starter 
notch,  the  crack  aspect  ratio  at  zero  load  is  effectively 
standardized,  and  the  aspect  ratio  becomes  a  function  of  the 
stress  intensity  only.  Therefore,  if  the  crack  length  to 
specimen  width  ratio  is  in  the  range  where  the  'effective' 
crack  length  remains  constant,  the  effect  of  the  stress  field 
on  ion  transport  is  determined  by  the  applied  stress  inten¬ 
sity  and  the  variables  describing  the  crack  geometry  do  not 
behave  as  independent  variables. 

These  remarks  are  not  meant  to  imply  that  the  stress 
field  has  no  influence  other  than  its  influence  on  ion  trans¬ 
port.  The  applied  stress  intensity  has  a  strong  influence  on 
fracture  surface  morphology  and  may  affect  the  kinetics  of 
charge-transfer  reactions.  The  fact  that  the  stress  inten¬ 
sity  influences  all  these  factors  is  one  of  the  more  serious 
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obstacles  to  determination  of  the  underlying  mechanisms  of 
cracking. 

The  need  for  accurate  kinetic  data  has  been  stressed 
throughout  this  work.  Qualitative  observations  and  order  of 
magnitude  estimates  such  as  those  presented  above  provide 
some  insight  into  the  stress  corrosion  process.  But  these 
are  really  not  adequate  for  understanding  or  controlling 
stress  corrosion  cracking  and,  if  history  is  a  useful  guide, 
they  are  just  as  likely  to  be  wrong  as  right. 

Consider  some  of  the  interpretations  made  of  the  ob¬ 
served  low  values  of  the  crack  pH.  This  observation  has 
been  viewed  as  providing  support  for  at  least  three  differ¬ 
ent  mechanisms  for  crack  growth.  While  it  would  be  unfair 
and  misleading  to  imply  that  these  proposals  are  based  on  no 
other  experimental  data,  a  low  pH  is  a  critical  factor  in 
each  of  them. 

T.  P.  Hoar  [26]  has  argued  that  a  How  pH  inside  the 
crack  is  consistent  with  a  metal-dissolution  mechanism. 

Under  such  circumstances,  the  rate  of  oxiue-film  formation 
should  be  suppressed  and  consequently  the  rate  of  metal  dis¬ 
solution  highly  favored  near  the  crack  tip.  Others  argue 
that  a  low  pH  should  accelerate  the  rote  of  hydrogen  ion 
discharge  on  the  bare  metal  surface  ab  the  crack  tip,  and 
that  fracture  must  th.n  be  the  consequence  of  hydrogen,  em¬ 
brittlement  (see  Reference  27).  Still  others  find  support 
for  a  film  rupture  mechanism.  Here  it  is  argued  that  the 
low  pH  confirms  the  hypothesis  that  oxide  film  growth  is  the 
primary  reaction  occurring  in  cracks.  If  this  is  the  case, 
rupture  of  the  film  under  applied  stress  may  then  cause  the 
crack  to  propagate  a  3hort  distance  into  the  metal  until  it 
again  comes  to  rest  and  the  process  repeats  itself. 

Obviously  a  low  pH  value  neither  confirms  nor  refutes 
any  of  these  theories.  Together  with  other  d' ta,  the  mea¬ 
surements  reported  in  Reference  3  do  imply  that  hydrogen 
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ions  are  generated  either  by  the  hydrolysis  of  metal  ions  or 
by  oxide  film  growth.  Whether  this  reaction  is  even  critical 
to  crack  propagation  remains  to  be  determined.  For  many  sys¬ 
tems,  a  low  crack  pH  is  expected  simply  from  the  exposure  of 
clean  and  highly  reactive  metal  surface  to  an  oxygen  depleted 
electrolyte.  When  this  is  the  case,  a  low  pH  could  be  the 
consequence  of  crack  growth  and  not  the  underlying  cause  of 
it. 

Regardless  of  the  actual  mechanism  of  crack  growth,  it 
is  clear  that  all  require  a  fine  balance  between  the  rates 
of  several  different  reactions.  For  example,  if  crack  prop¬ 
agation  is  due  to  hydrogen  embrittlement  in  a  near  neutral 
electrolyte,  there  must  first  be  a  source  of  hydrogen  ions. 
The  rate  of  generation  of  these  ions  must  be  approximately 
equal  to  their  rate  of  discharge  on  the  metal  surface.  The 
rate  of  hydrogen  absorption  must  be  fast  with  respect  to  the 
rate  of  formation  of  molecular  hydrogen  and  its  rate  of 
transport  away  from  the  crack  tip.  Finally,  the  rate  of 
mefcal  dissolution  must  be  slow  with  respect  to  the  rate  of 
hydrogen  absorption. 

Because  there  is  qualitative  agreement  between  the  ex¬ 
perimental  data  and  several  different  failure  mechanisms,  it 
is  not  possible  to  determine  the  actual  mechanism  in  a  given 
system  by  a  single  simple  experiment.  In  order  to  subject 
any  of  the  mechanistic  theories  of  stress  corrosion  to  a 
critical  test,  at  least  three  elements  will  be  required. 
First,  the  mechanism  must  be  cast  in  a  quantitative  form. 

This  will  probably  require  independent  evaluation  of  quali¬ 
tative  relationships  which  have  not  yet  been  quantified.  In 
the  case  of  hydrogen  embrittlement,  for  example,  the  rela¬ 
tionship  between  Kjq,  the  critical  stress  intensity  for  rapid 
fracture,  and  the  concentration  of  occluded  hydrogen  is  re¬ 
quired.  Next,  accurate  and  independently  measured  electro¬ 
chemical  kinetic  data  must  be  obtained.  Using  this  data, 
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the  analytic  techniques  presented  here  may  be  used  to  char¬ 
acterize  the  environment  and  to  determine  whether  the  actual  || 

reaction  rates  are  consistent  with  the  mechanism  postulated. 

Finally,  the  quantitative  predictions  generated  in  this  way  r 

must  be  subjected  to  rigorous  experimental  verification. 

I 

E 

E 

E 

E 

E 

E 

0 

E 

E 

E 

r 

L 

l 

l 

1 


-139- 


I 

1 

I 

I 
3 
0 
a 
0 
0 

II 
a 
a 

ii 


Conclusions 


1 )  A  method  has  been  developed  for  the  numerical  solution 
of  the  electrochemical  transport  equations  based  on 
dilute  solution  theory  and  subject  to  an  arbitrary  set 
of  charge-transfer  boundary  conditions.  It  has  been 
applied  to  problems  involving  dependence  on  time  and 
two  spatial  variables. 

2)  The  factors  of  primary  .  mportance  in  determining  the 
electrolyte  composition  within  a  crack-like  region  are 
the  average  values  of  the  species  fluxes  normal  to  the 
solid-liquid  interface,  the  crack  length,  and  the  crack 
aspect  ratio.  Other  factors,  such  as  the  crack  shape, 
the  form  of  the  boundary  conditions,  and  the  transport 
mode  can  affect  the  composition.  However,  these  fac¬ 
tors  are  of  lesser  importance. 

3)  When  the  crack  aspect  ratio  is  greater  than  three  and 
the  elapsed  time  is  greater  than  the  diffusion  time 

p 

(a  /D),  the  effects  of  transport  in  the  transverse 
crack  direction  are  negligible. 

4)  The  average  concentration  of  an  ion  within  a  crack-like 
region  is  given  to  a  good  first  approximation  by  equa¬ 
tion  (164) . 
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The  accuracy  of  this  approximation  may  be  determined 
from  the  ratio  of  the  calculated  value  to  the  ionic 
strength  of  the  bulk  electrolyte. 

5)  As  shown  by  equation  (164),  the  electrolyte  composition 
within  a  crack  may  differ  substantially  from  that  in 
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the  bulk  solution. 

6)  An  accurate  representation  of  reaction  kinetics  is 
required  to  determine  the  detailed  behavior  of  an 
electrolyte  within  a  crack-like  region. 

7)  Many  individual  results  have  been  presented  which  can¬ 
not  be  conveniently  summarized  but  have  direct  appli¬ 
cation  to  pitting  corrosion,  crevice  corrosion,  and 
stress  corrosion  cracking. 


GHMSMM 


-141- 

Recommendations  and  Future  Work 

1 )  For  application  to  stress  corrosion  cracking,  pitting 
corrosion  and  crevice  corrosion,  a  method  should  be 
developed  for  the  numerical  solution  of  the  one-dimen¬ 
sional  formulation  of  the  transport  equations  presented 
in  Appendix  C.  This  should  permit  improved  resolution 
in  regions  with  large  gradients  and  result  in  at  least 
an  order  of  magnitude  reduction  in  core  storage  re¬ 
quirements  and  computation  time. 

2)  The  methods  and  results  developed  in  this  study  should 
be  applied  to  the  study  of  stress  corrosion  cracking. 
Particular  emphasis  should  be  placed  on  the  identifi¬ 
cation  of  allowable  reactions  in  the  system  under  in¬ 
vestigation  and  on  a  determination  of  the  kinetics  of 
these  reactions. 

3)  The  many  individual  results  presented  here  provide 
tools  for  the  study  of  such  corrosion  phenomena  as 
pitting,  crevice  corrosion,  and  intergranular  corrosion. 
Application  of  these  results  to  such  corrosion  processes 
should  be  undertaken. 
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APPENDIX  A 


Finite-Difference  Analogue  for  the  Ion  Conservation  Equation 

A  key  element  in  the  numerical  solution  of  the  electro¬ 
chemical  transport  equations  is  the  method  used  to  solve  the 
ion  conservation  equation,  equation  (75)  on  page  38.  In  the 
notation  of  equations  (77)  >  the  equation  may  he  written 


X 


t  2iC; 
dt 


+  BiV-Ci  -  VCi-co^A^I  +  \[lQri 


(A-1) 


In  the  computer  program  developed  for  this  study,  equa¬ 
tion  (A-1)  is  solved  using  an  Alternating-Direction  Implicit 
(ADI)  technique.  The  term  ADI  describes  a  general  concept 
or  approach  rather  than  a  particular  algorithm.  It  may  be 
applied  to  any  algorithm  using  a  computational  sequence  in 
which  the  equations  are  first  made  implicit  in  one  coordin¬ 
ate  direction  and  then  in  another.  Such  methods  have  the 
advantages  of  unconditional  stability  and  computational 
speed.  This  appendix  is  included  to  show  how  this  concept 
has  been  applied  to  the  solution  of  equation  (A-1). 

Let 

I  Li)  (A-2) 


where  f  is  any  function  of  the  coordinates  tt  %  and  t.  Then 
by  analogy  with  differential  operators,  we  may  define  the 
following  difference  operators 
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A^(j,k,JO  =  2^[ Vi+l‘k.^  " 


(A-3) 


•  eA7L^j,k+i,i)  _  ^(j,k.-ij.;] 


4$^)  ~  *  r4<j>U)] 

AtV,i ,M)  =--  *  \i\MM)] 


(A-4) 


A\S(j,k,tt  =  +  $(j-i,M)]  (a-5) 


(A— 6 ) 


(A-7) 


The  central-difference  operators  defined  by  equations 
(A-3  through  A-6)  are  second-order  correct.  That  is,  the 
error  introduced  by  using  the  difference  operator  to  approx¬ 
imate  the  analogous  differential  operator  is  proportional  to 
(or  A7?).  The  backward-difference  operator  defined  by 
equation  (A-7)  is  first-order  correct. 


The  finite- difference  analogue  to  equation  (A-1 )  is 
written  as  two  equations.  In  the  first,  derivatives  of 
concentration  in  the  ^-direction  aro  approximated  at 
while  those  in  the  7-direction  are  approximated  at  S.t\t .  In 
the  second,  the  sequence  is  reversed  and  the  calculation  is 
advanced  another  half  time-step.  Derivatives  of  concentra¬ 
tion  in  the  7-directi.on  are  approximated  at  ({+l)At  and  de¬ 
rivatives  in  the  ^-direction  are  approximated  at  ((+  |;)At . 
There  are  several  ways  to  treat  the  products  of  the  deriva¬ 
tives  of  Ct  and  (|>.  In  the  present  version  of  the  program, 
these  are  approximated  at  (S+£)a£  . 
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The  difference  equations  are 

MV%C(j,ia4)  D^Ctj.k ,X*{)  +  B^C(j,kJ0  (A-8) 

+  2.)  +  ^  G*( j , , A *"2^ 


Note  that  in  writing  equations  (A-8  through  A-11)  the  species 
subscript  i  has  been  omitted  and  the  stream  function  has  been 
written  as  a  time  invariant  quantity.  This  was  done  for  sim¬ 
plicity.  Values  of  the  functions  and  G<  at  the  advanced 
time  are  obtained  by  iteration  as  explained  in  the  text. 

The  net  effect  of  alternately  applying  equations  (A-8) 
and  (A-9)  can  be  seen  from  the  sum  of  the  two  equations. 


LCCj.K.A*!',  -  CCi.k^r,]  -  XH.)  +  {b[k\C[^)+ 


(A-12) 

+  {[L,c:j1k1r  ■  &7c(j,klhi)]^(],^ 


The  term  on  the  left-hand  side  is  just  the  central  dif¬ 
ference  approximation  for  the  partial  time  derivative  of  C 
at  time  ($  +  i)At.  It  is  analogous  to  equations  (A-3)  and 
(A-4)  and  may  be  shown  to  be  second-order  correct.  There¬ 
fore,  equation  (A-12)  is  second-order  correct  in  both  time 
and  space  when  and  0>  are  correctly  evaluated  at  time 
(-£+{)&£.  When  the  functions  q  and  G*  are  evaluated  at  time 
i&t,  the  equation  remains  second-order  correct  in  space  but 
is  reduced  to  first-order  correctness  in  time. 
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APPENDIX  B 
Reaction  Kinetics 

Before  the  techniques  developed  here  can  be  applied  in 
the  study  of  a  particular  system,  one  must  first  decide  what 
reactions  are  relevant  to  the  behavior  of  that  system,  and 
then  obtain  quantitative  expressions  for  the  rates  of  those 
reactions.  For  example,  having  determined  that  a  redox  re¬ 
action  of  the  form  of  equation  (40)  is  important,  one  must 
decide  whether  an  equation  of  the  form  of  equation  (59)  is 
suitable  for  expressing  its  reaction  kinetics.  If  so,  the 
parameters  I*  and  o'  must  be  evaluated.  If  not,  another  ex¬ 
pression  must  be  found  and  its  parameters  evaluated.  The 
same  thing  applies  to  the  other  reactions  in  the  system. 

For  the  work  presented  here,  it  was  decided  to  select 
a  simple  system  of  reactions  which  would  exhibit  a  fairly 
wide  range  of  behavior  without  being  excessively  complicated. 
It  was  also  desired  that  the  system  of  reactions  be  one  for 
which  some,  if  not  all,  of  the  kinetic  constants  could  be 
evaluated  from  published  data. 

A  system  of  reactions  modeled  after  the  oxide-film 
formation  reactions  in  the  titanium-water  system  met  these 
requirements.  A  key  factor  in  this  choice  was  that  the 
oxides  of  titanium  have  an  extremely  low  electronic  conduc¬ 
tivity.  Virtually  all  of  the  current  passing  through  the 
film  must  be  transported  by  ions  or  vacancies  rather  than  by 
electrons.  In  the  limjt  of  zero  electronic  conductivity,  no 
redox  reaction  can  occur  on  an  oxide-covered  metal  surface, 
since  the  electrons  generated  cannot  pass  through  the  film. 

In  this  case,  the  only  allowable  reactions  are  oxide  film 
growth  by  ion  or  vacancy  migration  through  the  film,  and 
ion-exchange  reactions  at  the  oxide-electrolyte  interface. 

By  further  limiting  the  ion-exchange  reactions  to  a  single 


dissolution-precipitation  reaction,  a  relatively  simple 
system  of  reactions  is  obtained  which  nevertheless  retains 
some  physical  relevance. 

The  reactions  considered,  therefore,  may  be  written 
M  +  2H20  =  M02  +  4H+  +  4e“  (B-1) 

M02  +  2H+  ~  M02+  +  H20  (B— 2 ) 


Reaction  (B-1 )  may  be  viewed  as  an  overall  reaction 
consisting  of  three  steps. 


M(metal)  ^  M^+(oxide) 

+  4e“ 

Cii] 

(B-3a) 

migration  of  ions  or  vacancies 

[I,] 

(B-3b) 

through  oxide  film 

2[H20  ^  02“( oxide)  + 

2H+ ( elect )  ] 

[2I2] 

(B-3c) 

Steps  (B-3a)  and  (B-3c)  are  charge-transfer  reactions.  Step 
(B-3b)  is  a  transport  requirement  imposed  by  the  finite 
thickness  of  the  oxide  film.  The  current  density  correspond¬ 
ing  to  each  step  is  shown  in  brackets  for  later  reference. 

If  no  other  reaction  occurs,  stoichiometry  of  the  oxide  re¬ 
quires  that  the  current  density  I2  be  one-half  1^ . 

Reaction  (B-2)  may  be  viewed  as  consisting  of  two 
partial  reactions. 

M02+(oxide)  ^  M02+(elect)  [i^]  (B-4a) 
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O2”(oxide)  +  2H+( elect)  ^  H20  [-Ig]  (B-4b) 


Reactions  (B-4)  are  both  charge -transfer  reactions  occurring 
at  the  oxide -electrolyte  interface.  Reactions  (B-4b)  and 
(B-3c)  are  identical  but  are  written  in  the  opposite  direc¬ 
tions.  Again  the  current  density  corresponding  to  each  re¬ 
action  is  shown  in  brackets.  When  reaction  (B-2)  is  the 
only  reaction,  stoichiometry  of  the  oxide  requires  that  I2 
equal  minus  I-. 

The  overall  system  of  two  reactions  consists  of  three 
partial  reactions.  The  requirement  of  stoichiometry  may  be 
written  as 
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I4  -  Ii  -  U  -  o 


(B-5) 


To  obtain  a  complete  description  of  the  reaction  kinet¬ 
ics  it  is  necessary  to  express  the  current  densities,  1^ , 

I2,  and  1^,  as  functions  of  concentration  and  potential. 
Following  Vetter  [5]  it  is  assumed  that  reaction  (B-3a)  is 
near  equilibrium  and  transient  transport  effects  in  the 
oxide  film  are  ignored.  The  current  density  1^  may  then  be 
expressed  using  the  high  field  conduction  equation.  The 
current  densities  are 


I,  =  i;[2xp[f(<^,-<|ve,-70;  -  cyp[- t(<t>„-<Ve,,-70] 

la  -  i;^::.p[2rt7,]  -  (i^  Cxp[-*(!- *V£]j 


(B— 6 ) 


(B-7) 


(B-8) 
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Equatlons  (B-5  through  B-8)  provide  the  complete  ki¬ 
netic  description  of  reactions  (B-1)  and  (B-2).  Because 
the  solution  of  these  equations  for  ?$  is  lengthy  and  in¬ 
volved  and  a  parametric  analysis  of  the  equations  indicated 
that  is  small  over  a  wide  range  of  pH  and  potential,  it 
was  decided  to  treat  reaction  (B-2)  as  being  close  to  equi¬ 
librium.  Equation  (99)  in  the  text  is  the  result  of  this 
assumption.  A  simple  algorithm  for  obtaining  7*  from  this 
system  of  equations  has  since  been  developed,  so  a  simpli¬ 
fying  assumption  of  the  kind  used  here  need  not  be  invoked 
in  future  work. 

Once  7S  has  been  determined,  the  ion  fluxes  and  oxide- 
film  growth-rate  may  be  expressed  as 


(B-9a) 

JjM?uVri  -  "  l!/E3' 

(B-9b) 

8  =  Xl  '  21»)'4P3. 

(B-9c) 
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APPENDIX  C 

Electrochemical  Transport  in  a  Crack-like  Region 

When  the  aspect  ratio  of  a  notch  or  crack  is  large,  the 
effects  of  transport  in  the  transverse  direction  can  be  ne¬ 
glected  and  the  region  treated  as  a  one-dimensional  contin¬ 
uum.  The  equations  describing  transport  in  such  a  region 
are  developed  in  this  appendix  and  the  solutions  for  two  im¬ 
portant  cases  are  presented.  A  solution  technique  based  on 
developing  an  integral  equation  for  the  electrostatic  poten¬ 
tial  is  also  discussed.  An  example  of  the  use  of  this  tech¬ 
nique  is  presented  in  Appendix  D. 

Consider  a  region  in  the  x-y  plane  bounded  by  the 
x-axis  and  the  line  y  =  f (x)  as  shown  in  Figure  41 .  The 
region  corresponds  to  one-half  the  section  of  a  crack  pro¬ 
duced  by  a  cutting-plane  perpendicular  to  both  the  crack 
plane  and  the  crack  front.  The  x-axis  coincides  with  the 
crack  centerline;  the  line  y  =  f(x)  coincides  with  one  crack 
face.  The  region  is  assumed  to  be  symmetric  about  the  crack 
centerline  so  the  second  crack  face  corresponds  to  the  line 
y  =-f(x).  For  the  control  volume  shown  in  Figure  41,  the 
equation  of  mass  conservation  may  be  written 


J(!y  l^CxUi  -  «-)(* -  bKht  +  Ji 

t  ♦  e..'t 


I*  r;.t 


AMt  = 


[C;|  .  ”  C-!  lAV'TqUl  t 

iMt  I  t 


(c-1) 


Taking  the  limit  of  equation  (C-1)  as  ax  and  At  approach 
zero,  the  differential  form  of  the  mass  conservation  equa¬ 
tion  is 


0 


(C-2) 


where 

&  =  VT+ «w) 

When  the  mouth  of  the  crack  is  located  at  x  equal  to  a,  and 
h  is  the  crack  half -width  measured  at  the  crack  mouth, 

l*  ~  b,  ^  «  0  (C-3a) 

for  a  rectangular  crack  and 

1-1  -  fy  =  Vl  +  (C-3b) 

for  a  pie-shaped  crack.  Note  that  the  fluxes  and  in 
equation  (C-2)  are  flux  densities  and  not  total  fluxe3. 

Also,  the  quantities  appearing  in  the  equation  correspond 
to  average  values  across  the  crack  width. 

To  complete  the  description,  the  flux  equation  and  the 
electroneutrality  equation  are  required. 


r  ^  rs 

•J 


(0-4) 


( 0—5 ) 


When  the  average  solvent  velocity , V,  has 
equations  (C-2,  -4,  and  -5)  represent  the 
of  equations  describing  ma33  transport  in 
region. 


been  determined, 
complete  system 
a  crack-like 
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When  the  line  y  =  f(x)  is  the  generator  of  a  surface  of 
revolution  about  the  x-axis  rather  than  a  sheet,  the  enclosed 
region  is  pit-like  rather  than  crack-like.  Using  the  same 
procedure  as  that  used  to  obtain  equation  (C-2),  the  conser¬ 
vation  equation  for  such  a  region  is 

&(w)  ~  +  $  W  =  o  (o-e) 

When  equation  ( C— 6 )  is  used  instead  of  equation  (C-2),  the 
results  are  modified  in  detail  but  not  in  principle.  There¬ 
fore,  it  will  not  be  considered  further. 

At  steady  state,  integration  of  the  conservation  equa¬ 
tion  is  straightforward.  Since  the  boundary  conditions  typi¬ 
cally  require  specifying  the  flux  at  the  crack  tip  and  the 
concentration  at  the  crack  mouth,  it  i3  convenient  to  express 
the  result  as 

x 

J;  ~  4  %(pWp)  +  \  ]  (C-7) 

Op 

where  the  crack  tip  is  located  at  x  equal  to  p  and  the  prime 
denotes  differentiation  with  respect  to  the  argument  of  a 
function. 

As  indicated  in  the  text,  there  are  two  important  limits 
to  the  behavior  of  an  electrolyte.  The  first  is  when  the 
total  normal  current  density  is  small.  The  ions  then  behave 
as  uncharged  species  obeying  the  equation  of  transport  by 
simple  diffusion.  The  second  is  when  the  current  density  is 
large.  The  concentrations  of  the  two  principal  ions  closely 
approximate  those  of  the  ions  in  a  binary  electrolyte.  The 
potential  gradient  is  determined  by  the  need  to  satisfy  both 
the  boundary  conditions  of  the  principal  ions  and  the  elec¬ 
troneutrality  equation.  Th?  concentrations  of  the  minor 
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ionic  species  are  strongly  affected  by  the  potential  gradient 
but  they  themselves  have  virtually  no  effect  on  it. 

Consider  the  case  of  steady-state  transport  by  simple 
diffusion.  Prom  equations  (C-4)  and  (C-7)  with  (j/  and  V 
equal  to  zero 

Ct  - 
or 


X 


where 

X 

iw  =  =-•  5^  (c-9) 

JP 


The  solution  is  almost  as  simple  for  the  case  of  steady- 
state  transport  in  a  binary  e.  ctrolyte.  Prom  equations 
(C-4)  and  (C-7)  with  V  equal  to  zero 


ZiCiCt)7  +  Cl  —  -  (C-tOa) 

ZA<t>'  +  C[  =  -•  -^(x)  (C-lOb) 


where  f^  end  f^  are  defined  according  to  equation  (C-9). 
Adding  equations  (C-10)  and  using  the  electroneutrality 
equation  to  eliminate  both  the  electrostatic  potential  and 
C^,  it  follows  that 


or 
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=  C? 


(c-11) 


The  concentration  Cg  may  be  obtained  from  equation  (C-11) 
by  simply  interchanging  the  subscripts  one  and  two. 


When  the  normal  flux  of  one  of  the  two  species,  say 
species  two,  is  everywhere  zero,  it  follows  from  equation 
(C-9)  that  fg  is  also  zero.  Since  the  ith  diffusion  co¬ 
efficient  enters  the  result  only  through  f^  it  is  clear 
that  the  solution  docs  not  depend  on  the  diffusion  coeffi¬ 
cient  of  a  species  with  zero  flux.  Comparing  equations 
(C-S)  and  (C-11)  with  fg  equal  to  zero,  it  may  be  seen  that 


[ct  -  cS 


1  -  OH  KM 
i  -  Cfwr 


(C-1 2) 


For  the  problem  presented  in  the  text  this  ratio  was  one- 
half.  It  will  typically  lie  between  one-quarter  and  three- 
quarters  . 

The  electrostatic  potential  may  be  obtained  from  equa¬ 
tions  (C-lOa)  and  (C-11).  However,  we  will  use  a  different 
method  in  order  to  show  how  an  integral  representation  may 
be  used  to  solve  certain  electrochemical  transport  problems. 

Consider  a  steady  state  system  for  which  the  solvent 
velocity  is  zero.  From  equations  (C-4,  -7,  and  -9),  the 
differential  equation  governing  ion  transport  may  be  written 
in  the  form 


r  «/ 

[Gt-2Xp(Z;$)j  —  '  $tCXp(Zi<|>) 


(C-13) 


Integrating, 
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a 

Ci  =-•  1C?  +  \  Sc^ixptZ^C^de^xpt-Z^tx)]  (c-14) 

°x 

where,  for  convenience,  (J)(a)  has  been  set  equal  to  zero. 

Note  that,  as  before,  the  kernel  hi  contains  all  the  in¬ 
formation  regarding  the  normal  flux  and  the  crack  shape. 

When  the  potential,  0,  is  everywhere  zero,  equation  (C-8) 
is  recovered.  When  (£)  is  known  or  otherwise  specified,  as 
for  a  minor  ionic  species  in  a  near-binary  electrolyte, 
equation  (C— 1 4)  expresses  the  solution  and  no  additional 
information  is  required. 

In  general,  the  electrostatic  potential  is  not  known 
and  the  electroneutrality  equation  must  be  used.  Thus,  from 
equations  (C-5)  and  (C-14), 

-  o  (c-15) 

<  v  X 


Every  steady-state  one-dimensional  electrochemical 
transport  problem  is  reducible  to  the  problem  of  finding  a 
value  of  (J)  that  satisfies  equation  (C-15).  In  practice, 
there  are  two  ways  to  do  this.  The  first  is  direct  numeri¬ 
cal  solution  of  the  integral  equation.  The  second  is, 
through  repeated  differentiation,  to  derive  an  equivalent 
differential  equation.  It  may  be  possible  to  solve  this 
equation  analytically.  Otherwise,  it  can  be  solved  numeri¬ 
cally. 

The  derivation  of  an  equivalent  differential  equation 
is  straightforward  but  laborious.  For  example,  in  the  case 
of  a  binary  electrolyte,  equation  (C-15)  may  be  written 


U 


.a 

x  w  * 


( C— 1 6 ) 
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where 


cxp[(j);x)] 


(C— 1 7 ) 


Differentiating  once  and  rearranging, 


^  1 1^  ^ 

-fe.+  ^yy,  1  +  22«i-20C?  *  yvy'^eo'.w^  -  o 

U  X 


(C-1 8) 


Differentiating  again  and  rearranging  yields  the  desired 
differential  equation. 


[(2i^+z2y^/]  +  tit.- v  :l  j  y  ~  o 


(C-1 9) 


Integration  of  equation  (C-1 9)  is  straightforward. 
After  evaluating  the  constants  of  integration  using  equa¬ 
tions  (C-17)  and  (C-18),  it  follows  that 


§  —  !n: 


{  . . lIiSiCC  ' 


(C-20) 


The  correctness  of  this  result  may  be  verified  by  com¬ 
parison  with  the  result  obtained  directly  from  equations 
(C-lOa)  and  (C-11).  The  labor  involved  does  not  justify 
using  equation  (C-1 5)  for  a  binary  system.  For  more  com¬ 
plex  electrolytes,  however,  there  is  no  other  straightforward 
method  for  solving  the  steady-state  transport  equations,  and 
a  method  based  on  the  solution  of  equation  (C-1 5)  is  neces¬ 
sary  whether  it  is  recognizable  as  such  or  not. 
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APPENDIX  D 

Electrochemical  Transport  in  a  1-1-1  Ternary  Electrolyte 

Determination  of  the  implicit  function  (])  defined  hy 
equation  ( C— 1 5 )  in  Appendix  C  is  the  key  step  in  solving 
one-dimensional  steady-state  electrochemical  transport  prob¬ 
lems.  Since  the  equation  is  non-linear,  it  usually  is  im¬ 
possible  to  obtain  a  solution  in  closed-form  and  approximate 
methods  must  be  used.  In  several  cases,  however,  closed- 
form  solutions  can  be  displayed  and  one  of  these  is  presented 
in  this  appendix. 

Consider  a  1-1-1  ternary  electrolyte  for  which  the  nor¬ 
mal  flux  of  one  species  is  everywhere  zero.  In  the  termi¬ 
nology  of  Appendix  C 

Zi  =  '2*  *■=  Z;  +i  (D-la) 


=  0 


(D-lb) 


With  these  values,  equation  ( C— i 5 )  can  be  written  in 


the  form 


a 

f  t  (f\  f 

C!  ~  c|u*  -  uf- \ -'uco'-H  +  c|  +  = 


(D-2) 


where  u(x)  is  defined  according  to  equation  (C-17) 
entiating  once  and  rearranging 


Differ- 


r  ,  1  f 

IV  Si]  ux  -  2C*S  -  z\  U(Od*  rv 


0 


(D-3) 


Differentiating  and  rearranging  again,  the  desired  differ¬ 
ential  equation  is  obtained. 


[«.-y 

=  0 

(D-4) 

By  the  definition  of  u(x) 

and  equation  (D-3), 

U(a)  = 

exp[(j)(a)]  =  i 

(D-5a) 

U(a)  « 

ac*/(i«0-yal) 

(D-5b) 

Therefore,  integrating  equation  (D-4)  twice  and  using  equa¬ 
tions  (D-5),  it  follows  that 


[  to  -  :t(y,]de, _ 

M  *  M&i\ 


i!  zc ;  ■'  "(is 


(D-6) 


The  ion  concentrations  may  be  determined  from  equation 
(D-6)  and  equation  (C-H).  For  example,  for  species  two 


(D-7) 


To  carry  out  the  integration,  it  is  convenient  to  define  a 
new  variable.  Note  that  equation  (D-6)  may  be  written  in 
the  form 

a 
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or 


u. 

<t>  ln(i  (  +  y^]<m  ~  '-Via 


where 


HV,  = 


f  _  _____  $*<0<U 

\  +  4Uuru 


H +  vrtd?] 


(D-8) 


(D-9) 


Using  equations  (D-8)  and  (D-9),  equation  (D-7)  becomes 

a 

C2  =  C«{i  -  Vxp^V]  ^?d£5l^p[-?]axp[4)  +  4>] 

or 

,  ra 

c2  =  c?1i  +  (D-io) 

1  x 

Since  f^(x)  was  assumed  to  be  zero, 

Ct  =  ci?Ap[-'t)(^3  (D'1 1 ) 

and  by  electroneutrality 

C,  =  Cl"  Cl  (D-12) 

Equations  (D-6,  -10,  -11,  and  -12)  express  the  general 
solution  for  a  1-1-1  ternary  electrolyte  for  which  one  of  the 
positively  charged  species  has  zero  normal  flux.  Note  that 
for  a  1-1-1  ternary  containing  two  negatively  charged  ions 


and  one  positively  charged  ion  and  for  which  one  of  the 
negative  ions  has  zero  normal  fiux,  the  same  solution  is 
obtained  except  that  (J)  is  the  negative  of  the  value  given 
by  equation  (D-6).  Also,  the  signs  of  the  charge  numbers 
must  then  be  taken  opposite  to  those  shown  in  equation 
(D-la) . 
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Fig.  2  Coordinate  transformations  from  the  physical 
(x-y)  plane  onto  the  descriptive  (‘  -v)  plane 
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Fig.  5  Longitudinal  concentration  profiles  in 
a  hyperbolic  crack  a z  three  different 
times  (transport  by  simple  diffusion). 


Fig.  7  Variation  of  the  crack-tip  pH  with  time  for  three  dif¬ 
ferent  crack  shapes  (transport  by  simple  diffusion). 
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Fig.  9  Longitudinal  concentration  profiles  in 
a  pie-shaped  crack  at  three  different 
times  (transport,  ty  simple  diffusion). 


Fig. 


Change  with  time  in  the  ratio  of  the  longitudinal  to  the 
transverse  concentration  gradient  at  three  points  in  a 
rectangular  crack  (transport  by  simple  diffusion). 
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Fig.  15  Dependence  on  crack  length  of  the  steady-state 

£ri?k*"tlp  1  pH  for  F  .OS-centimeier-Viid/  byper- 
b--.jc  cracks  (transport  by  simple  diffusion). 
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Variation  of  the  normal  flux  (current  density)  along 
the  crack  wall  with  concentration-dependent  boundary 
conditions  (transport  by  simple  diffusion). 
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Fig.  21  Variation  of  the  crack-tip  pH  and  oxide-film  thickness 
with  time  for  a  hyperbolic  crack  with  charge-transfer 
boundary  conditions  (transport  by  simple  diffusion). 


Fig.  23  Convection  in  a  pie-shaped  crack. 
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Fig.  24  Double  edge-notched  plate  tension  specimen  used 
to  investigate  the  effects  of  solvent  convection. 
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Fig.  27  '  Steady-state  '  lcngi  tudi r.ai  concentrate  Jn 

profile  in  a  hj'perbolic  crack  showing 
the  influence  of  convective  transport. 
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Fig.  28  Dependence  on  crack  width  cf  the  steady-state 
crack-iip  concentration  for  C. 20-centimeter-long 
hyperbolic  cracks  (transport  by  simple  diffusion). 
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Fig.  31  Variation  cf  the  ion  concentrations  with  time  in  a  rec¬ 
tangular  crack  (transport  by  diffusion  and  migration). 
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Fig-  35  Steady-state  longitudinal  concentration  profiles 
in  a  rectangular  crack  for  several  values  of  the 
flux  parameter  /\  (transport  by  simple  diffusion) 
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Fig.  37  Dependence  of  the  concentration  and  flux 
in  a  rectangular  orach  on  the  flux  param¬ 
eter  (transport  by  simple  diffusion). 
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Fig.  38  Variation  of  the  steady-state  crack-tip  pH  with 
normal  current  density  for  a  rectangular  crack 
shewing  the  effect  of  transport  by  migration. 
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Fig.  4C  Variation  with  normal  current  density  of 
the  potential  drop  for  a  ternary  electro¬ 
lyte  contained  in  a  rectangular  crack. 


